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§ 1. [yTRODUCTION 


THE present review article aims at bringing out essential features of the 
development of the theory of antiferromagnetism together with some 
of the experimental facts hitherto obtained. In 1950, at the International 
Conference on Ferromagnetism and Antiferromagnetism at Grenoble, 
van Vleck gave an excellent review of the theoretical side of the same 
subject and Bizette of the experimental side and these have been published 
in Le Journal de Physique et le Radium (1951). Since then, however, a 
considerable advance has been made in this field, notably in the study 
of the resonance absorption of microwaves, the problem of anisotropy 
energy, spin-wave theory, structure determination by means of neutron 
diffraction and of proton resonance absorption, and so on. The present 
article inevitably overlaps to some extent with those of van Vleck and 
Bizette, but is concerned more with these recent advances. The same 
field is also covered, though in less detail, by Lidiard (1954). 

The concept of antiferromagnetism was first put forward by Néel 
(1932 a) in connection with his study of the paramagnetic susceptibility 
of metals and alloys of transition elements. He noticed that metals 
such as Pt, Pd, Mn, Cr or alloys such as Pt-rich Pt-Co and Pd-rich 
Pd-Ni show an almost temperature-independent susceptibility which 
is too large to be explicable in terms of Pauli’s free electron paramagne- 
tism. He supposed that the atoms of these substances possess magnetic 
moments and that a negative exchange coupling is operating among 
them. As a consequence, he guessed that the arrangement of the atomic 
magnetic moments in the crystal lattice is such that they point upward 
and downward in an alternating manner when the temperature is suffi- 
ciently low. He showed through a simple calculation that the suscepti- 
bility of such a system is nearly constant at low temperatures and obeys 
at high temperatures the Curie-Weiss law 


ee 
X= 76° .. 33°52) en 


with a negative paramagnetic Curie temperature—@. He also (Néel 
1932 b) made measurements with Mn and Cr diluted in noble metals and 
found that these atoms actually possess a magnetic moment of 4 to 5 
/y (Hp is a Bohr magneton) and suggested for pure Mn and Cr values of 
@ higher than 1000°K. 

Néel (1936 a and b) extended his experiments to other systems of 
alloys and worked out his theory more in detail. In these works the 
following points may be noticed : (1) He studied the magnetic properties 
of binary alloys of two transition elements, using the Weiss molecular 
field approximation. This problem was studied by himself in greater 
detail in connection with ferrites (Néel 1948). (2) A formula for the 
antiferromagnetic Curie point was derived with the same approximation. 
(3) By plotting the coefficient of the Weiss field deduced from experiment 
against d—6, namely the difference between the atomic distance d in the 
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metals, or the alloys, and the mean diameter of the d-shells of the 
constituent atoms 5, he found that the points fall on a smooth curve 
whose positive part corresponded to ferromagnets and the negative 
part to antiferromagnets. (4) The anisotropy energy in antiferromagnetic 
substances was suggested to be of the same order of magnitude as in 
ferromagnetic ones, and the susceptibility at absolute zero was computed 
for the case of uniaxial anisotropy as a function of the intensity and 
direction of the applied magnetic field. An interesting result of this 
analysis is the existence of a certain critical field strength, H,, for the 
case where the field is applied along the axis of easy magnetization. 
The antiparallel magnetizations of the two sublattices turn from the 
direction of the easy axis to that perpendicular to it when the field 
exceeds this critical value. This topic will be dealt with in detail in § 5 
of the present article, but a preliminary representation of it might 
conveniently be given here as follows :— 

Suppose we have an anisotropy energy (per unit volume) of the form 
K sin? ¢, where K is a positive constant and ¢ the angle between the easy 
axis and the direction of the magnetization. We shall use the notation 
A for the latter. If the external field is parallel to this direction J, the 
induced magnetic moment of the system will also be parallel to it. We 
shall denote the corresponding susceptibility as y,. If the magnetic 
field is perpendicular, the induced moment will likewise be perpendicular ; 
the corresponding susceptibility will be denoted as y,. It can be shown 
that y, is independent of temperature to a good approximation, that 
x decreases with decreasing temperature, vanishing at absolute zero, 
and that y,=x, at the antiferromagnetic Curie temperature. Now 
suppose an external magnetic field H is applied along the easy axis. 
Initially the direction 4 coincides with this axis, and the change in free 
energy per unit volume of the system is given by —y,H?/2. With in- 
creasing field strength, however, this change in free energy becomes 
equal to K—y,H?/2, that is, the free energy for 4 perpendicular to the 
easy axis, and then exceeds it, since in general y,>y,. The critical 
field is obtained by solving the equation 


—4y,H?=K —}ty,H? 
or 
A= {2K/(x1— Xu) }?. a hi.) Iestactiag, Beaten re) 


This remarkable phenomenon of the ‘flopping’ of the direction 4 
was actually observed by Gorter and his co-workers (1951-53) in single 
erystals of CuCl, .2H,O. This substance has a Néel point—another name 
for the antiferromagnetic Curie point initiated by these authors—of 
4-3°K and is orthorhombic. In the orthorhombic case, the magnetization 
vectors turn from the direction of the easy axis to that of the inter- 
mediate axis. Up to the present, however, in no other substance has 
this phenomenon been observed, presumably because in most substances 
the value of the critical field strength is too high. 
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In the early days Néel employed the expression paramagneétisme 
constant for antiferromagnetism. The latter term can be traced back to 
Hulthén’s work (1936, 1938) in which the first spin-wave treatment of 
antiferromagnetism was given. (Prior to that, an approximate quantum- 
mechanical solution of antiferromagnetism for the case of one dimension 
had been given by Slater (1930) and the exact solution of it by Bethe 
(1931); the latter treats one-dimensional ferromagnetism, but the sign 
of the exchange integral can be reversed without modifying the formalism. ) 
Considerable advances have since been made in the spin-wave theory by 
Anderson (1951, 1952), Ziman (1952-53), Kubo (1952, 1953), Nakamura 
(1952), Tessman (1952), Keffer, Kaplan and Yafet (1953) and Kasteleijn 
(1952). On the other hand, Néel’s original theory, based on the Weiss 
approximation, has been extended by Bitter (1937), van Vleck (1941), Néel 
himself (1948), Anderson (1950 b), Garrett (1951), Nagamiya (1951-54), 
Yosida (1951-53), Kittel (1951), Keffer and Kittel (1952), Gorter and 
Haantijes (1952), Ubbink (1953 a, b), Smart (1953), Wangsness (1952-53), 
and possibly by others. Among these authors, Néel and Anderson, as 
well as Smart, recognized the importance of the next-to-nearest neighbours 
for the Néel point and the susceptibility, Kittel, Nagamiya, Yosida, 
Keffer and Kittel, Gorter and Haantjes, Ubbink and Wangsness treated 
the problems connected with the anisotropy energy and the resonance 
absorption of microwaves by antiferromagnetic substances, and Bitter, 
Garrett and Leiden workers under Gorter considered, among others, 
the case in which the applied magnetic field is of the order of, or greater 
than, the Weiss molecular field. At present, only a few substances are 
known in which the last mentioned case is really important. 

A number of papers also appeared which tried to improve the statistical- 
mechanical treatment of the problem. Li (1951) applied the Bethe— 
Peierls—Weiss method to antiferromagnetism, Nakamura (1953) applied 
another similar method, and Kubo and his co-workers (1951—52) and also 
Oguchi and Obata (1952) developed a moment expansion method and 
other expansion methods in antiferromagnetism. : 

The resonance absorption of a microwave by antiferromagnets was 
first attacked experimentally by Trounson, Bleil, Wangsness and Maxwell 
(1950) with a polyerystal of Cr,O,. They observed no remarkable effect 
below the Néel point, because the resonance frequency in this substance 
should lie in the sub-millimeter region according to theory. Ubbink 
(1951-53) at Leiden, however, discovered peculiar resonance phenomena 
in the single crystal of CuCl, . 2H,O, which have been explained by his 
own theory and that of Nagamiya and Yosida. 

A fundamental problem of antiferromagnetism is the origin of the 
exchange interaction. Bizette (1946) and Néel (1948) recognized for 
antiferromagnetic compounds the importance of superexchange, that is, 
the exchange coupling between magnetic atoms through the intermediary 
of a negative non-magnetic ion. This sort of interaction was first con- 
sidered by Kramers (1934) and its realistic treatment was given by 
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Anderson (1950a). For metallic antiferromagnets a band-theoretical 
treatment might be useful, and attempts of this kind have been made 
by Matsubara (1953) and Lidiard (1953) following a suggestion made by 
Slater (1951), who with Pratt (1953), also considered another mechanism 
of superexchange for compounds. 

Another fundamental problem of antiferromagnetism is the origin of 
the anisotropy energy, which was studied by Keffer (1952), Yosida (1952), 
Moriya and Yosida (1953), Niira and Oguchi (1954) and Nakamura and 
Taketa (1954) for a few substances. 

So far we have cited mainly theoretical works, but there are great 
many valuable experimental studies. It is not convenient to cite them 
in this introduction, but we shall summarize them in the next two sections 
and refer to them where appropriate elsewhere. 
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Part I. FUNDAMENTAL Facts AND ELEMENTARY THEORY 


§ 2, ANTIFERROMAGNETIC ARRANGEMENT AS REVEALED BY 
Neutron DirrRactioN AND OTHER EXPERIMENTS 


It may be useful first to present some experimental results on the 
magnetic structure of a number of antiferromagnetic substances. 

The most direct method of determining the arrangement of spins in an 
antiferromagnetic crystal is that of neutron diffraction. An ordinary 
Debye-Scherrer pattern is obtained by passing a monochromatized 
neutron beam through a powdered sample, the neutrons being scattered 
by the nuclei and by the magnetic moments of the atoms. From the 
latter the direction of the magnetization can be definitely determined in 
most cases, but the superstructure arrangement of positive and negative 
spins seems to be less definitely determined since one has to depend on the 
method of trial and error. X-ray analysis also gives an indication of the 
superstructure, though indirectly. In the case of CuCl, . 2H,O, proton 
magnetic resonance was useful for the superstructure determination. 

MnO, FeO, CoO, NiO, MnS, MnSe are all of the NaCl-type. Shull 
and his co-workers (1949, 1951 b) found for them by neutron diffraction 
that the magnetic atoms on each (111) plane have the same direction of 
spins and that the direction alternates as one goes from one (111) plane to 
the adjacent one. In MnO, MnS, MnSe, CoO and NiO, the direction of 
spins was found to be one of the cubic principal axes, so that the symmetry 
of the spin arrangement is monoclinic. In FeO, however, the spin 
direction was found to be [111], the symmetry thus being rhombohedral. 
X-ray analysis of these compounds, except MnSe, has been carried out by 
Rooksby (1943, 1948), Tombs and Rooksby (1950, 1951), Shimomura and 
Nishiyama (1948, 1949), Shimomura (1953), Greenwald and Smart (1950, 
' 1951) and Greenwald (1953). They found that all these substances, 
except CoO, deform slightly below the Néel point and become rhombo- 
hedral, a change (contraction except in FeO) of a fraction of 1°, occurring 
along the [111] direction. The atoms with opposite spins thus appear to 
attract each other in most substances, while magnetostriction appears 
to be of minor importance. In CoO, the lattice deforms tetragonally, a 
contraction occurring along one of the original cubic axes. There is at 
present no reasonable explanation for this if one assumes the evidence of 
neutron diffraction to be conclusive. Shimomura (1953) dissolved copper 
oxide in NiO and observed the change in lattice constant as a function 
of the copper concentration ; above 21-5 atomic per cent Cu, the lattice 
was found to be perfectly cubic, but at 24°, he still found antiferro- 
magnetism. It may, however, be premature to conclude that in the 
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latter case the exchange force is inactive in deforming the lattice, since we 
do not know the magnetic structure of this solid solution. 

Greenwald (1953) observed also in hexagonal antiferromagnetic 
compounds, Cr,0, and MnTe, a contraction along the principal axis, 
which increased as the temperature was lowered from the Néel point. 

MnO, is another antiferromagnetic substance. It is tetragonal body- 
centred, and the suggested magnetic structure on the basis of neutron 
diffraction (Krickson 1952) is such that the corner atoms and the body- 
centred atoms form separate antiferromagnetic superstructures, the 
directions of the magnetizations of the two being normal to the c-axis and 
at right angles to each other, possibly one pointing in [110] and the 
other in [110]. 


Fig. 1 


MAGNETIC Renew t 


UNIT CELL UNIT CELL 


Mn ATOMS IN MnO 
Antiferromagnetic structure existing in MnO below its Néel temperature of 
120°x. The magnetic unit cell has twice the linear dimensions of the 
chemical unit cell. Only Mn ions are shown in the diagram. (After 
Shull, Strauser and Wollan 1951 b.) 


The magnetic structure of CrSb (NiAs-type) was investigated by Snow 
(1952) and the spins were found to point parallel to the c-axis, the atoms in 
the same (001) plane pointing in the same direction and the atoms in the 

adjacent (001) plane in the opposite one. 

MnF,, FeF,, CoF, and NiF, are all of the rutile type. Except in Nif,, 
the metallic atoms at the corners of the tetragonal unit cell point down- 
wards and those at the body-centres point upwards (Erickson and Shull 
1951, Erickson 1953). In NiF, the spin direction does not appear to 
coincide exactly with the c-axis but is inclined to it at an angle of 10°, 
possibly pointing to one of the neighbouring fluorine ions. Erickson 
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determined the magnitude of the spontaneous magnetization of each 
sublattice as a function of temperature from the observed intensity of the 
100 line. It was well represented by the Brillouin function (see § 4), with 
S=5/2, 2, 3/2 and 1, respectively, for the fluorides of Mn, Fe, Co and Ni. 
Figures 2 and 3 show his results. At the same time he found the Landé 
factors to be 1-99, 2:32, 2-00 (2:03 from the diffuse scattering above the 
Néel point) and 2-04, respectively. 

The oxide, x-hematite, Fe,O,, is interesting in that it shows a feeble 
ferromagnetism, possibly due to the presence in it of parasitic ferro- 
magnetic substances. The main structure is without doubt antiferro- 
magnetic, and its structure was determined by Shull, Strauser and Wollan 
(1951 b) as shown in fig. 4. The figure also shows the structure found in 
Cr,0, (Brockhouse 1953). Below —20°c, the direction of the spins is 
parallel to the crystallographic principal axis, but perpendicular above. 


Fig. 2 
MnF, (100) FeF, (100) 


2 2 -24 2 
Fiog = 5435 Bgyaly) #10 “CM 


FZ = 4.55 Baty) #107? cM? 


T, =75°K T, = 90°K 


20 40 60 80 100 20 40 60 80 100 
TEMPERATURE, °K 


Temperature dependence of the 100 peak intensity of the diffracted neutron 
beam, in terms of the crystal structure factor, for MnF, and FeF,. The 


points are fitted with an appropriate Brillouin saturation curve. (After 
Brickson 1953.) 7’, stands for our Ty. 


At this temperature an anomaly of the susceptibility has been observed, 
and will be referred to later (§ 3). The neutron diffraction measurements 
were done between 80°K and 1000°k, that is, up to a temperature a little 
over the Néel point. It was found that the intensities of the diffraction 
lines, which should be proportional to the square of the spontaneous 
magnetization, decrease linearly with temperature over a wide range from 
250°K to 1000°K. More detailed measurements in the neighbourhood of 
— 20°C have been made by Corliss, Hastings and Goldman (1954). 
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Some of the ferrites, such as ZnOFe,O3, seem to be antiferromagnetic. 
The neutron diffraction experiment reveals not only whether the ferrite 
is normal or inverse but also the direction of the spontaneous magnetiz- 
ation of the sublattices. In all cases hitherto investigated (ferrites of Fe, 


Fig. 3 
CoF, (100) NiF, (100) 


24 2 2 - i 
cm? Figg = 070 8, ty) x 10724 om 


alps 2 = 
Figo = 1.90 BS ,oly) x 10 


T, = 50°k T, =83°K 


* 10 20 30 40 50 20 40 60 80 100 
TEMPERATURE, °K 


Temperature dependence of the (100) crystal structure factors in CoF’, and NiF,. 
The saturation curves indicate Néel temperature of 50°K for CoF, 
and 83°K for NiF,. (After Erickson 1953.) TT’, stands for our T'y. 


Fig. 4 


Magnetic structure of o-Fe,0, at room temperature. The thombohedral unit 
cell containing four unique iron atoms is shown. (After Shull, Strauser 


and Wollan 1951 b.) 


Ni, Mg and Zn), except in ZnOFe,O;, the spins on one of the two sub- 
lattices were found to point in the opposite directions to those on the 
other sublattice, as predicted by Néel (1948), and the preferred direction 
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to be one of the cubic axes (Shull, Wollan and Strauser 1951 a, Hastings 
and Corliss 1953, Corliss, Hastings and Brockman 1953). In ZnOFe,0; 
no coherent magnetic scattering was observed at room temperature, 
possibly because the Néel temperature is at liquid air temperatures. 

A number of transition metals have also been investigated by Shull and 
Wilkinson (1953). In V, Nb, Mo and W, no measureable magnetic 
scattering lines, neither diffuse nor sharp lines, were observed. In Cr 
and «-manganese, however, an antiferromagnetic arrangement was 
observed, which varied as a Brillouin function, and the suggested Néel 
points are 473°K and 100°K respectively. The atomic magnetic moment 
deduced from the intensity of lines for aligned moments in Cr and that 
deduced from the diffuse lines above the Néel point in Mn are 0-40 4, and 
about’ 0-5 uz respectively, not 4 and 5, as one would expect. Néel 
(1932 b) obtained the latter values from magnetic measurements on these 
metals dissolved in noble metals. A test was done for Cr at another 
neutron velocity 25°, higher, as there was a suspicion that a rapidly 
changing polarity of the moment might have affected the observed value of 
the moment. No effect was found, however. Néel (1953 b) points out a 
possibility that the observed moment is concerned with the moment of 
those Cr atoms which are situated in the regions where impurity atoms 
stabilize the antiferromagnetic sublattices in such a way as to fix their 
polarities during the passage time of neutrons. A band-theoretical 
treatment of antiferromagnetism, or a collective electron antiferro- 
magnetism, worked out by Matsubara (1953) and also by Lidiard (1953) 
predicts a small effective moment and thus agrees with experiment, but 
Anderson (1953) remarks that a theory of this type is inappropriate in the 
case where a form-factor scattering is observed above the Néel point, as in 
the case of Mn. His criticism may be valid if the theory is understood in 
its naive sense. The low values of the Néel point are also contradictory 
to those suggested by Néel from magnetic measurements. We shall 
return to these topics in the next section. 

So far we have summarized the results obtained by neutron and x-ray 
diffraction experiments. There is a proton resonance experiment in the 
case of CuCl, . 2H,0 (Poulis and Hardeman 1952 a, b). This substance 
has a Néel point of 4-3°K and is most interesting from various points of 
view, as we shall see in later sections. The orderéd arrangement of the 
cupric ions will produce different internal magnetic fields at different 
proton positions and thus will cause a splitting of the proton resonance 
line. The magnetic unit cell obtained by Poulis is shown in fig. 5 (a) 
which is twice as large as the chemical unit cell along the c-axis (fig. 5 (a) 
and (b)). The spins point parallel to the a-axis if there is no external 
field, but they turn to the b-axis when a field greater than H,, the critical 
field strength, is applied along the a-axis. Figure 6 shows the magneti- 
zation as a function of temperature, deduced from proton and microwave 
resonance experiments, from specific heat measurement, and that caleu- 
lated on the basis of Weiss molecular field approximation, 
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(0) (c) 
The crystal structure of CuCl, .2H,O. (a) Projection on (001), (b) Projection 
on (010). Large open circles represent Cl atoms, small open circles O 
atoms and black circles Cu atoms. (After Harker 1936.) 


Fig. 6 


QYOOE 


The square of the magnetization per sublattice as a function of temperature for 
CuCl, .2H,O. Curve 1: from the proton resonance measurements of 
Poulis, curve 2: from the molecular field theory with S=1/2, curve 3: 
from the specific heat measurements of Friedberg, curve 4; from H, and 
a; see§6, (After Ubbink 1953 a.) 
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§ 3. Suscepripitity, Speciric HEAT AND OTHER PROPERTIES 


Of the many experiments on susceptibility and other properties only a 
few typical ones will be quoted here. 

The earliest measurements of the susceptibilities characteristic of anti- 
ferromagnetic materials are, as far as the writers could find, those done in 
Japan by Honda and his co-workers (Honda and Soné 1914, Ishiwara 
1914, Honda and Ishiwara 1915), together with Weiss and Foéx’s (1911) 
measurements of the susceptibility of y-Fe which exists in a limited range 
of temperature with no Néel point in that range. They investigated MnO, 
Cr,03, CuBr,, CuO and NiO between liquid air temperature and 1000°c, 


Fig. 7 


Magnetic susceptibility of MnO. White circles are for H=5000 oe and black 
circles are for H=24000 oe. (After Bizette, Squire and Tsai 1938.) 


and later (Shimizu 1930) Mn between room temperature and 1250°c, <A 
well-defined peak in the susceptibility-temperature curve was found in 
MnO at 108°K, a fairly well-defined peak in Cr,O, in the neighbourhood of 
40°o and a broad peak in other substances. Bizette, Squire and Tsai 
(1938) extended the range of temperature to 14°K and worked possibly 
with a better sample of MnO, because they obtained lower susceptibility 
values. The Néel point was found to be 116°K and the susceptibilit 

above the Néel point was found to obey the Curie-Weiss law BES HE with 
@=-610°K and ¢ j= 4-40 deg. cm*/mole, or an effective atomic moment of 
5:95 pp (theoretically, 5-924), for spin 5/2). Figure 7 shows their result 
Similar results were obtained for Cr,0, by Foéx and Graff (1939) Foax 
and Wucher (1951) and McGuire (quoted in the paper of Trounson et an 


Antiferromagnetism 13 


(1950)). There are many other substances which show similar properties ; 
we summarize the data in table 1. 

Antiferromagnetic substances usually show a A-type anomaly at the 
Néel point in their specific heat curves. Figure 8 shows an example 
obtained by Millar (1928) for manganese oxides. The thermal expansion 


Heat capacity, calories per mole. 
y ) 
5} 


80 100 120 140 160 180 200 -220 240 260 280 300 
Temperature degrees Kelvin. 


The molal heat capacities of MnO, Mn,O, and MnO,. (After Millar 1928. 


Fig. 9 
rx 10° 
70 
a MnO Coefficient de dilatation 
d'aprés M. Foex 
50 


O70 -100 0 +100 +200 +300 +400 °C 
Thermal expansion coefficients of MnO, FeO, CoO and NiO. (After Foéx 1948.) 


coefficient behaves in a similar way ; fig. 9 is due to Foéx (1948). Street 
and Lewis (1951) and Fine (1953) observed a sudden increase in Young’s 
modulus above the Néel point for NiO and CoO, the temperature coefficient 
varying somewhat like that of the thermal expansion, 
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Table 1 


Substance 


9 MnO, 
10 a-Fe,O, 


11 Crs 
12 MnS 


13 FeS 
14 MnSe 
15 MnTe 


16 FeTe 
17 MnF, 


18 FeF, 
19 CoF, 
20 NiF, 


21 VCl, 
22 VCl, 
23 CrCl, 
24 FeCl, 
25 FeCl, 
26 CoCl, 


(1) 


crystal Lattice constant a, b, ¢ ; « 
type 

NaCl* 4-4345 

str. trans. to trig. for < —100°c 
NaCl* 4-332 

str. trans. to trig. «<60° for << —70°C 
NaCl* | 4-2495 (20°c) 

tetr. 4:2552, 4:2058 (—180°c) 
NaCl* 4-1684 

trig. 4-16768 ; 90° 3-8’ (9°c) 
monocl.*| 4-653, 3-410, 5-108 ; 99° 29’ 
Cr;0, 5-43 ; 53° 53’ 
rutile 4-54, 2-88 u=? 
trig. 5:38 ; 54° 50’ 

str. trans. at 307 ~318°K 
rutile* 4:44, 2:89 w=? 
Cr,0,* | 5:4135; 55° 17’ 
NiAs 3-448, 5-754. 
NaCl* 5:212 

str. trans. to trig. for <f, 
NiAs 3-453, 5-670 
NaCl* 5-448 
NiAs 4-124, 6-698 sharp contraction of ¢ 

below 329°K 

NiAs 3-800, 5-651 
rutile* 4-8734, 3:3103 w=—0-310 
rutile* 4-670, 3:297 w=0-31 
rutile* 4-69, 3:19 u=0°31 
rutile* 4-6505, 3-0837 w=0:31 
CdCl, 6-20; 33° 36’ w=? 
trig. 6-69 ; 52° 30’ 
CdCl, 6-16; 33° 26’ 1=0-25 


(2) 


Néel temp. 7'y °K 


x 


122 


198* 
186? 


293 


492 ~647 


~453 
173 


343" 
335~355? 
311° 


323>:¢ 
3504 
84 
~950 


165* 


147 
613 


323%¢ 
3075 


728 
70» 
719 


30 
40 
24 


25 


~220 ? 


~306-5° 
~950 


(3) 


sp. heat 
1162 
117-8° 
183* 
188-5" 


289°7 
520 


168-8 


305% 
304-5 
92-1 


140 


66-5 
78:3 
37-7 
73-2 
104-9 
23°50 


24-90 
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ee ee ea 


(4) (5) (6) References and remarks 
ratio of (The numbers in brackets refer to the 
from XT=0 columns, and the numbers which 
a 2) Cno1 to follow refer to the corresponding 
other meas. XT=Ty references.) 
116 D 610 4-40 0-69 | (1) 97 (2, 5, 6) 11 (3) 966 ©96 (4) 34 
s.s. CoO 7, FeO 17, MgO 16, 17 
186 D 5702 | 6-248 0-788 | (1) 97, (2, 5, 6) #15 +°63 (3) 967 >96 (4) 34 
190» 4-60> 0-76" | (5 5) 196<T'<460°K 0973<2 7 = 1320°K 
195° 2-42° s.s. MnO 17, MgO 16, 17 
292 D 280 30546 (1) 97, 38’ (2, 5) 61 (3) 108 (4) 34, 95, 71 
289 E s.s. MnO 7, MgO 17, NiO, CuO 108 
2918 
523 D 0-67 | (1) 86 (2, 6) 61 (3) 112 (4) 34, 95 
507 E 
~230 N 0-60 | (2, 6) 53 see also 60, (3) 57 (4) 110 
170 D (2) 54 (3) 3 (4) 33 
168 C Nee ores in oS eh 31 
720 0-55 (2) 954 ®74-(5) 7 
305-8 D 1070° 256% 0:76? | (1) 38 (2, 5, 6) #100 ©30 °52 432 (3) 94 P56 
| (4) 56, 71 
3188 550» 1-82 0-85” | res. 100, 64 
4934 
0:93 | (2, 6) 18 (3) 59 s.c. x 18 
~950 C 2000 4-4 (2, 3, 4) 26 (5) 69 res. 5 anis. change at 
—10~—30°c accomp. by weak ferro. 23, 
70, 19 
CrS,+, 7>0-14 shows ferro. 46, 48, 106 
5284 430° 02527) (1) (79).(2, 06) SLT. '88 81, (3) 2. (6)% at 
H=32 000 ¢ 
454° 3-485° 0:76" | res. 64 
857 | 3-44 (2, 5) 20 trans. at 120°c 
FeS,4,, e>0-1 shows ferro. 47, 73, 10, 106 
361 4-01 (3) 58 (5) 81 large therm. hysteresis 13, 62 
C 102 
6902 4-59 0-68 | (1) 38’ (2, 5, 6) 281 ©88 °104 (3) 58 C 89, 102 
=>0-5» 
220 0-92 (5) 103 
75 N 113-24. 4-082 0-72 | (1) 39 (2, 5, 6) #12 ©25 ¢22 (3) 91 (4) 27 
97-0°| 4:-47¢ 0-78» | s.s. ZnF, 22 yn—x, 92, 408.c. x 109 
90 N 117 3°88 0-72 | (2, 5, 6) 14 (3) 94 (4) 27 xu—x, 93 
50 N 52-7! 3-29 (3) 94 (4) 27 (5) 17 see also 24 yy—x, 93 
83 N 115-64 1-528% (1) 45 (3) 94 (4) 27.(5) 917 25 
100° 1-3? weak ferro. appears below Tn 111 
565 2-13 (5) 90 
- 30-1] 1-005 0-72 | (2, 5, 6) 90 (3) 84 
149 3:26 0:84 | (2, 5, 6) 90 
—48-0| 3-59 <0:2 | (2, 5, 6) 90 (3) 98 
11-5) 4:07 (5) 90 
—38-1| 3:46 ~0:6 | (2, 5, 6) 90 (3) 99 


P.M. SUPPL.—JAN. 1955 


oy 
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Table 1—cont. 
aS eee 
(1) (2) (3) 
tal Néel temp. 7'y °K 
thks Lattice constant a, b, ¢ ; « _—<§<——<—<— 
Substance type i ; nae 
27 NiCl, CdCl, | 613; 33°36’ w=? 50 | 52-35 
28 CuCl, nae 
29 CuBr, 0) 
30 CrSb NiAs 4-107, 5-468 ~673 
31 MnAs NiAs 3-710, 5-691 (<45°c) 399 | ~45 
3-659, 5-691 (>45°c) lat. heat 
32 MnBi NiAs 4-30, 6-12 621 ~633 
c—6-11 > 5-82 at 350°C 718 
\ (lat. heat) 
ys5.Cr O° |OKOKE. 2:8796 ~1673 
if 34 a-Mn complex 95 
, 35 MnCl, . 4H,O ~1-68 1-622 
36 MnBr, . 41,0 ~2-2 
{| 37 CuCl, . 2H,O rhomb.* | 7-38, 8-04, 3-72 4-31 
oF 38 Co(NH,)» monocl. | 9:23, 12:49, 6-23; 106° 56’ 0-084 
(SONA ae) 
39 FeCO, NaNO, | 5-754; 47° 25’ uw=0-27 


40 FeCO, . 2MgCO, 


Those marked with * in column ‘crystal type’ have known spin super structures ; 


for the detail see text § 2. 


} Phase diagram in H—T-plane for antiferromagnetism = paramagnetism (or = ferro- 


magnetism) is measured. 
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(4) (5) (6) References and remarks 
ratio of (The numbers in brackets refer to the 
from XT=0 columns and the numbers which 
——— () Col to follow refer to the corresponding 
other meas. XT=Ty references.) 
—68-2| 1-36 (2, 5, 6) 90 (3) 21 
109 0-536 0-64 | (2, 5, 6) 90 
(2) 53 see also 60 
723 D ~1000 ~t (2, 5, 6) 30 (4) 87 
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K, —0-050 
35 Op 14 | 3-49 (4, 5) 17 Xy—%, 17 
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Halides of the iron group are also typical antiferromagnetics. They 
show similar properties to the oxides cited above. We will cite three 
experimental facts : 


1. As in the case of MnO (fig. 7), the susceptibility of MnF, has been 
found to depend on the measuring field strength below the Néel point of 
72°x (Bizette and Tsai 1939, de Haas, Schulz and Koolhaas 1940). Theory 
predicts (see § 5) a dependence on H? of the susceptibility, but the measure- 
ments are not accurate enough to decide whether it changes proportionally 
to H or H?. 

2. Stout and Griffel (1949), Griffel and Stout (1950) and Stout and 
Matarrese (1953) observed the parallel and perpendicular susceptibilities 
of MnF,, FeF, and CoF, in single crystals. The results are as seen in 
figs. 10-12. The susceptibility y,, with spins and field parallel to the tetra- 
gonal axis, goes down to zero as the temperature decreases to absolute 
zero and the susceptibility y,, with spins parallel but field perpendicular 
to the tetragonal axis, is nearly constant, though it varies slowly with 
temperature. Experiments give their difference, and the separate values 
were obtained by combining it with the susceptibility of the powdered 
specimens. There is a pronounced anomaly below and above the Néel 
point. Yosida (1951) gave a rough theoretical estimate of the anisotropy, 
Xi—Xy, above the Néel point for MnF, and Niira and Oguchi (1954) 
analysed the data for FeF, to a certain extent (see § 11). Evidently, we 
have here (except in MnF,) to do with the splitting of the electronic energy 
levels of the magnetic ions in the crystalline electric field, and for the 
complete understanding of the susceptibility curves it would be necessary 
to do an elaborate statistical-mechanical treatment. 

3. The specific heat of MnF, measured by Stout and Adams (1942) 
shows a sharp peak at the Néel point as in the case of manganese oxides 
(fig. 13). The change of entropy corresponding to excess specific heat 
over its normal values has been estimated to be 1-2 cal/deg. mole, which is 
close to R In 2=1-377, This fact shows that here we are really concerned 
with an order—disorder phase change between two orientations of the spins. 
In other cases, such as VCl,, the values of the excess entropy reported are 
smaller. 

Recent susceptibility measurements by Bizette and Tsai (1954) with a 
single crystal and a powder specimen of MnF, give the following interesting 
results. The molal perpendicular and parallel susceptibilities and 


powder susceptibility between liquid helium temperature and room 
temperature are : 


T (°K) 4 14; 20-9) 988 ema 77 90 195 293 
10*yx7, 244 24-4 24:4 24-2 24:1 23-5 292 14:5 10-8 
10Oyun = ~— 1 3:2 21-5 Ses eee 

10°yy, 163 16:6 17:3 23:3) (oo 


Xi—Xy was found to agree with the results of Stout and Griffel, while 
x, and x, differ considerably from those given by these authors; in 
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particular, X, Temains nearly constant below 7'y—70°K. The powder 
specimen was such that the powders were fixed in paraffin, and it was 
found that the susceptibility did not depend on applied field between 5000 


The molal magnetic anisotropy of MnF,. yx, and x, are molal susceptibilities 
parallel and perpendicular to the c-axis. The arrow indicates the 
temperature of the heat capacity maximum. (After Griffel and Stout 
1950.) 


_ (0) 
Molal magnetic susceptibilities of MnF’, parallel and perpendicular to the c-axis 
of the crystal. (After Griffel and Stout 1950.) 


and 23 000 oe, contrary to the previous measurements. ‘This is surprising 
and its theoretical interpretation is hardly possible at present. 

Matarrese and Stout (1954) observed in NiF, a weak ferromagnetism 
below its Néel temperature of 73-2°k. They measured the torque on a 
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oe 


(2) 
Molal magnetic anisotropy of FeF,. x, and x, are molal susceptibilities parallel 
and perpendicular to the claxis. (After Stout and Matarrese 1953.) 


Molal magnetic susceptibilities of FeF, parallel and perpendicular to the c-axis 
of the crystal. (After Stout and Matarrese 1953.) 
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Fig. 12 


Molal magnetic anisotropy of CoF,. y, and y, are molal susceptibilities 
eae and perpendicular to the c-axis. (After Stout and Matarrese 


fo) 100 200 300 
(0) 
Molal magnetic susceptibilities of CoF, parallel and perpendicular to the c-axis 
of the crystal. (After Stout and Matarrese 1953.) 
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crystal suspended perpendicular to the (110) plane in one case and parallel 
to the tetragonal axis in another case. In the former case they observed 
that, unlike the other isomorphous fluorides, the susceptibility of NiF, at 
room temperature was greater perpendicular to the tetragonal axis than 
parallel to it and the difference between the perpendicularand parallel molal 
susceptibilities rose gradually from 1-102 x 10~* at 301-15°K to 1-890 x 10~* 
at 90-07°k. Below 73-2°K their measurements indicate that the crystal is 
spontaneously magnetized along the (100) directions, with a permanent 
moment whose magnitude is about 350 erg gauss-+ mole~’. Neutron 
diffraction measurements by Erickson, referred to in the preceding section, 
have indicated that in the ordered alignment the spins are inclined at an 


Fig. 13 


EB 
Saeae 
Scaen 


Temperature, °K. 


Heat capacity in calories per degree per mole of manganese fluoride. (After 
Stout and Adams 1942.) 


U 1 


angle of 10° to the tetragonal axis. If one assumes that the cocking of 
one sublattice is in the [110] direction and that of the other sublattice at 
right angles in the [110] direction, then there is a net moment in the [100] 
direction. The authors report that the moment of 350 erg gauss~1 mole=}, 
which is 3% of the saturation moment of the nickel ions, would on such a 
model correspond to a cocking of 2°5°. 

We now turn to CuCl,.2H,O. Detailed measurements of the 
susceptibility of this substance have been made by van den Handel, 
Gijsman and Poulis (1952). As mentioned in § 1, there is the phenomenon 
of ‘flopping’ of spins in this substance and a consequent change in 
susceptibility. These were actually observed. Friedberg (1952) 
measured the specific heat of this substance and showed that only two- 
thirds of the total excess entropy, namely (2/3) R In 2, is obtained below 
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the Néel point, the remaining third being distributed above it. The 
Weiss approximation should therefore be poor for this substance. 

A number of metals have also been investigated, but the experimental 
results are somewhat unrelated among themselves. In a-manganese, 
for instance, a small hump in the specific heat curve was observed by 
Shomate (1945) between 90° and 100°x, in agreement with the result of 
neutron diffraction, namely that the Néel point is at 100°x, though the 
hump is apparently too small for an order—disorder transition. On the 
other hand, the susceptibility measured by Shimizu (1930), Grube and 
Winkler (1936) and Bates and Pantulu (1935, amorphous manganese) 
above room temperature and by Serres (1938) between —200°c and 600°c 
showed no anomaly except in one example reported by Serres which 
showed a maximum at —100°c (173°K). These measurements show a 
slowly decreasing susceptibility above room temperature and thus indicate 
that the Néel point is in the range of low temperatures, if it exists at all, 
contrary to Néel’s suggestion referred to in § 1. Ashworth (1936) observed 
a A-anomaly in the specific heat curve at 350°c, while Armstrong (1950) 
observed no anomaly. In Cr, the Néel point predicted from neutron data 
is 200°C (473°K), but the susceptibility measurements due to Bates and 
Baqi (1936) between 100°K and 600°K and those due to Séchtig (1940) and 
McGuire and Kriessman (1952) at liquid oxygen or nitrogen temperature 
and above room temperature show that the susceptibility increases 
(rather than decreases) slowly with rising temperature, indicating that 
the Néel point is in the high temperature region, though the data are 
apparently not sufficient to give a definite conclusion. Armstrong (1950) 
has, however, observed a slight change in slope of the specific heat curve 
in the neighbourhood of 200°c. Fine, Greiner and Ellis (1951) observed 
an anomaly in Young’s modulus at 121°K, and in this and in other 
properties at 37°c. Weertman, Burk and Goldman (1952), however, found 
no anomaly in the specific heat at the former temperature. Pursey 
(1952) also observed an anomaly in Young’s modulus and rigidity at the 
latter temperature, which is close to the Néel temperature of Cr,O3, 
38°c. The former authors (Fine, etc.) report, however, that they used 
99-8% pure Cr which was annealed at 1400°c in purified helium; the 
oxygen content was believed to be 0-:05%. 

At this place, it may be convenient to refer briefly to the band-theoreti- 
cal treatment of antiferromagnetism put forward by Matsubara (1953) 
and Lidiard (1953, 1954). Following Slater’s idea (Slater 1951), Matsubara 
starts with the assumption that the Hartree-Fock potential for electrons 
with -+-spin is different from that for electrons with —spin owing to the 
presence of an antiferromagnetic superstructure and of exchange 
interaction. The energy band for each spin is then shown to split into 
two, the difference between the two kinds of Hartree-Fock potential 
being the measure of this splitting. The electrons with positive spins 
in the lower band of the positive spins prefer to be on atoms of one 
sublattice, and those in the higher band prefer to be on atoms of the other 
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sublattice, while the preference is reversed for electrons with negative 
spin in the band of the negative spins. The bands are then filled with 
electrons and density matrices calculated for the electrons with positive 
and negative spins respectively. It is proved that, under a certain 
condition for the exchange integral, which is not a sufficient condition 
to eliminate the possibility of ferromagnetism appearing, the picture for 
such antiferromagnetism is consistent. Interesting results obtained 
are that the mean magnetic moment of the atoms can be predicted to be 
small and that the specific heat anomaly at the Néel point is also small. 
In fact, Matsubara could show that the observed values of these quantities 
in Mn are related to each other by his theory in a consistent way. 
Lidiard follows the same line and remarks further that the susceptibility 
does not necessarily show any anomaly at the Néel point. 

We shall once more return to the experiments. Of solid solutions of 
antiferromagnetic compounds with diamagnetic compounds, four 
examples, MnO+Mg0O, FeO+MgO and CoO-+MgO (Bizette and Tsai 
1943, Bizette 1946) and MnF,+ZnF, (Corliss, Delabarre and Elliott 
1950), are reported. In the MnO and MnF, examples, the Curie constant 
per Mn atom is always constant but © is found to change proportionally 
to the Mn concentration, as can be expected, since 9 should be proportional 
to the Weiss molecular field. The Néel point is displaced to lower 
temperatures and becomes less definite as the concentration of MnO 
decreases. In FeO-+Mg0O, the data appear to be not so definite but the 
general tendency is the same. Solid solutions of two antiferromagnetic 
compounds have been investigated by Bizette (1946) with FeOQ+MnO 
and by Bacon, Street and Tredgold (1953) with MnO+CoO. The Néel 
temperature varies linearly with the concentration of one of the consti- 
tuents, as one would expect. But © does not so vary (table 2). The 
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latter authors have made a reasonable assumption for the three types of 
superexchange coupling constant (A,,2=A,,A,,, where A,,, for 
instance, is the Weiss field constant for atoms a arising from atoms 6, 
or vice versa) and derived the following formula for the susceptibility : 


a leGat(1=0)0,]+o(1~o)D/T 
T+[oO,+(1 —o) 0, ] wet o(1 4 o) E/T 5) 
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Susceptibility of alpha hematite. Curve 1: along a direction in the plane 
perpendicular to the trigonal axis, curve 2: along the trigonal axis. 
(After Néel and Pauthenet 1952.) 
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Spontaneous magnetization of alpha hematite. Curve |: along a direction in 
the plane perpendicular to the trigonal axis, curve 2: along the trigonal 
axis. (After Néel and Pauthenet 1952.) oy is written as M, in the text. 
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where o is the concentration of the a atoms, C, and C, are the Curie 
constants for the pure compounds and D and # are constants. They . 
found, however, that this formula does not accord well with experiment. 

There are other substances which show more or less antiferromagnetic 
characteristics as well as ferromagnetic characteristics. Pyrrhotite, 
FeS, and CrS are antiferromagnetic when the composition is stoichio- 
metric, but becomes ferromagnetic when an excess of sulphur is added ; 
a-Fe,0, would probably be antiferromagnetic if the crystal were perfect, 
but it shows a feeble ferromagnetism in actuality ; MnBi and MnAs 
change from ferromagnetic to antiferromagnetic at a certain temperature ; 
alkaline earth metal manganites and cobaltites are antiferromagnetic 
but become ferromagnetic when quadrivalent Mn-ions are introduced, 
by adding, for instance, LaMnO,, in a certain range of concentration 
of the latter. There are substances which have a Néel point in the 
neighbourhood of 1°K and show interesting behaviours in strong magnetic 
fields. We shall return to the latter substances in § 7. 

We might say a few words about the properties of «-hematite, which 
have been reviewed by Néel (1953 a) together with those of pyrrhotite and 
other substances. Néel (1949) suggested a ‘ parasitic ferromagnetism ” 
as the origin of its weak ferromagnetism, though quite what is meant by 
‘ parasitic ’ is not entirely clear at present. The magnetization of this 
substance under high field strengths has been shown by Pauthenet (1952) 
to follow the law M=M,+yxH. When the field is applied perpendicular 
to the trigonal main axis, the susceptibility y is almost constant over a 
wide range of temperature from 20°K to 950°K, while, when the field is 
applied along the trigonal axis, it drops to very small values below 20°c. 
These could be well understood if we assumed a high anisotropy energy 
along the trigonal axis, a negligible anisotropy energy in the basal plane 
above 20°c and a change of the sign of the anisotropy en-rgy b:low 20° ; 
the spins would then turn freely in the basal plane above 20°c and would 
become perpendicular to a field applied perpendicular to the trigonal axis, 
as well as to the field applied parallel ; this would give rise to a constant 
susceptibility y,. Below 20°c they would point along the trigonal axis 
and give rise to x, for the former field but x, for the latter field. According 
to Pauthenet’s experiments the saturation magnetization M, for the 
parasitic ferromagnetism appears to consist of two parts: one is independent 
of the anisotropy energy of the antiferromagnetism and the other makes 
its appearance in the basal plane above 20°c for a perpendicular field, 
but both have the same Curie temperature which is common to the main 
antiferromagnetism, namely, 1000°K. Recent measurements of ferro- 
magnetic resonance by Anderson and others (1954) confirm the second 
parasitic ferromagnetism but give no evidence of the first one. It is, 
however, very curious that the antiferromagnetism and the parasitic 
ferromagnetism possess the same Curie point and yet no common 
anisotropy energy for the first parasitic one nor any common direction 
of magnetization for both. 
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Alums of Fe, Mn and Cr have a Curie or Néel point at a fraction of one 
degree absolute, but it is not clear whether they are ferromagnetic or 
antiferromagnetic. Magnetic dipole-dipole coupling seems to play here 
a predominant role. We shall not discuss them in this article. Molecular 
oxygen crystals may perhaps be antiferromagnetic, but if so it is obscured 
by the presence of two polymorphic phase transitions in its solid state. 
Recent measurements of susceptibility by Kanda at Tohoku University 
(Kanda 1954) show that the susceptibility changes discontinuously at 
these transition points. 


§4. THEORY OF THE SUSCEPTIBILITY AT Low Frevp STRENGTHS 


The theory of the susceptibility of antiferromagnetic substances based 
on the Weiss molecular field approximation was founded by Néel (1932, 
1936, a, b), Bitter (1937) and van Vleck (1941). An account of this 
theory will be given here, together with the extensions made by subse- 
quent workers. 

The Dirac formalism shows that exchange interactions of the Heisenberg 
type are equivalent to an interatomic potential 

V,;=—iJ[1+4S, . S,], ° . . . . . (4.1) 
where S, and S, are, respectively, the spin angular momentum vectors 
of atoms 7 and j measured in units of % and J is the exchange integral. 
In ferromagnetism J is positive and denotes the coefficient of direct 
exchange coupling between neighbouring atoms whose electron clouds 
overlap with each other to some extent, if we assume Heisenberg’s theory 
of ferromagnetism to be applicable. In antiferromagnetism, however, J 
is negative and denotes, in the case of compounds, predominantly the 
coefficient of indirect exchange coupling or superexchange coupling, that 
is, the coupling via the intermediary of a negative ion such as O—— 
(Anderson 1950 a); however in this case too we have the interaction of the 
same form as (4.1) (see §8). The superexchange is the strongest between 
those atoms which lie on opposite sides of the negative ion, since then one 
of the p-orbitals (2p in the case of O——) contributes to it. 

We assume that J is the same for all the interacting neighbours, z in 
number, and replace the sum over them, 2S,, by z( S;), where (S,) means 
the statistical average of S,; over the sublattice to which j belongs. Then 
the potential energy of the atom # is given as follows, apart from an additive 
constant : 

Ve AV a eSy) Sir bass, » 21 (4.2) 
We assume further that the Landé g-factor is isotropic, that is, it does not 
depend on the direction of the spin with respect to the crystal orientation ; 
this assumption will be removed later. The magnetic moment of an 
individual atom is then t.=—gppS, and we have 

V = —2J2g-*yy*X py) = —Hey th 6 (4.3) 


where 
He eyed tate uy) (EO) E ner er wases\ | (4,4) 
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this quantity being the exchange ‘ magnetic field’ acting on atom 7 
which results from the sublattice 7. Another assumption we make in 
the following is that the whole lattice is divided into two similar sublattices, 
of which one possesses on each of its lattice points the spin of positive 
direction and the other the spin of negative direction. The existence of 
antiferromagnetic ordering has not yet been proved from first principles, 
except for the special model of a two-dimensional Ising system with a 
suitable lattice type (Onsager 1944, Husimi and Syozi 1950 and Syozi 
1950), but all the experimental evidence shows that it does exist. There 
are cases where we have four or more sublattices, but we shall at first 
confine ourselves to the case of two sublattices, more complicated struc- 
tures being reviewed at the end of this paragraph. 

The exchange field acting on each spin of the positive lattice results 
from the negative lattice and also from other spins of the positive lattice ; 
the importance of the latter has been emphasized by Néel (1948) and 
Anderson (1950 b). So we may put 

H,t=—a(u-)—( p+), We ee a ES 
and similarly, 
H-=—o(nt)— A 


where « and y are constants. If N is the number of magnetic atoms per 
unit volume, the magnetizations of the two sublattices are M*=(3.V)(u*) 
and M-~=(4N)<p.-) so that (4.5) and (4.6) can be written 


Ht=-"AM-—IM# (0025 © ee ee 

H-=—AMt—FM-, Wi. sae Ray eee 
where 

A=2a/N; D[=2y[N. 0s.) nn hay oe 


In the case where there is no external field, M* and M~ are antiparallel 
to each other, so long as A is positive and J"is not too large. (If [>A>0, 
each sublattice will divide itself into further sublattices with antiparallel 
spins.) We shall denote the common direction of M+ and M~ by 4. 
If there is an external field H, they will be tilted to each other, and a 
convenient way of defining 4 may be to take it as the direction of 
M*+/M+— M~/M-, as done by Gorter and Haantjes (1952). This direction 
coincides with that of M+t—M> so long as the angle between M+ and 
— Mis small or the lengths of these two vectors are about the same. 

It is evident from (4.1) that the direction 4 cannot be determined when 
there is no external field. The situation is the same as in ferromagnetism. 
We need the anisotropy energy to fix the direction of the magnetizations. 
In the ferromagnetic case, the direction of the magnetization is fixed along 
the preferred axis and a strong field is necessary to make it deviate from 
it. In the antiferromagnetic case, we expect the same situation, but the 
necessary field is even stronger, for the energy relative to the external 
field, —}yH?, is much smaller than that for the ferromagnetic case 
namely, —HM. Another difference in these two cases is that, if there 
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is no anisotropy energy, the magnetization vector is parallel to the applied 
field in the ferromagnetic case while the magnetization vectors of the 
two sublattices in the antiferromagnetic case are perpendicular to it. 
The latter follows from the fact that the susceptibility for the perpendicu- 
lar position is the greatest and consequently the free energy is the lowest. 
We postpone the problem of the anisotropy energy to the next section, 
and wish here to discuss two special cases. 


4.1. A parallel toH 
In this case, the magnitudes of the fields acting on the plus and minus 
spins are |H-+H,+] and |H-+H,~ |, respectively. The magnitudes of 
the thermal averages of w+ and p-, or the magnitudes of M+ and M- 
when multiplied by 41, are given by the standard theory of the molecular 
field approximation, and are 


M+=4NgppSB,(y*), with y+=|H+H,* |Sgup/kT, . (4.10) 
with a similar expression for the minus spin. Here 


B(y)= ane coth ae y— 5 coth Es Rae hee (ebb) 
is the Brillouin function, which reduces to tanh y for S=4 and to the 
Langevin function coth y—1/y for S=oo. When there is no external 
field, (4.10) gives the same magnitude of the saturation magnetization, 
M,, for the two sublattices, which can be obtained from 

M,=tN or pSBs(yo), with y=(A—L)M Sopp/kT. . . . (4,12) 

Here (4.7) and (4.8) have been used. MM, is a decreasing function of 

temperature as in the case of ferromagnetism and vanishes at the critical 
point, or the Néel point, given by 

CAR mal Pol eae hs eo jan ss (489) 


where 
C=Ng?up?S(S+1)/3k caer inal era (4,14) 
is the Curie constant. 
To derive (4.13), we assume y, to be small and replace SB,(y9) by 
4(S-+-1)y) and have from (4.12) 
M,=(C/2T)(AM,—Iy,), ke Oe Ry 
from which (4.13) follows at once. 
The susceptibility above the Néel point can be written in the form 


C : 1 
SS pe =%5 me (4,16 
X= FL? with @=140(A+T) (4.16) 
To derive this, we assume y+ and y~ to be small and have from (4.10) 
Ai a (A—-AM-—M-), etc: ©. 3... (4.17) 


since M+ and M~ are both parallel to H in the paramagnetic region, 
Solving this equation, we obtain M+=M-, and, putting M++M-=yH, 
we obtain (4.16). 

P.M. SUPPL.—JAN. 1955 D 
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The susceptibility below the Néel point can be calculated by expanding 
SBs(y*) in powers of H and retaining the first order term. Tt is essential 
here to assume that H,*+ and H,- are respectively parallel and antiparallel 
to H in order to obtain a consistent solution ; otherwise, except in the 
case where they are equal in magnitude and make the same angle with 
H, as we shall see below, (4.10) cannot be solved. We have M*—M)= 
M,—M- and 

— oon / / a 
ee a ma 5% (vo) ol | P+H4+rC see ve) (4.18) 
This susceptibility vanishes at 7'=0, increases with rising temperature, 
and reduces to (4.16) above the Néel point. 


4,2. A perpendicular toH 


Tn this case M+ and M~ have the same magnitude and incline symmetri- 
cally to H. The torque on, say, M*+ must vanish : 


[M+ x (H-+H,*) |= |M*xH—AM+x M-|=M*+H cos $¢—AM* sin 26=0, 


from which we have 2M+ sin 6=H/A, 4 he) a! a et ee 


where ¢ is the angle which M+ makes with 4. The left-hand side of 
(4.19) is the magnitude of the total magnetization, so that the suscepti- 


bility 1S xX. =constant= 1/A. a Viele te iat es ts tance (4.20) 


(4.20) and (4.16) give the same susceptibility at 7=T'y. 

Comparison with experiment. The general features of (4.16), (4.18) 
and (4.20) agree well with the experiments mentioned in § 3, if we under- 
stand by parallel and perpendicular susceptibilities those which correspond 
to H being parallel and perpendicular, respectively, to the axis of easy 
magnetization. Two quantities here are of particular importance : the 
ratio of x at 7’=0 to x at 7’=7'y in powdered specimens and the ratio of 
@ to Ty. 

In a powdered specimen of a substance with one preferred axis of 
magnetization, this axis is distributed in random directions in the specimen. 
One can decompose the external field acting on each individual crystallite 
into components parallel and perpendicular to the corresponding preferred 
axis and apply the result obtained above. The resulting averaged 
susceptibility is 

Xp=F Xut+3X oo he Gehlts fsa ne eee 
Therefore the ratio y,(7’=0)/y,(2=Ty) is } because y, vanishes at 
T=0, x, is constant and yy;=y, at Z=T7'y. Experimental values for 
this ratio are in the neighbourhood of 3 (0.69 for MnO) but usually 
greater (see table 1). This may be attributed to the following reasons : 
incorrectness of the Weiss approximation for y at 7'=T'y, existence of 
more than two sublattices, effect of anisotropy energy and anisotropy of 
the g-factor, etc. Li (1951) investigated the first point for the simple 
cubic lattice with spin 4 by taking account of the local order of the spin 
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arrangement but neglecting the interaction among atoms of the same 
sublattice (see §13). He obtained 0.769 for the ratio and algo the result 
that Ty=2-00 |J |/k, while our result (4.13) yields Ty—=1< [J |/k=3 | J |/k 
for P=0. Kanamori and Motizuki (1953, unpublished) applied the same 
method to the body-centred cubic lattice and obtained 0-723. Li gave 
for this case the Néel point only, which was Ty=3-18 |J |/k, while (4.13) 
gives Ty=4 |J |/k. We see that the value of the ratio is increased and 
the Néel temperature lowered when we go to a higher approximation. 
As regards the existence of more than two sublattices, the case of MnO, 
might be cited. The powder susceptibility of this substance measured 
by Bizette and Tsai (quoted in Bizette 1951) shows the ratio of about 
0-9. The neutron diffraction experiment (see § 2) reveals that the spins 
are in the plane perpendicular to the tetragonal axis and form two 
uncorrelated sets of antiparallel arrangement, one being formed with 
corner atoms only and the other with body-centred atoms only, and that 
the direction of the magnetization in one set is possibly at right angles 
to that of the other set. If we assume an interaction between these 
two sets, 1 and 2, and denote the coefficient of the corresponding molecular 
field as I”, we have in place of (4.17) the following equations : 


M*= (HAM, -TM,* I'M, I'M), 
ji ee = (ews i 1M + _r'M=): 


Summing up these four equations we obtain the total magnetization 
M,~++M,-+M,t+M,-. Then dividing it by H we obtain for the 
susceptibility above the Néel temperature 
C 
fer OA Teal 4 

The Néel temperature can be obtained by putting all the M’s equal in 
magnitude to each other but taking the signs of the plus and minus 
magnetizations to be opposite. Then we have T w=O(A—Tl)/4, in place 
of (4.13). The susceptibility at the Néel temperature is thus 2/(A-+I”). 
This gives at the same time the susceptibility for a field parallel to the 
tetragonal axis below the Néel temperature. To discuss the susceptibility 
for a field perpendicular to the tetragonal axis, we may suppose the field 
to be applied along the direction of the magnetization of one of the two 
sets, thus perpendicular to the other set. At absolute zero, the first 
get will contribute nothing to the susceptibility, while the second set 
will give rise to a susceptibility of magnitude 1/A, just as in the simple 
case discussed before. The powder susceptibility at absolute zero is 
therefore given by 


Le 2a 2 5) 
ee eeT! aie 844217) (24 A 
The ratio in question is thus (2+J”/A)/3, which is greater than 2 if I” 
is positive. The argument here given is due to Yosida (1953). 
D2 
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To discuss the effect of the anisotropy energy on the ratio in question, 
we might call attention to the formula (5.10) to be derived later. If the 
anisotropy constant K is large, the perpendicular susceptibility is no longer 
constant with respect to the direction of the applied field and the simple 
formula (4.21) does not rigorously hold. If nevertheless we apply this 
formula and take as x, the value given by (5.8), which applies when 4 
is along the easy axis and the field is perpendicular to it, and in which 
K/2M,2 should be independent of temperature, we obtain a ratio smaller 
than 2, since x, is unaffected by the anisotropy energy. In the majority 
of actual substances, however, in which the exchange field is much greater 
than the anisotropy field, such an effect should be negligible. The last 
of the possible reasons, that is, the effect of the anisotropic g-factor, will 
be discussed at the end of this section. 

We now come to the ratio @/T'y. (4.13) and (4.16) give 


O/T g=(A4-T)(A—Dje ee ee 


This is equal to 1 if [=0 but greater than 1 if [>0. Experimental 
values for this ratio can be obtained from table 1 and are usually greater 
than 1. According to (4.22) the ratio @/T'y could have any value however 
large, but this is not actually the case, since we are not always allowed to 
confine ourselves to the case of two-sublattices. In principle, we have to 
consider a number of possible superstructures and to compare the corre- 
sponding free energy values before applying (4.22). The model of two- 
sublattices would be valid only for a limited range of the relative values 


of A and I. Arguments of this kind have been developed by Anderson | 


(1950 b, see also van Vleck 1951). The result, for instance, for the case of 
MnO is as follows. The structure of MnO found by neutron diffraction 
experiments (fig. 1) is stable only when A>3?I, if we understand by A 
the coefficient of the coupling between those neighbours which lie at 
opposite sides of an oxygen ion and by I’ that between geometrically 
nearest neighbours. Otherwise another structure becomes more stable. 
For the former structure, which possesses four sets of antiparallel 
alignment, we can derive 


O=}(A+48I), Ty=}4; O/Ty=(A+3D)/A. 
These relations can be obtained by noticing that (4.15) and (4,17) are 
replaced in this case by 


M. AM, and Mt= pu—AM+—3rM, ete. 


eae Wy 
The maximum value of the ratio @/T'y is therefore obtained for A=3r 
and is 5, which is close to that observed in MnO. For other details of 
this kind of consideration, the reader may consult Anderson’s or van 
Vleck’s original articles and a contribution by Smart (1953). 

Also the crudeness of the molecular field theory might be responsible 
for the difference between the actual @/T7'y and its value 1 predicted for 
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I’'=0 from that theory. More elaborate statistical theory, such as the 
Bethe—Peierls—-Weiss method employed by Li (1951) predicts generally a 
curvature for the plot of 1/y versus temperature above the Néel point. 
But the theory cannot be trusted to separate this statistical effect from 
the effect of the next nearest neighbour interaction. (See also § L7.) 

Anisotropic g-factor. We have so far assumed that g is a scalar. In 
some substances, the g-factor depends on the orientation of the spins 
relative to the crystalline axes. This occurs when there is a spin-orbit 
interaction and a splitting of the electronic energy level in the crystalline 
electric field, the theory of which will be given in §9. In such a case, 
the g-factor is a tensor and moreover may depend on temperature. The 
principal axes of the g-tensor may be different from the crystallographic 
axes for individual atoms, since the symmetry axes of the crystalline 
field for a given atom does not necessarily coincide with the latter axes. 
If there is more than one kind of g-tensor in the crystal, the problem 
becomes complicated, since we have to treat more than two sublattices 
and to quantize the spins separately on each of the sublattices in a given 
external field and in the exchange field, which itself may be different for 
the different sublattices of the plus or minus spins only. This situation 
has not been clearly analysed even in the approximation of the molecular 
field, though actually the substance CuCl, .2H,O corresponds to this 
case. 

We shall simply assume that there is only one g-tensor in the lattice. 
(If not, the different g-tensors may be thought to be replaced by an 
average g-tensor.) Then the principal axes of this tensor must coincide 
with the crystalline axes in the case of orthorhombic or tetragonal or 
cubic lattices. We shall denote the corresponding principal values of 
the g-tensor by g,, Jo, J3. The reduction of this g-tensor to a scalar 
g-factor is then a simple matter. The interaction between the external 
field H and the spins is represented by the Hamiltonian p,2S,gH, and 
this can be replaced by ppg2,S;H’ if we define H’ by gH=gH’, or, in terms 
of the principal values of g, by 


9, ,=9Hy’, Golip==Gil pe aGslis=Gile,.) « *. (4:23) 


where H,, H,’, etc. mean the components of H and H’ along the principal 
axes. In other words, we can use the fictitious external field H’ in place 
of the real external field H and treat the g-factor as a scalar. 

In particular, if we are interested in the susceptibility, for instance, 
along the axis 1, that is, the susceptibility for the field applied along the 
axis 1, we have only to multiply the susceptibility x(x or x, below the 
Néel point) obtained in the fictitious (g, H’)-system by (9,/g)? to obtain 
the real susceptibility, because the susceptibility can be obtained by 
differentiating the free energy twice with respect to the component of 
the field and each time when the free energy is differentiated we have the 
factor g, for the real field but g for the fictitious field. This is because the 
free energy is a function of the interaction Hamiltonian given above. 
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We may take the axis 1 as the preferred axis of magnetization and 
axes 2 and 3 perpendicular to it. Then the powder susceptibility (the 
real one) is given by : 

Xp=4F{G1/9)?xu t+ (gal9)?xi+(9sl9)?xi}s - > + + (424) 
instead of 4+ x,. The ratio of x,(7=0) to x,(7=Ty) is therefore 
2 


(9o2-+-932)/(9x2-+922+95"), which is not necessarily equal to 3. 


§ 5. THEORY oF THE SuscHPTipitiry AT MopERATELY HicH FIELps 


We now proceed to the case where the applied field is high but not so 
high as to decouple the exchange interaction between magnetic atoms. 
We conveniently divide the subject under the following headings. 


5.1. Anisotropy Energy 


We shall introduce here the anisotropy energy in a phenomenological 
way. Let the direction cosines of «S+), or of M+ ,be «,, 8,, yi, with a 
similar notation for negative sign, and those of 4, the direction of Mt— M-, 
be «,8,y. Weshall always suppose an isotropic g-factor and consequently 
M+ to be constant multiples of (S*). If otherwise, M’s and H are to be 
primed and the methods mentioned at the end of the last section are to 
be applied. We introduce the anisotropy energy in the simplest form : 


Uniaxial case: H,=—4K(y2+y_2)=—Ky?. «SA ee 

Biaxial or orthorhombic case : 
B,=3K(B2+B2) 44K oy 2ty 2 E+ Key. (5.2) 

Cubic case: H,=—2,4K(«,4+8,4+-y,\—$K(at+B4+y4). . (5.3) 


In the uniaxial case, one can think of other expressions such as 


K'(a¢_+B,B_)+K"ysy_, 
but this reduces to (5.1) when («,-+-«_)?, ete. are neglected compared with 
unity, since Mt and M~ are almost antiparallel to each other. <A similar 
remark applies to the other cases. As in ferromagnetism, we may 
assume that the anisotropy energy is small compared with the exchange 
energy. We therefore assume that the magnitudes of the plus and minus 
magnetization vectors are not influenced by the anisotropy energy. 
Then we can introduce an ‘anisotropy magnetic field’ through the 
equation 
d£,=—H,* .S5M*—H, .dM, printer ki 

where 6 denotes an arbitrary variation of the corresponding quantity 
and it is understood that the magnitudes of the plus and minus magneti- 
zation vectors, M* and M™, are kept constant, though they can change 
their orientations. Applying (5.4) to the uniaxial case, we have 

oF —K (v45y4. —- y_oy_) 

— aia. (Hy,*80,+ Hoy t8B,+Haptdy,) Mt, 

so that 

H,,*=H,,*=0, Hy,* =(K/M*)y, =(K/My)y.., eats (5.5) 


Antiferromagnetism aye! 


where we have put M+=M-= WM, at the last equality sign, M, being the 
value of M+ for vanishing external field. In deriving (5.5) we have 
tacitly assumed that the variations of the direction cosines can take 
arbitrary values, but actually this is not the case because « pas oe 
This condition introduces, however, an arbitrary additive anisotropy 
field parallel to each of the magnetization vectors, which is not important 
in constraining the direction of the latter. 

In the biaxial and cubic cases, we obtain in a similar way the following 
anisotropy fields : 


Ha,"=0, Hy,*=—(Ky/My)B., Hat*=—(K2/My)ys, . . (5.6) 


ax 


Ag,*=(K/My)a3, gy” =(K/My)B .°, Fgg*=(K/Mo)y 4. - (5.7) 


We shall investigate two special cases. 


5.2. The Case of Uniaxial Anisotropy with Positive K 

In this case the magnetization vectors M* and M~ are normally in the 
direction of the axis of anisotropy, namely, the z-axis. When a field H is 
applied along this direction, they change their magnitudes but do not 
deviate from that direction as long as the field strength is not higher 
than the critical field strength, H,. In this case, the anisotropy energy 
evidently plays no role. 

When, however, H is applied perpendicular to that axis, we have a 
small correction to the perpendicular susceptibility due to the anisotropy 
energy. The torque acting on, say, M* is 


|M+x (H-+H,-++H,) |=M*+H cos 4—AM* sin 26—M*(K/M,)sindcos 4, 
(cos d6=y,=—y_), which must vanish, so that 

H 
A+ K/2M,?° 
The left-hand side of the equation is the magnitude of the induced 
magnetic moment, so that the susceptibility is given as 


y= 1/(A+K/2M,*). mete ae ee. (5.8) 


Similarly, we obtain the susceptibility 7, when 4 is perpendicular and 
H is parallel to the z-axis as follows : 


X= 1/(A—K/2M,”). 5. ud Me AOR) 


For the general case, we may assume H to be applied in the ae-plane 
and J also to be in this plane. It may be seen from the arguments given 
later in §5.3 that 4 cannot deviate from this plane. We define angles 
62, ¥, 0,, 0_ and ¢ as indicated in fig. 15. We decompose H into compon- 
ents parallel and perpendicular to 4, H,—H cosy and H,=H siny. The 
parallel component causes the magnitude of M* to change to M,+4y,H, 
and the perpendicular component makes it deflect from 4 by an equal 
angle ¢. We assume that the square of this angle can be neglected 
‘compared with 1, that is that M o> t(x,4,)2, and also neglect 


2M+ sin d= 
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3(x,H,)2 compared with JZ,°. In this approximation the equations of the 
balance of torque can be written as 
(M,+4xnHy)H sin (b—¢)—AM,? sin 2¢-+-K cos 0, sin 6,=0, 


(My—4xuHy)H sin (-+-4)—AM,? sin 26—K cos 6_ sin @_=0, 
where 
6,=b—04—4, 6_=—Oy+¢. 
When we again neglect x, x,,H,7/4 against M,? and put 2M, sin d=y,H,, 
we obtain 


¥1=1/[A+(K/2M,2) cos {2(¢—Oy)}] . . . « (5.10) 

and 
| (x1—xu) H? sin ¢ cos p=K sin {2(¢—Oy)}. =. «~~ (5.11) 
Fig. 15 , 


Definition of angles. 
Fig. 16 


2k A(X —Xy) 


Geometrical representation of eqn. (5.12). 


(5.10) agrees with (5.8) and (5.9) for the corresponding special cases, and 
(5.11) can be written in another form 


sin 204% 


va = $$ 
ae COS 204 —(x1— X,)H?/2K° © 


(5.12) 


(5.12) gives ys as a function of 6, and H if the term in (5.10) dependent on — 
K is neglected. A geometrical representation of this relation is shown 


in fig. 16. 
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An interesting result of the above calculation is a sudden change of the 
directions of magnetization in the case where H is applied along the easy 
axis. From fig. 16 it can be seen that ¢ is 0 or 47 according as H<H, or 
H>H, when 0,=0, where 


dei yi yy huee 2 3) el hae (Dek) 


Here the K-dependent term in (5.10) is (and should be) neglected since 
it is of the order of (y,—y,)H2/4M,2 and must be considered small 
compared with A, or the reciprocal of y,, in our approximation when H 
is of the order of H,. A simpler derivation of (5.13) was given in § 1, 
comparing the values of the free energy at the parallel and perpendicular 
positions. Thus the direction of the magnetization turns at H=H, 
from the direction of the easy axis to that perpendicular to it. This 
phenomenon of flopping of spins was first predicted by Néel (1936 a) 
and was found by Gorter and his collaborators (1951) to occur actually 
in the crystal of CuCl, . 2H,O. 

Another result which follows from the arguments above is the field- 
dependence of the susceptibility. We shall now assume H?<H,?. The 
components of the magnetization induced in the directions parallel and 
perpendicular to A are y,Hcos% and y,H sing, respectively. 
depends on H? through eqn. (5.12), so that the induced magnetization also 
depends on H?. (5.12) can be written in another form 

(X1—Xu)H? sin 20y7 
2K —(x,—x)H? cos 207 Ae 


H 
tan 26= sin 20(1— He cos 20) : 
(5.14) 


where 6=y—6, is the angle between 4 and the easy axis, and thus 
for H?<H,? we have tan 26=(H?/H,’) sin 20,. The parallel and per- 
pendicular components of the induced magnetization are therefore given 
approximately as 

xy {cos 6,—(H?/2H,?) sin On sin 207}, 

x {sin 6, —(H?/2H ,?) cos Oa sin 207}; . ° ° (5.15) 
respectively. Here we have to take (5.8) for x,, if the anisotropy energy 
is not at all negligible, because 4—O, is small. Now if we are interested 
in the component of the induced magnetization along the applied field, 
we have to multiply each of (5.15) by cos % and sin #, respectively, and 
add them 

xu cos? b-+ x, sin? ¢ 
=H [x cos? Og-+ x1 sin? O_g+2 (x1 — Xu) sin* Oy cos? 6, (H?/H,”)]. 
The average of this expression over all possible crystal orientations yields 
the powder susceptibility 


1 Dy 4 1 2 a a 
Xp= 3 xt 3 Xi TF (yi — xn) (H?7/H 7) = 3 xt 5 xt TERK (vi 
(5.16) 
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This susceptibility formula was first published by Nagamiya (1951 a) 
and Yosida (1951), though Néel appears to have known it earlier. 

Bizette and Tsai (1939) observed in MnF, the field-dependence of the 
powder susceptibility, namely, y,—180x10~-* for H=8000 oe and 
¥,=190x 10-® for H=26 000 oe at the temperature of liquid hydrogen. 
If one assumes y,—0, one obtains by applying (5.16) a value of K= 
6-0 x 10° ergs/em3 for the absolute zero of temperature. The corresponding 
value of the critical field strength is H,=6-6x10* oe. These values 
roughly agree with those which can be derived from experimental data 
for the anisotropy of the susceptibility above the Néel point (see § 11 
where H,=1 10° oe is obtained). See, however, the remark on p. 21. 

For MnO, one estimates from the measurements of Bizette, Squire 
and Tsai (1938) that K=5-7 x 104 ergs/cm® and H,=3-6 x 10* oe. 

Néel (1953 b) proposes another interpretation of the field-dependency 
of the susceptibility for crystals of cubic symmetry. He supposes 
domain walls such as found in ferromagnets, though little is known 
about the domain structure in antiferromagnets. If one side of a domain 
wall has the magnetization direction 4 parallel to the wall and the other 
side that perpendicular to the wall, both being assumed to be the easy 
directions of the crystal, one has on applying an external magnetic field 
a displacement of the wall; the domain in which 4 is more perpendicular 
to the applied field increases its volume and the other domain contracts. 
If c is the elastic. constant, per unit area of the wall, of the elastic force 
which pulls back the wall to its original position, and if one assumes the 
existence of parallel domains of a constant thickness Z with alternating 
direction of 4, one can easily derive the following H-dependent part of 
the powder susceptibility 


1 
3Lec (x1 — xn)? H?, . . . . . . . (5.17) 


which has the same form as (5.16). Néel suggests by analogy with 


ferromagnetism that the domain effect would override the effect of | 


anisotropy energy. 


5.3. T'he Case of Orthorhombic Anisotropy 
We now take the anisotropy energy to be of the form 
Ba=2K, (BP +B?)+4K9 (yy +y?), Ke>Ky>0. . (5.18) 


If H and J are confined in the ay- or xz-plane, the effect of the anisotropy 
energy is exactly the same as in the uniaxial case, provided we now take 
the x-axis as the preferred axis. For the general case the problem can be 
formulated by writing down the free energy of the system as follows : 


F=—}x Hy’—4x, H,2+-anisotropy energy 
=— 2X A? (ag +PBy+yyn)?—}x. A? [1—(a%g-+PBut+yyu)?] 
TAP + Koy?, PR rt ye 
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where tH, Po, Yq are the direction-cosines of H and the anisotropy energy 
is simplified as shown. We look for the minimum of this free energy 
expression subject to the condition «2+62+ y2—1. The problem is the 
same as that of looking for the longest principal axis of the ellipsoid 
F=1 in the «, 8B, y-space. We have only to solve the equation 


oe — A toy “HY A 
tnx Bo?+ki—A Baye == 0) 1. Ae eon (20) 
nV Ha Puyn Va tKe—A 


to obtain the stationary values of F, i.e. 
A°— (1K y+K2) A? +[(1ey+K9) Oe + KPA +Ky Yu? +k, Ko] 


KN es Opt — Oe OO. us. (D- 21) 
where 


y= 2K 4 /[(x1—Xu) 47], ke=2K9/[(xu— Xu) H?]. . . .  (5.22) 


The lowest root of (5.21) is the desired one, and the corresponding direction 
a, 8, y can be calculated from the elements of the matrix of (5.20) after 
substituting this lowest root for 2. 

By a simple calculation one can show that, when H is in one of the 
coordinate planes, there is such a root which does not depend on the 
direction of H. For instance, when H is in the xz-plane, one gets a root 
A=k,. We have therefore the possibility of a crossing of two of the roots 
when the direction of H is varied in that plane. In the case of the 
xz-plane, it can be shown easily that the crossing occurs for the two 
lower roots along the hyperbola 

i ee He 1 

ote 2 ko) yi 
and that the magnetization direction «, 8, y changes abruptly when the 
field strength crosses over this hyperbola; this direction is confined 
in the xz-plane for field strengths inside the hyperbola but coincides with 
the y-direction for those outside the hyperbola. We shall call (5.23) 
the critical hyperbola. Its vertex evidently coincides with the critical 
value of H given by (5.18) if we put K,=K. When the applied field 
is in the xy-plane, we have a crossing of the two higher roots along an 
ellipse, but this does not interest us. For the yz-plane, there is no 
crossing. It can be said in general that for the applied field outside the 
xz-plane no abrupt change of the magnetization direction can occur. 

The matter can be understood in another way by using a geometrical 
representation. We introduce oblique coordinate axes y, z, H and their 
reciprocal axes 7, £,2. (The reciprocal axes are those with which we are 
familiar in crystallography.) The ellipsoid F=constant can be written 
with respect to these reciprocal axes as 


Ky? + KC+E (xy — xn) WP? = 1. Geta oes (0.24) 
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The longest principal axis of this ellipsoid is that required, and it deter- 
mines the direction of the magnetization, 4. Figure 17 represents the 
ellipsoid when H is in the xz-plane. The semi-axial lengths, a, ¢, indicated 
in the figure and that along the y-axis, 6, are, respectively the inverse 
square roots of 3(y,—x,) H?, K, and K,. Therefore an increase in H 
flattens the ellipsoid along the a-axis and makes the principal axes 
rotate in some way. If H is in the zz-plane, the ellipsoid is symmetrical 
with respect to this plane so that two of the principal axes lie in this 
plane and the remaining one is along the y-axis. With varying H the 
former two rotate in the wz- or x«é-plane and change their lengths at the 
same time, while the latter one remains unchanged. For small H the 
longest principal axis is evidently that which is close to the x-axis, but 
as H increases a value is reached for which the principal axis along the 
y-axis becomes the longest, since K,>K, and hence c<6. This explains 
the sudden flopping of the spin direction mentioned above. If H is in 
any other direction, nothing like this can occur, the longest principal axis 
changing its direction continuously with changing field intensity. 


The ellipsoid given by eqn. (5.24). The ellipse A represents the cross section of 
this ellipsoid with the xz-plane for H>H, and the ellipse B represents 
that for H<H,. a and ¢ represent semi-axial lengths, and the dotted 
circle represents the sphere of radius 6 which touches the ellipsoid at 
the j-axis. The field is applied in the xz-plane. 


The critical hyperbola (5.23) was first derived by Gorter and Haantjes 
(1952) for absolute zero (y,=0) by a.different method, and its existence 
in actual substances was proved by susceptibility measurements by 
van den Handel, Gijsman and Poulis (1952) for CuCl, .2H,0. 
A simplified derivation described above, with a generalization to finite 
temperatures, is due to Nagamiya (1953, 1954). Gorter and Haantijes 
assumed an anisotropic Weiss field coefficient, namely (—A,M Be, 
—A,M,~, —A,M,>) for the exchange field H,*, but this is equivalent 
to our formulation as long as the applied field is small compared with 
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the exchange field so that the plus and minus magnetization, vectors are 
almost antiparallel to each other. To see this equivalence, we write the 
free energy due to the exchange field as 


A eM uM ame i AM 
=A, (MM, +N M,-+M5M,)+(4,—A,) UM, 
The first term on the right-hand side represents the free energy corre- 


sponding to the isotropic Weiss field coefficient and the remaining terms 
represent the anisotropy energy, which can be written as 


—K1B,8_—K oy sy_, 
with K, and K, proportional to M,?, if the plus and minus magnetizations 
are almost antiparallel and their magnitudes do not differ greatly from 
each other. Gorter and Haantjes treated also the case where the applied 
field is of the order of, or greater than, the exchange field. We do not 
consider this case here (but see § 7). 


§ 6. ANTIFERROMAGNETIC RESONANCE 


6.1. Theory 


We have so far confined ourselves to static problems. To go over to 
the dynamical case, we have to set up the equation of motion. This can 
be done in an approximate way by extending Kittel’s equation of motion, 
familiar in ferromagnetic resonance problems, which states that the rate 
of change in the angular momentum of the whole system is equal to the 
torque acting on the magnetization vector of that system. We apply 
this principle to each of the two sublattices of an antiferromagnetic 
crystal and have 


id 
panda ic aie hte) Y=9bp/h. Sey» * (6.1) 
We define 


M—M+-.M-, M’/=M+—M-. ee (6-2) 


Then, adding and subtracting the two equations (6.1) and utilizing the 
relation H,t+=—AM-—I™M?, etc. for the exchange fields, we obtain 


aS M=MxH+Mx4 (H,++H,-)+M’x3 (H,*—H,),. (6.8) 
v6 

+f M’'=M'x(H—AM)+M x4 (H,t—H-)+MxHHt+H,-). (64) 
Y 


In the ferromagnetic. case we have only one equation of the form of 
(6.3) except that a single term MH, replaces the anisotropy terms, 
provided we neglect the demagnetizing field and also the damping term. 

‘Moreover, save for the anisotropy field term, the equation in the ferro- 
magnetic case is a rigorous quantum-mechanical one, though the anisotropy 
term is more or less of a phenomenological nature (van Vleck 1950). 
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In contrast to this, one of our equations, (6.4), depends explicitly on the 
approximation of a molecular field, as one sees from the appearance in 
it of the constant A. One may therefore feel uncertain of its validity, 
particularly near the Néel point. To avoid that difficulty as far as 
possible, it may be advisable to make use of the relation A=1/y, and of 
experimental values for y, and x), which will appear in the following 
in a combined form 1—x,/x,. The way thus followed actually gave very 
satisfactory results when applied to a substance, CuCl, . 2H,0, the only 
substance hitherto investigated experimentally in detail, as we shall see 
later. 

We shall confine ourselves to the case of orthorhombic anisotropy 
energy ; the case of uniaxial anisotropy will follow as a special case and 
the case of cubic symmetry is equivalent to the uniaxial case when the 
direction of M’ is in the neighbourhood of one of the preferred axes. 
In the orthorhombic case, the anisotropy terms appearing in (6.3) and 
(6.4) can be written in components by use of (5.6) as follows: 


K mais / "6 K , , 
(6.3) al sae (MMH, M, ); ry ie at Ue Mu, if 
WK 
ap eM), 
VG ss , , K , ’ 
(6.4) are 7 (M,M, ae 2 M,), ze Mel » +M, M ,), 
ie 
~ ay (Mell +UsM,), . . . . . . s (68) 


where we have put (8,+8_) M)»=WM,, (8,—B_) Myp>=WM,/ and (y,+y_) My 
=M,, (y,—y_) Mo>=M,'. We notice that those terms in (6.5) which are 
quadratic in the components of M can be neglected compared with 
those quadratic in the components of M’, since we assume here, as in 
the preceding section, that the square of the induced magnetization 
M is small compared with the square of the spontaneous magnetization 
of each sublattice or of M’, 

Now we define 5M and 5M’ as arbitrary small variations in M and M’ 
and wish to preserve the latter letters for their static values, namely, 
their values at the equilibrium configuration in a given static field H. 
Then, taking the variations of eqns. (6.3) and (6.4), we obtain the equations 
for small oscillations. In doing so, we notice that H—AM(M: the 
static value) can be written as 


H—AM=H,+H,—4 (xu + xt) =(1— xu/x1) Ay 
=«M’ (M’H)/42/,2, er eee a 
a=1—y/x,, : . . . . . . * (6.8) 
where we neglected the anisotropy field compared with the exchange 
field and put simply A=1/y, (see (5.10), which is the general formula for 


X1 in the uniaxial case, including the anisotropy field). In this approxi- 
mation of neglecting H,, compared with H,, we see also that the anisotropy 
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terms in (6.4), namely (6.6), can be entirely neglected compared with the 
first term in (6.4), since they are in the ratio H,: H 2 a8 to their order 
of magnitude. If this term is not neglected, the treatment becomes 
complicated, except in one case, (6.1.1(a)) below, where the calculation is 
still quite simple. This term is therefore neglected except in this case. 

The equation of small oscillation now comes out by use of (6.5) and 
(6.7) to be 


2 ‘Se (M,'3M,'+M,'3M,’), 
1 asm ——HxSM+ +42 (M,/3M/+M,3M,), - + (6.9) 
y dt 2M? 

x oh (M1, '3M,'+M,/3M,'), 
- (aa x8’ AM XM. genes Sates <= 6.10) 


In these equations we have to substitute the static values of M and 
M’ derived in §5. We shall apply these equations to several special cases. 


6.1.1. The Applied Field is Parallel to the x-axis 

In this case we have to distinguish three cases: H<H,, H>H,, and 
H=H,. At H=H, the axis of magnetization, that is, the direction of 
M’ turns in the xy-plane and the frequencies of the eigen-oscillations of 
the system vary with the angle through which it turns. For a given 
frequency of an applied oscillating field we can therefore expect a resonance 
to occur at H=H, if the range of those frequency variations cover this 
given frequency. We shall come to this case later. 

(a) H<H,. The vectors M and M’ both point in the z-direction, so 
that (6.9) and (6.10) can be written as 


y 8M, =HdM,+(K,/M,) 5M, 
y 8, a —HbM,—(K,/M) OM 5, 
y 8 M,/=cHSM,/+2AM,$M,, +(K,/2M,2) x,H5M,'+(Ko/M,) 8M, 
y-18.M,,' =——¢HbM,’—2AM,8M,. | —(K,/2M,*) x, H5M,'—(K,/M,)sM,. 
Terms written to the right of the vertical bar are due to (6.6) which has 
been neglected. If we take account of them, we cannot put A=1/y,, 
so that « should be 1—Ax,. Neglecting them for a moment, and taking 
. 6M,, 5M,, 5M,', 5M,’ proportional to exp (wt), we obtain the deter- 
minantal equation 

—H —iwly —K,/M, 0 

0 2AM, —iw/y oH 
—2AM, 0 —tH —ta/y 
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This reduces to 
(w/y)4—(w/y)?[(1-+02)H?-+¢,+¢9]-+ 02H 4 — H?(acy +g) +602=9, « (6.11) 


where 


C;=2AK,, ¢g=2AK,. oh vo otn ee eee 
We therefore have 
(2)! spate peytegk (a8 22 (1+), 09) 
4 “t-(C4—Ca)*} 48]. se) 5 ee eee 


For the special case of c;=c, and absolute zero (where «=1), (6.13) reduces 
to 
w/y=(2AK)?+H, (K,;=K,=K) ce et ees 


We shall now take account of the neglected terms. It is easy to see 
that we have only to replace «2, «c,+-ac, and c,+c¢, in (6.11) by «1%, 
%oC4-+0%,C, and c,+c.—(K,—K,)*/M,?, respectively, where 


%=1—(A—K,/2M,") xy, %-=1—(A—K,/2M 9”) xu 
and c, and c, are now 
¢4=2AK,+(K,/M,)?, ¢-=24K,+(K,/M,)?. 
We shall write down here only the formula that corresponds to (6.14) 
w/y=[24K+(K/M,)?})?+A=[(28,+4,)0,)*+H, Be, (88 sy 


where H,=AM, and H,=K/M,. (6.15) was first derived by Kittel 
(1951). 

(b) H>H,. In this case M is parallel to the z-axis but M’ is perpen- 
dicular to it and points parallel to the y-axis. The equations are 


y8M ,=—(K,—K;)/M, .8M,', -y18M,'=— 2AM, SM,, 
y18M, =H8M,, y 18M,’ =0, 
y 8M, =—H8M,—(K,/M,)sM,', y-8M,/ =2AV, SM, 


These can be separated into two sets, one for 5M, and 6M,’ and the other 


for 6M,,5M,andé6M,’. (If we did not neglect (6.6), we should not have 
got such a separation.) The corresponding frequencies are 


(w/yP=cg—ey ww ww wt ww (OG) 
and (w/y)?=H?—c,. ee a yf! eh 


(6.16) corresponds to a mode in which the projections of M+ and M- on 
the xy-plane oscillate symmetrically with respect to the x-axis and the 
vector M’ does a rotational oscillation with the x-axis as the axis of 
rotation. (6.17) corresponds to a mode in which M+ and M~ are rigidly 
bound to each other and the vector M does an elliptic rotation about the 
«-axis, keeping M’ parallel to the y-axis. The former can be excited by an 
oscillating magnetic field parallel to the z-axis and the latter by an 
oscillating field in the yz-plane, particularly one parallel to the z-axis. 


EE ee 
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6.1.2. The Field is Applied in the xy-plane 


This case is fairly complicated because the angle between H and M’, 
that is, , varies with the intensity of the applied field. Yosida (1952) 
obtained the following equation : 


w\4 w 
(*) eo (“)' [ H?(«2 cos? p+1)+¢,+¢, {cos? (4—67)—2 sin?(#s—O77) }] 


+ H*a2 cos? — H? [e, {a cos* } cos 2 (b—07) +a cos yb sin 0, sin (s—O) 
+sin ¢ cos Oy sin (Y—6z)}+c, (« cos? 4—sin? %)]-+c, cos 2 (s—Oq) 
X [¢o—c, sin? (J—6,)|=0. Some Mere da oi. Uke (OELG) 


In this equation and 6, are related to each other by (5.12), K being 
replaced by K,. 

An approximate form of (6.18) has been obtained by Nagamiya 
(1951 b) and Keffer and Kittel (1952) neglecting (H?/c)? compared with 1 
but retaining H?/c, where c=c,=c, because they treated the case of uni- 
axial anisotropy energy. (The H‘-terms given in the former’s paper are not 
correct except in the cases of 0,=0 and 6y=7/2.) Eliminating % by 
(5.12) with K=K,, and retaining terms up to the order of H*4, we obtain 

w 


4 2 
(<) —(2) [c,+c,+ H?(1-+-«? cos? 6,)—«?(3-+2«) c,-1H4 cos? Og sin? 67] 


+¢,¢.—H? [(c, +c.) « cos? 6_—C, sin? 6g|+ H* [«?+-«(2c,c,;-1—1) 
xX sin? 6,,| cos? 0,=0. Ss A eae SO ee SR ee de 8) 


This equation is valid when H is small compared with H,. 


6.1.3. The Field is in the xz-plane 

If the field strength is inside the critical hyperbola (5.23), the same 
equation as (6.18) applies with suffixes 1 and 2 interchanged. On the 
other hand, if we are outside the critical hyperbola, we obtain after a 


simple calculation 
H,2 H2 


SF Oooo 6.20 
elo? —)—-el vite) 
6.1.4. The Field is in the yz-plane 
This gives 
ee lO 
—_— — + ——,-—_ =, 6.21 
(aye, * (oly es ee 


6.1.5. The Critical Field Resonance 

When the field is applied along the w-axis and its strength is just equal 
to H,, we have to put 0,;=0 and let > be arbitrary in (6.18). The eigen- 
frequencies vary with the latter angle. If w, is the frequency of the 
exciting oscillating field, the resonance condition w=az is satisfied at a 
certain value of y% for a certain range of the values of the parameters 
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contained in the frequency equation. We shall not give here the mathe- 
matical details, but rather indicate the result by fig. 18, which was 
obtained by Nagamiya (1954). 


Fig. 18 


(w/a)? x10“ 
4 


2 

Flak 

Diagram for obtaining the range of temperature in which resonance takes place 
at the critical field when the field is applied along the easy axis a. Thick 
parts of horizontal lines show that range. The heights of these horizontal 
lines correspond to assumed different values of (wo/y)?, where w, is the 


applied circular frequency. This diagram is applicable quantitatively to 
CuCl, . 2H,0. 


6.2. Application to CuCl, . 2H,0 
Detailed measurements of the antiferromagnetic resonance in this 
substance have been made by Ubbink (1952, 1953). We shall compare 
the theoretical results obtained above with his results. 
In this substance the g-factor possesses a pronounced anisotropy. 
Paramagnetic resonance measurements by Itoh, Fujimoto and Ibamoto 
(1951) at room temperature yielded 


Ja=2195, g,=2:075, g,=2-260 


while susceptibility measurements by van den Handel, Gijsman and 
Poulis (1952) above the Néel temperature give 


Jag=2:20, gy,=2:03, g,=2-24. 
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Applying the argument given at the end of § 4 and making use of the data 
on the susceptibility below the Néel temperature, given by the latter 
authors, one obtains a curve for « (fig. 19). Theoretically the ratio 
C,/c, should be independent of temperature, and its value, determined in 
the course of comparison between theory and experiment for the resonance 
absorption, is 3-3. It then remains to determine the parameter cj. 
This can be obtained from the value of the critical field strength along the 
“-axis (cy=aH,?) for the range in which it is measured by antiferro- 
magnetic resonance, susceptibility and proton resonance measurements. 
At lower temperatures, where resonance occurs at a field strength different 
from H, formula (6.17) affords its value. Figure 20 gives c, as a function 
of temperature. 


Fig. 19 
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Curve for «=1—y,/y, versus temperature for CuCl, .2H,O. Experimental 
points are due to van den Handel, Gijsman and Poulis (1952); x and + 
from the susceptibilities y,,, and xoq with the formula «=1—(95/9,)?Xaq/Xba 
with measuring field strength of 3600 oe and 6000 oe, respectively, and 
A from x,, and x with the formula «a=1—x,,/xa, where xoa, for 
instance, means that the field is applied along the b-axis and the spin 
direction is along the a-axis. 


We can now proceed to the comparison between other details of the 
theory and experiment. Figure 21 shows two resonance values of H as 
functions of temperature measured by Ubbink with a 3 cm wave. The 
upper curve has been used to determine the parameters as mentioned 
above and the lower curve has been calculated by (6.13) with the minus 
sign for its main part and with the plus sign for the rapidly ascending 
part near the Néel point. In the upper curve, one should see a kink at 
2-25°K, above which temperature resonance takes place along the critical 
field curve and below which it follows (6.13), but Ubbink’s experiment 
is not accurate enough to confirm the existence of thiskink. Aremarkable 
fact that theory predicts and experiment supports is that at 4:22,°K 
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Resonance field strength as functions of temperature when the field is applied 
along the a-axis ; experimental points are after Ubbink (1953 a). 


Curve of c, (=2K,A) versus temperature for CuCl, . 2H,0. 
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the condition (6.16) is satisfied and that above this temperature resonance 
along the critical curve vanishes. We therefore have at 4:22,°K a 
resonance curve semi-infinitely broad, being cut off for H<H, and extend- 
ing indefinitely for H>H,. Such a peculiar curve has actually been 
observed. 

Figure 22 shows the case where the static field is applied in the xz-plane. 
The curves lying in the lower part have been calculated from (6.18), 
interchanging the suffixes 1 and 2. In the upper part we have two types 
of curves. For 7’<2-25°K the critical hyperbola lies lower than the 
hyperbola given by (6.20), so that resonance takes place along the latter. 


— HH,» kilo-oersted 


3.95 5° AW 


aia 3760 K 
x 1.895 °K a 3.883 °K 


12 10 8 6 4 2 
Hz, kilo-oersted 


Resonance field curves in the ac-plane ; experimental points are after Ubbink 
(1953 a). 


For 7>2-25°K the situation is reversed near the x-axis, so that resonance 
takes place at first along the critical hyperbola and then follows (6.20) 
from the point where these two hyperbolas cross each other. At 4-22,°K 
the latter hyperbola degenerates into the z-axis and one then, has the 
situation mentioned above. Figure 23 shows the case where the static field 
is applied in the zy-plane. All these curves are for a 3 cm wave. We 
might also mention explicitly that in all these figures the field strengths 
are actual ones, not H’, although the resonance formulae have been applied 


with H’ in place of H. 
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6.3. Survey of Theory and Experiment on Antiferromagnetic Resonance 


The first experimental work on antiferromagnetic resonance is that due 
to Trounson, Bleil, Wangsness and Maxwell (1950). They measured the 
absorption of a microwave of 9300 Mc in Cr,O, and observed that the 
height of the paramagnetic resonance peak above the Néel point, 47 °9, 
increased when the temperature was lowered towards it. Below the 


Fig. 23 
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Resonance field curves in the ab-plane ; experimental points are after Ubbink 
(1953 a). 


Néel point, resonance suddenly vanished and only a broad residual 
absorption was observed. Maxwell and McGuire (1953) carried out 
further measurements with a number of antiferromagnetics, including 
Cr,03, MnF,, MnTe, MnS, MnO and other oxides, and found similar 
results. In the manganese compounds other than MnF, a maximum of 
the peak height was observed at a certain temperature above the Néel 
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point, but this temperature practically coincided with Ty in Cr,O0, and 
MnF,. The absorption peak broadened as the peak height decreased. 
They observed even less residual absorption below the Néel point than 
in the first experiment. The theory of the line broadening was suggested 
by Tsuya and Ichikawa (1951) and the behaviour of the height of the 
absorption peak, h, as a function of temperature, above the Néel point 
was shown by Wangsness (1953 a) to follow the rule 


A(T) [1—(oloy)I[x(Z)/x(T'y)I, 


at least in the case of MnO, where o is the short-range order parameter 
and oy its value for 7=Ty. 

The theory of antiferromagnetic resonance described in this section 
was initiated by Kittel (1951) and Nagamiya (1951 b) almost at the 
same time. Both assumed a uniaxial anisotropy energy and confined 
themselves to the case of H<H,. Kittel treated the special case of H 
parallel to the preferred axis at absolute zero, while Nagamiya gave an 
equation for the frequency for an arbitrary direction of the applied static 
field and arbitrary temperature. Keffer and Kittel (1952) obtained the 
same result as that of Nagamiya and went further on to the case of 
orthorhombic symmetry. They considered the demagnetizing effect, 
but this may be of minor importance in common antiferromagnets. 
The general theory for the case of orthorhombic symmetry, applicable 
to higher field strengths, is due to Yosida (1952 b) whose analysis appeared 
almost at the same time as Keffer and Kittel’s paper. He missed, how- 
ever, the resonance at the critical field, which was supplied by Nagamiya 
(1953, 1954). As seen inthe preceding subsection, the Yosida—Nagamiya 
theory has been successful in describing the experiments on CuCl, . 2H,0 
done by Ubbink, Poulis, Gerritsen and Gorter (1952, 1953). Ubbink 
himself explained a part of this experiment by the use of his own theory 
(Ubbink 1953 a, b, Ubbink, Poulis, Gerritsen and Gorter 1953 and Gorter 
1953) for absolute zero and also making use of Yosida’s theory. The 
- first report of this stimulating experiment appeared in 1951 (Poulis, 
van den Handel, Ubbink and Gorter 1951). The most complete account 
of it is found in Ubbink’s thesis (1953 a), which is largely duplicated in a 
paper in Physica (Ubbink, Poulis, Gerritsen and Gorter 1953). The 
Gorter-Ubbink theory is identical with ours as long as the exchange 
field is assumed to be much greater than the applied field and anisotropy 
field. They, however, treated also the case where the applied field is 
of the order of, or greater than, the exchange field, though they confined 
their arguments to absolute zero. We have so far not considered this 
case, but a few remarks will be given in the next section concerning the 
static problem at such high field strengths. 

All these theories predict a strong dependence of the eigen-frequencies 
upon the direction of the applied field relative to the orientation of the 
crystal. Therefore, in experiments with powdered samples, we must 
expect a broad averaged-out pattern of the resonance. One difficulty 
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encountered at the present time is that the resonance frequencies lie in 
general in the sub-millimeter region for most antiferromagnetics. An 
exception is CuCl, . 2H,O, which has a Néel point as low as 4:3°K and 
consequently the exchange field in it is very small. The resonance thus 
comes into the ordinary microwave region. There are a few more sub- 
stances which have a Néel point even lower, but no resonance measure- 
ments have yet been made. 

The spin-wave theory also gives resonance frequencies at absolute 
zero for a field applied parallel to the preferred axis (see Part III). 


§ 7. TRANSITION FROM ANTIFERROMAGNETIC TO PARAMAGNETIC 
Srate Unprer A Hieu FIELD 


The expression (4.13) for the Néel temperature was derived on the 
assumption of vanishing external field. For non-vanishing external 
field one will have non-vanishing sublattice magnetizations and conse- 
quently a different value of the Néel temperature. We shall give here 
the theory of this for three special cases and review the corresponding 
experiments. 

(a) Field parallel to the preferred axis x« and H<H,.—In this case, 
M+ and M~ are inline with H ; M~ decreases in magnitude with increasing 
field strength, then becomes parallel to H (and so to M*), and finally 
coincides with M+. This coincidence occurs continuously, and the 
corresponding antiferro-paramagnetic transition is of the second kind 
except when J’<<0 and the temperature is sufficiently low (Garrett 
(1951 b) for P>0, Kanamori, Motizuki and Yosida (1953) for the general 
case, Sauer and Temperley (1940) for a particular case, and Ziman (1951) for 
I'=0 with the Ising model and the Bethe approximation). We shall not 
enter into the mathematical detail of the nature of the transition but will 
derive the relation between the applied field and the transition tempera- 
ture. The latter may be called the Néel temperature under an applied 
field H and will be denoted as 7'y(H). 

If M+, M~ are the algebraic lengths of the vectors M+, M~>, taken 
positive when the latter are parallel to H and negative when antiparallel, 
we have from (4.7), (4.8) and (4.10) 


= Mo [((H—AM-—I'M*+) Sguz/kT], etc.,' Myy=NSouz/2. . (7.1) 
Therefore 
Bs-(M+/Mg.)—Bs“(M-|Myo)=(A—I)(M*—M-) Sgup/kP, 


and in the limit of /+—» M- we have 


d 
aye Bs) (M+/Mog)=(A—F) Soup/hT. =. . . (7.2) 
(7.2), together with (7.1) for M+=M-, namely 

M+=M yy Bs [(H—(A+T) M+) Soup/kT] Sys ee) 


ee 
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determines the Néel temperature as a function of H. An approximate 
solution can be found by expanding the Brillouin function in powers of its 
argument, and we obtain, to order H?, 


[By(H)—T'y(0)]/T'y(0)=[88/(S+1)PBg'’(0)(8Mo?)-? yx2H?,. 2. (7.4) 
By"(0)= —(S+1) [(S+-1)?+-8?]/9088, 
where y,=A7! and 
PN Op = OAT Zeer. ee CT) 
For S=1/2 the right-hand side of (7.4) is —(4M/q92) x,2H?. The field 
corresponding to 7'y(H)=0 is H=(A+T)M,,. Further, we know that 
the spontaneous magnetization M, of each sublattice for H=0 obeys the 
equation 
My=M yy Bs [((A—L) My Soup/kT] 
= MP y(0)/T+4M yy Bs'(0) [(A—D)My Soup[kTP+ ..., 
so that (7.4) can be written in another form 
[Py(H)—P yO) VL x(0)=— $ [(Ly(0)—T)/T Mg lpr yo - Xx? H?. (7.6) 
(6) Field parallel to x but M+— M~ points along the next preferred axis y ; 
H>H,.—In this case the sublattice magnetization vectors M+ and M- 
make equal angle with the x-axis and lie in the xzy-plane for lower field 
strengths but they will coincide with the z-axis above a certain high field. 
This limiting field will be a function of temperature and, conversely, the 
temperature of this transition will be a function of field strength. Thus 


we obtain 7',(H). 
We assume the effective field acting on M+ and M~ to be given by 
{0 
H.gt=H—AM-—IM+t-+ < kM ,-+k,'M,* fy eee re 


| k,M,---k,'M,*., ete. 
When the applied field is not so strong as to change the magnitudes of 
M+ and M~ appreciably from their static value M), the coefficients 
ky, ky’, kg, ky’ are connected with the anisotropy constants K, and Kg, 
defined in § 5, by the following relations : 
(hy I Ko (oak Mg. es (7.8) 

The effective field must be parallel to the corresponding sublattice magneti- 
zation, so that we put AM*+=H,,", etc. Now, if we define M and M’ as 
the sum and difference of M+ and M-, we obviously have only M, and 
MM,’ as non-vanishing components in the present case, so that the propor- 
tionality between M* and H,,*, etc. can be written 

AM ,=2H—(A+Tl)M,, 

AM =(A—I—kyt+ky)M,’. 
If M,'40; we obtain A=A—I"—k,+ ky’ and 

M =H) A— (kb, \/2) 2... 7s (7.9) 
(compare this result with (5.9)), while the magnitude of M* is given by 

M+=Mo. Bg [AM* (Squp/kT)], - . . ~ (7.10) 
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so that M+ is independent of H. If, on the other hand, M,’=0, that is 
M+—M-=WM/2, M+ follows the equation (7.3). The transition between 
these two states gives the Néel temperature as a function of H. That is, 
T'x(H) can be obtained by eliminating M* from (7.9) and (7.10) using the 
relation M+=WM,/2. Thus the relation between 7'y and H is somewhat 
similar to that between 7’ and M). 

At absolute zero, the value of the critical field is 2A —(k,—ky’)/2] Moo, 
which is greater than that of case (a) as long as we assume I"< A (if ">A the 
two-sublattice model fails). When the field is not too strong, we obtain 
the following approximate relation : 


[7 y(H)—T y(0)]/T y(0)=[8S/(S+1)PBs’"(0)(24M 97)? x1°H?, . (7.11) 
T.,{0)=C(A—I'—ky-L ke ne 


where a simplification is made by putting A>A—JI and M,—y,H, with 
y,=1/A on the right-hand side of (7.11). The coefficient of H? in (7.11) 
is one third that in (7.4) and the Néel temperature 7',(0), (7.12), is a 
little lower than (7.5). 

(c) Field is parallel to the y-axis ——We obtain the same relation as 
(7.11) but the Néel temperature for zero field is the same as (7.5). 


7.1. Comparison with Experiment 


Poulis and Hardeman (1952) measured, in CuCl, . 2H,0, the change of 
Néel temperature with applied field, by the method of proton resonance. 
The field was applied along the a- or b-axis. In fig. 24, parabola A was 
calculated from (7.6) using experimental values for M/,? derived from 
specific heat measurements (fig. 6) ; A’ was calculated from (7.5), putting 
S=1/2; B was obtained by taking one third the curvature of A and 
adjusting its height so as to cross A at the experimental value of H, at 
Ty, and B’ was obtained in a similar way relative to A’. Parabolas C and 
C’ correspond to the case of a field parallel to the b-axis and were drawn 
parallel to B and B’, respectively. In all these curves a correction for the 
anisotropy of the g-factor has been made. From the difference of the 
heights either of A and B or of A’ and B’ at zero field we can estimate the 
value of the anisotropy constant AK, at absolute zero, which is found to be 
either 0-0044 cm~? per ion or 0-0067 cm™! per ion, respectively ; the value 
estimated from antiferromagnetic resonance (see §6 and §10) is 
0-0051 cm~! per ion, which is thus in good agreement. 

Garrett (1951 a) made adiabatic measurements of a number of charac- 
teristic quantities of Co(NH,).(SO,4)..6H,O whose Néel point is 0-084°x 
and in which the exchange field is of the order of 400 oe. He derived two 
sets of curves representing M=constant and entropy—constant in the 
7'-H-plane and obtained from them a curve that separates the antiferro- 
magnetic region from the paramagnetic region, that is, the relation between 
Ty and H. The field was applied along the preferred axis, so that we have 
to distinguish two cases a and 6. In the first case we have, theoretically, 
increasing M/H for constant H as a function of 7 in the antiferromagnetic 
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region, with a maximum at the Néel point, whereas in the second case we 
have a constant value of M/H up to the Néel temperature and a decreasing 
value above. Garrett’s experiments appear to favour the first case, so 
that it seems that the spins cannot deviate from the preferred axis in the 
antiferromagnetic state even at a field comparable to the exchange 
field; the crystalline field acting on each atomic magnetic moment 
appears to be too strong. In such a case, one would naturally have to 
incorporate the crystalline field into the effective field from the very 
beginning, but one would then obtain a function similar to the Brillouin 
function in its general behaviour, if the molecular field type of approxi- 
mation is made. 


Fig. 24 
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Decrease of the Néel temperature with applied field strength for CuCl, . 2H,0. 
+ represent the observed values when the field is applied along the 
b-axis and xX represent those when the field is applied along the G-aXxis 
(Poulis and Hardeman 1952). Curves A, B, Cand A’, B’, C’ are explained 
in the text. 


Henry (1953, 1954) measured the isothermal magnetizations of spherical 
samples of MnCl, . 4H,0 and MnBr, . 4H,0, at liquid helium temperatures 
as functions of applied field strength, up to 58000 0e. He confirmed that 
the relation (7.3) holds very well with S=5/2, where 2M* is the magnitude 
of the magnetization, and obtained at saturation 


(A+T)M.=14 000 oe for MnCl, . 4H,0, 18 000 oe for MnBr, . 4H,0. 


He also measured the magnetization, at constant field strength, as a 
function of temperature, and observed a maximum at a temperature, 
T sax, Which varied with the field strength H, though the lower limit of 
temperature in his experiment was 1-3°K and the corresponding highest 
value of H for which the maximum was observed was in the neighbourhood 
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of 20000 0e. If this 7',,,, could be identified with 7'y(H), his result 
should be explained by a theory of the type described above. Unfor- 
tunately, the anisotropy energy is not known in these substances and, 
moreover, the simple theory above does not apply to polycrystalline 
specimens. Nevertheless, if we assume the anisotropy energy to be 
negligible in his experiments and apply the case b, putting k,=k,'’=0, the 
general features of his results appear to conform to the theory. The 
7 ,-H curve begins with a parabolic shape at lower field strengths and 


appears to end at a field strength of the order of (A+J")Mo 9, given above. . 


The theoretical value for the terminal field strength is 2A Mo) and should 
thus be greater than that observed (but, at most, twice as large, so long as 
we assume A>I'). Quantitative comparison cannot be made since the 
value of x, isnot known. By extrapolating 7',(H) to H=0, one obtains 


T y(0)=1-67°K for the chloride, 2-2°K for the bromide. 


The former is in good agreement with the value obtained from specific 
heat measurements by Friedberg and Wasscher (1953), which yielded 
1-622°K. 
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Part II. Inprtrect Excuance INTERACTION AND 
ANISOTROPY ENERGY 


§ 8. SUPEREXCHANGE AND DouBLE ExcHancE INTERACTIONS 
8.1. Theory of Kramers, Anderson and van Vleck 

The majority of the substances which show antiferromagnetic properties 
are ionic crystals such as MnO, MnS, MnF,, etc. In these substances the 
magnetic ions are separated from each other by the anions O, 8, F, etc. 
Neutron diffraction studies on MnO, FeO, CoO and NiO, by Shull and 
others (1951), revealed that there exists a strong antiferromagnetic 
exchange coupling between next-to-nearest magnetic ions, which are 
separated by the intervening large oxygen ions. The Néel temperature, 
which is a measure of the strength of this exchange coupling, is 122°x for 
MnO, for instance, which is not much lower than the Curie temperatures of 
ordinary ferromagnetic substances. This fact indicates that the anti- 
ferromagnetic coupling in such substances arises from a mechanism 
definitely different from that of the usual direct exchange interaction. 

The existence of such an indirect exchange interaction is also found 
in magnetically diluted salts, such as Tutton salts and alums. The 
analysis of the experimental data for the paramagnetic susceptibility 
and specific heat at low temperatures of these substances confirms the 
existence of an exchange coupling between magnetic ions separated from 
each other by a distance as large as 5~10 A. For instance, the para- 
magnetic Curie temperature, deduced from the temperature dependence 
of the susceptibility, is 0-7°K for CuSO,.5H,O, in which the nearest 
distance between magnetic ions amounts to 5-55 A. 

The mechanism of superexchange, different from the direct exchange 
familiar since Heisenberg’s original theory of ferromagnetism, was 
proposed by Kramers as early ag 1934 to explain the magnetic behaviour 
of such paramagnetic salts. His treatment was based on a perturbation 
calculation. Then, Bizette (1946) and Néel (1948) emphasized the 
importance of this superexchange in antiferromagnetic substances, and 
Anderson (1950 a) presented a detailed theoretical treatment of Kramers’ 
idea and clarified the mechanism of superexchange. The same problem 
was treated by van Vleck (1951) using a different method. 

Anderson isolates one oxygen ion and the two adjacent magnetic ions 
as shown in fig. 25, and considers, for simplicity, two electrons occupying 
the same p-orbital of the O-~ ion and one electron each in a d-orbital of 
the magnetic ions. In this state the two electrons in O— will form a 
singlet state, so that there will be no exchange coupling between magnetic 
ions M, and M,. One can think, however, of the possibility of one of 
the two electrons of the O-- ion being transferred to one of the two 
magnetic ions, say M,, and occupying its s- or d-orbital. In this excited 
state, the unpaired electron left in O-- can couple with the electron of 
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the other magneticion M,. There should also be a strong coupling between 
the electron transferred to M, and the original electron in M@,. Moreover, 
the two electrons originally belonging to the oxygen ion have opposite 
spins, so that there will appear an indirect spin coupling between M, 
and M, through this excited state. In this way, Anderson could explain 
the superexchange coupling between magnetic ions separated by a 
non-magnetic ion. 


Ground and excited states of (MnOMn)t+. 
Mn++ O-—- Mn++ Mn+ O- £Mnt+ 
ds pp ads d,d,’ p dy 

ground state excited state 


The mathematical details of this superexchange mechanism are given 
in the original articles of Anderson and van Vleck and we might not, 
perhaps, need to reproduce them here. However, we wish to go a little 
further, and for our purpose it is convenient to give here a short account 
of the mathematical method. 

Let us denote the ground configuration by A and the excited con- 
figuration by B. Then configurations A and B are represented by the 
following wave functions : 


A: ar(1) Yp(2) tp(3) Haxl4), «ee + (8:1) 
B: baiQ)) bar‘ (2) bp(3) baol(4). 6. 4p gp ihe erie we aoe 


The perturbing Hamiltonian of this four-electron system can be considered, 
for convenience, to consist of the following two parts : ’ 


V=Vjt4V, 2... ee 


where V, represents that part of the perturbation which connects the 
ground configuration A to the excited configuration B and is, for simplicity, 
assumed to be spin-independent, and V, represents the spin-dependent 
part, which is assumed to be diagonal with respect to the orbital states. 

As mentioned before, the p-orbital of the O-~ ion is occupied by paired 
electrons in the ground state and so there is no spin coupling between the 
O~ ion and its adjacent magnetic ions. Accordingly the first order 
perturbation energy vanishes. The second order energy does not depend 
upon, the spin configurations on account of the assumption just made 
concerning V,and V,. According to the ordinary perturbation theory 
the third order energy is given by 


BaD (g,t|V | i, ry. | V Ld. s)(9,8| V |g, u) 
ir 9,8 (L;, r— Hl, (Hy, oH, u) 
— yy GUVlen@rlV 190g, t1V 19, t) 
ir (H;, pete t)” ; 


where g andi, j denote the ground state and the excited states, respectively, 
and r, s, ¢ and uw denote spin configurations. Applying this formula to 


(8.3) 
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our case of two configurations, we obtain the following non-vanishing 

term : 

(At V.|BY(B* 1V, [BY \(B“ [V,| 4%) 
(Het a") h 3” — hy") ; 

where again ¢ and w denote spin configurations. 

In his original paper Anderson used Serber’s method of the Dirac vector 
model, extended so as to include different configurations. Here we shall, 
however, follow a more straightforward procedure employed by van 
Vleck. First, with van Vleck, we shall regard the two orbitals of O-- 
as different and denote them by p and p’. We denote by ¥ [(d,d,)3(pp’)+] 
the wave function of the triplet state for the four-electron system in which 
the spins of the two electrons occupying the d,-orbital of M, and the 
d,-orbital of Jf, are parallel and those of the electrons occupying the p 
and p’-orbitals of O-— are antiparallel. Similar notation will be used 
for other states. If we recombine (d,d.) (pp’) in such a for as (pd,) 
(p’d,), we find 


£4 [(pp')(dyd2)" |= VE 4" [(pdy)*(p'de)®]+-3 % 41 [(pdy)*(p'ds)4], 
4° [(pp’)*(dyde)* = 31.2 [(pd1)*(p'de)*] +b. [(pdy)3(p’de)"]} 
+4 ea ad) (ig) olsen. lee ee) = 1 « | (8.0) 


These relations can be proved by comparing spin functions for the two 
systems. For instance, the spin function of 1 [(pp’)1(d,d,)"] is repre- 
sented by 


By= 2 (8.4) 
Ueu’ 


4 {a(1)B(2)—a(2)8 (1) a(3)B(4) —a(4)B(3)}, 


where it is assumed that electrons 1 and 2 occupy the orbitals p and p’ 
and electrons 3 and 4 the orbitals d, and d,. This can be written in 
another form 


v/2[(1/+/2){«(1)B(4)+-(4)8(1) { a-(2)8(3)-+2(3)8(2)} 

ea (1)B(4)a(2)a(3 ae 4)B(2)B(3)}] 

+ 3L(1//4){ (1) (4) —ax(4)B(1) }{B(2)ae(3) —o-(2) (3) }]. 
In this expression the portion enclosed in fee brackets in the first 
term represents the spin function of 541 [(pd,)3(p’d,)*] and that in the 


second term the spin function of #41 [(pd)'(pd.)*]. 
Spin configurations for the excited state can be represented by 


ba,” (dds P(pde)"1, — hn,* U(A141')(pa)*], 
bp? (dada )(pde)*], deg (dady')(pds)*, =» + (8-6) 


where d,’ is the orbital of M, to which an electron has been transferred 


from O—. 
Tn the excited state, V, can be expressed as 


V = —t pq(1+4 (S1- $2)). A ry ah 
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Then, by use of (8.5), (8.6) and (8.7), we can calculate the matrix elements 
included in (8.4). The results are as follows : 


(A1|V,|BX)=—V2.6,  (At| Vi | Bs)=46, | 
(43|V,|B3)=—v4.6, (A®]V,[ B28 4, | ae (8.8) 
(A®| V, | Ba*)=—2 4, 

B17 SSL eae |! cf seer 
(Be|V.|BY)=Bi|Ve|BX)=t+4 oe J 


where b represents the transfer integral 
b= | ba(1)Vebp(1) dry Mvmt hk of .41) 


and J,,4 the exchange integral between the p-orbital of the O~~ and the 
d-orbital of the magnetic ion, 


Jee ty(Lbal2) VePp(2ba(l)drydry. 2: . . . (8.11) 
We insert (8.8) and (8.9) into (8.4) and observe that 
E; =H, 3=4H(t) and H;z=Hz3=ALK(s), . (8.12) 


where AH(t) and AE(s) represent respectively the energy of the triplet 
and singlet states of the metallic ion which has received an extra electron 
from O-~ in the excited configuration, thus we finally obtain 


B (triplet) =4{1/4H(t)?—1/4B(s)}402 J, g@—H{3/ 4B (2+ 1/4 Es) WT p a 
E (singlet) =4{1/4H(t)?—1/4E(s)?}{ —3b? J, a} 
—4{3/AE()2+1/4B(s)?}0? Jy a. 
These two equations can be expressed as a single equation, as follows : 
E=}(1/4E()?—1/4E(s))b7J pa (Sa, + Sa.) —3(3/4 E(t)? + 1/4 H(s)*)b? J, g. 
‘ (8.13) 


In deriving (8.13), it was assumed that the two p-orbitals of O——- were 
different from each other and also that only one of the metallic ions could 
receive an extra electron. Since, actually, either metallic ion can receive 
the electron, with equal probability, the right-hand side of (8.13) must be 
doubled. Anderson has shown that for a pair of identical p-orbitals the 
transfer integral b must be multiplied by 1/2, except when d, and d,’ are 
also identical, in which case this is not necessary. The other assumption 
used in deriving (8.13) is that all the orbitals are orthogonal to each other 
except for the expression of the exchange integral J,,,, which can be 
negative for non-orthogonal orbitals. Anderson has also generalized his 
calculation to the case of spin S greater than $ and shown that an expression 
of the form of (8.13) does still hold. 

If, in the excited configuration, the parallel spin state is excluded by the 
Pauli principle, we must put 4H(t)= 00 in (8.13). Then we see that M 1 
and M, are coupled by an exchange integral which has the same sign as 
Jyq- If, on the other hand, the antiparallel spin state is excluded by a 
strong Hund’s rule effect, we must put 4H(s)= oo. Then we see that the 
effective exchange integral has the sign opposite to that of J an 
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From the above calculation it may be concluded that there exists a 
coupling of the form J(s, . s.) between the spins of the two magnetic ions 
and that the order of magnitude of the effective exchange integral J is 
proportional to the product of the direct exchange integral J pa With the 
Square of the transfer integral b?. We also see that the superexchange has 
a directional property corresponding to that of the p-orbital. The 
p-orbital of the O-— ion which extends towards both magnetic ions will 
have a large overlap with them. We expect, therefore, a strong super- 
exchange coupling between those magnetic ions which are situated 
oppositely with respect to the oxygen ion. This directionality explains 
the fact that a strong exchange coupling exists between next-to-nearest 
Mn ions in MnO. 

There is a considerable overlap between the magnetic ions and the 
oxygen ion, so that it is reasonable to suppose the direct exchange integral 
Jq to be negative. If the magnetic ion has a more-than-half-filled shell, 
the transferred electron will set its spin antiparallel to that of the magnetic 
ion on account of the Pauli principle, while, if the magnetic ion has a, less- 
than-half-filled shell, the electron’s spin will become parallel on account of a 
strong Hund rule effect. We expect, therefore, an antiferromagnetic 
coupling between magnetic ions in the former case and a ferromagnetic 
coupling in the latter case. In fact, it has been experimentally recognized 
that oxides, sulphides and fluorides of Ni?*+, Co?+, Fe?+ and Mn?+ show 
antiferromagnetism. It is aremarkable fact that MnTe exhibits antiferro- 
magnetism and CrTe ferromagnetism. It has been pointed out by Shimizu 
(1952), however, that some compounds of magnetic ions with a less-than- 
half-filled shell are antiferromagnetic, contrary to the prediction (see 
table 1). Shimizu has emphasized that, if the Hund rule does not operate 
so strongly as to exclude antiparallel configurations, we can not omit the 
AE(s)?-dependent term in (8.13) and also that, when the many electrons 
belonging to one Mn ion are taken into account, there are other contri- 
butions arising from configurations not considered above. Even after 
taking these into consideration, one will see that the dominant term comes 
from the Hund rule effect. There are two possible ways for the anti- 
ferromagnetic terms to overcome the ferromagnetic ones. The first is 
that the large number of antiparallel configurations may compensate for 
their energies being higher than those of parallel configurations. This 
possibility was suggested by Anderson. The other possibility, which is 
due to Polder (1951), is that the crystalline field may lower the energy of 
the antiparallel configurations. 


8.2. Other Superexchange Effects 
Slater (1953) has proposed, as another mechanism of superexchange, 
that the polarization or deformation of the oxygen ion stabilizes the energy 
of the antiferromagnetic state of Mn-O-Mn. According to his idea, the 
potential energy of an electron with plus spin is, on account of the exchange 
effect or Hund’s rule effect, effectively lower when it stays in the charge 
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distribution of other electrons with plus spins than when it stays in the 
charge distribution of electrons with minus spins, if the distributions are 
the same. Therefore, when we have an antiferromagnetic configuration 
of Mn—O-Mn, in which the spins of one Mn ion point in one direction and 
those of the other Mn ion in the opposite direction, the electrons of the 
oxygen ion with plus and minus spins will be pulled each towards that 
Mn ion which has the same spin direction. In this way the oxygen ion is 
polarized and the antiferromagnetic state is stabilized. This effect does 
not occur in the ferromagnetic state. Similar effects of polarization can 
occur if there is a covalent bonding, where, however, electrons of plus 
spin are attracted by the Mn of minus spin and those of minus spin by 
the Mn of plus spin. 

Pratt (1953), following Slater's idea, has calculated in detail the energies 
of ferromagnetic and antiferromagnetic states of (Mn—O—Mn)?", using 
orthogonalized orbitals. He took into account the same excited con- 
figuration as that considered by Anderson and van Vleck. However, 
he obtained the unreasonable result that the ferromagnetic state should 
be more stable than the antiferromagnetic one. When he used non- 
orthogonal orbitals, somewhat simplifying the model, this result was 
improved. Yamashita (1954) has also discussed the effect of the non- 
orthogonality of the orbitals of Mn and O ions, in the Heitler—London-— 
Lowdin approximation. 

Slater and Pratt’s mechanism of superexchange is a little different 
from Anderson’s in that they included in the excited configurations a 
simultaneous transfer of the oxygen electrons to the right and left Mn 
ions. 

Anderson and Hasegawa (1954) have found another superexchange 
which arises from the same mechanism as that proposed by Slater. This 
acts antiferromagnetically in both of the cases of more-than and less- 


than-half-filled shells. The effect comes from the term 
(At| V,| B,*)\(By* PV Be" Bee ee | A’) 
H= J — te et tC : 8.14 
ie (Hy,"—E)\(Hp,” B49 (eta 


which has been ignored in (8.4). Here (B,“|V,|B,“) means the matrix 
element of V, between the excited state in which the right Mn ion has an 
excess electron and that in which the left Mn ion has an excess electron. 
The excited wave functions of B, are given by replacing d, and d,’ in 
(8.6) by d, and d,", and are as follows : 
Pog” [(date')*(pdy)*}, Hay, [(dade')(pdy)4), 
Pao, Udele')(pdy), PBog (dede')(pdy)*]. 
Here the suffixes c and d are interchanged for convenience. Then, the 
non-vanishing matrix elements of (B,"| V, |B.) are 
Bist | ae | Boq'*)=+*, (By,1 | ie | Byy1)=+d*, ] 
(B,. | Ve | B,,°)=—J*, (By @ | Ve B,)=—J*, | ae 


(8.15) 


————— ee ee 


Antiferromagnetism 65 
where J* means the integral 


J* =f a, )xb(1) V aby (’ 2) g (2 )dr, Ar». 5 6 A (8. 17) 


Noticing that the transfer integrals for B, are equal to those for B,, (8.8), 
except for the sign, namely, 


(At| Violen: j= V2.6, (A?| V,| B..)=— 40, | 
(A? | V,| B,.°)=/3 . 6, 6 1Ga8) 
(a4¢ | V, | By°)=40, (42 | V, | B,)=—4), 


we obtain the following results : 

E(triplet)= —{1/4E(t)2+1/AH(t) 4E(s)} Woes*, | 

(singlet) = —{3/4H(t)?+1/4E(s)?} 4627*. 
These equations are equivalent to 

H= 7b? J* {1/A4E(t)—1/4E(s)}(s 4, 
—(1/16)b2J* {3/4 (t)+1/4H(s es ae wee (oeaO) 

(8.20) shows that the sign of the effective exchange integral is always 
equal to that of J* irrespective of whether 3d shell is more-than or 
less-than-half-full. The integral J* may be approximated by 28,4), 
where S,,,, is the overlap integral between d, and p-orbitals. On the 
other hand, the ordinary exchange integral J, is of the order of 2S, 4,0. 
Since Sq, has a sign opposite to that of S,,z,, we see that J* and J, 4 
are equal in magnitude but opposite in sign and that the term furnished 
by (8.20) operates always antiferromagnetically. 

The inclusion of this term would make Anderson’s original treatment 
come nearer to Slater’s effect of polarization. We might add, however, 
that there is also an effect of polarization arising from the excited 
configuration in which an electron of the oxygen ion is excited to its 
3s state. Pratt (1953) also calculated this effect, basing his calculation 
on a simple model, and showed that it is less important than that arising 
from the excited states considered above. 


(8.19) 


UH 
. 


pd 


8.3. Double Exchange 


Zener (1951) proposed another kind of indirect spin-spin interaction 
to account for the ferromagnetic properties found in mixed crystals of 
LaMnO, and BaMnO, (or other alkaline earth manganites), neither of 
which is known to show ferromagnetism in itself. These mixed crystals 
include Mn3+ and Mn‘+ ions, and Zener considered that an indirect spin— 
spin interaction works between Mn°+ and Mn** by the following 
mechanism. Consider for simplicity one oxygen ion and its adjacent 
two Mn ions. In order to simplify the model we focus our attention, 
as before, on the two paired electrons of O--, and assume that one of 
these two Mn ions has one electron and that the other has two electrons 
with spins parallel owing to a strong Hund rule effect. We shall call 
the state in which the two Mn ions have parallel spins configuration I. 


F2 
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This can change into configuration II, shown in fig. 26, by the transfer 
of an electron of that spin direction from O— to the right Mn (which had 
been occupied by one electron). Further, the configuration IT can change 
into another configuration III by the transfer of one of the electrons of 
the left Mn to the vacant orbital of O-~ produced by the first transfer. 


Fig. 26 
Ground and excited configurations of (MnOMn)?°*. 
Moet .02=.0 eine Mn** O-"” Mne Mn*+ O2- Mn+ 
did,’ pp’ dy d,d,’ p dod,’ dy pp’ dd,’ 
configuration I configuration IT configuration ITT 


Configuration I has the same energy as that of configuration III, and 
therefore a strong resonance will take place between them through the 
excited configuration II. Thus the ferromagnetic state will be stabilized. 
For two Mn ions with antiparallel spins, such a resonance is not expected. 
Therefore, a ferromagnetic spin coupling appears between Mn** and 
Mn‘+. Zener called this type of coupling the ‘ double exchange ’. 

Anderson (1954) calculated the energy splitting due to this double 
exchange. Each of the configurations I and III has one quartet and two 
doublets whose wave functions can be written, in the notation employed 
in §§ 8.1 and 8.2, as 


by fo? [(dydy/P(d)], $1," (dydy')*(do)], 
Prrrg ” (ai)(dede’)*], ary,” [(4s)(do4')*)- | 
As in the case of superexchange, we have the following relations between 
two systems (d,d,')(d,) and (d,)(d,d,’) : i 
U(ddy')(ds)]=v*U(d,) (ded), 
PU(dydy’)°(da)] = V/$07E(d)(dady’)]— BPE, ) dad’, | (8.22) 
PU(dydy’)*(da)] =P" [(,) (dad! 1+ VIVA (dedy'9)- | 
We assume again that the effective Hamiltonian, “; 111, giving rise to 


the transfer is spin-independent ; then the non-vanishing matrix elements 
of #; 1;1 can easily be calculated by use of (8.22), i.e. 


(I, |x 111 | IIT ,*)=6, (I," | Hy 1311 | TIT ,?)=—}0, 

(Te | Hy y11 | IIT ,?)= +/ 20, A |“ y 111 | TIT ,?)= +/ 30, (8.23) 

(Lp” | 7,111 | IIL ,?) = 30. I 
Thus, the degeneracy of the configurations I and III is entirely removed 
by the transfer Hamiltonian #; 41;;. Denoting by AH(t) and AH(s) 
the energies of the states in which two electrons on the same Mn have 


spins parallel and antiparallel, respectively, we obtain energy eigenvalues 
for the quartet state 


H,=4E(t)+b (quartet) . . . .°. . (8.24) 
and for the doublets, assuming the strong Hund rule, 


Ey,=AE(t)+}) and E,4,=AE(s) (doublets). . (8.246) 


(8.21) 


a 


EE  ——— 
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From these results, we may conclude that double exchange shows 
rather strange properties, very different from those of superexchange ; 
thus, double exchange interaction cannot be expressed in the convenient 
form of a multiple of (s,.s,) and the average energy value of the state 
of higher multiplicity (quartet) is equal to that of one of the states of 
lower multiplicity (doublets). Anderson suggests that one cannot 
therefore expect the Curie-Weiss law to hold for the susceptibility of 
substances which involve double exchange and also that, in them, the 
approximation of the Weiss molecular field cannot be justified. 

The double exchange is intimately connected with electron transfers 
and accordingly with electric conductivity. Zener has shown that the 
resonance energy 6 is proportional to the diffusion coefficient of the 
Mn** ions and thus explained qualitatively the relation between the 
magnetic behaviour and the conductivity observed in manganites. 


§ 9. OrIGIN oF THE ANISOTROPY ENERGY 


9.1. Three Sources of the Anisotropy Energy ; Anisotropy Energy due to 
Crystalline Electric Field 

In § 6 we introduced anisotropy energy in a phenomenological way. 
We shall now investigate the atomic origin of the anisotropy energy. 
The ordinary exchange energy and superexchange energy, which have the 
form of the scalar product of two spins S, and S,, do not give rise to 
anisotropy energy. ‘Two kinds of anisotropic interactions can be con- 
sidered : one is the dipolar interaction between magnetic moments and 
the other is a somewhat complicated one, called usually ‘ anisotropic 
exchange interaction’ or ‘ pseudo-dipolar interaction’. The latter is 
a combined effect of LS coupling and exchange interaction. This 
interaction was proposed by van Vleck (1937 b, 1951) and is particularly 
important in ferromagnetic metals. Beside these two anisotropic 
interactions, the crystalline field acting on each magnetic ion produces 
an anisotropy energy, expressible as a function of the spin of the 
magnetic ion: The anisotropy of the g-factor usually found in para- 
magnetic salts originates also in the anisotropy of the crystalline field. 
Thus, the anisotropy energy in antiferromagnetic substances can be 
considered to stem from three sources: magnetic dipolar interaction, 
anisotropic exchange interaction and the anisotropy of the crystalline 
field. However, the effects of these three parts are quite different in 
different substances. They also give rise to the anisotropy of the 
susceptibility above the Néel temperature, just as in paramagnetic 
substances. . 

When a free magnetic ion is brought into the crystalline electric field 
arising from the surrounding plus and minus ions, the directional 
degeneracy of the orbital angular momentum is removed, completely 
or partially, because of the crystalline Stark effect. For simplicity we 
shall here assume that the ground state is not degenerate, namely, that 
in this state the orbital angular momentum is completely quenched. 
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Then we can suppose that the magnetic moment possessed by a magnetic 
ion in the ground state is the spin magnetic moment only, namely, 
2u,8 where py, is the Bohr magneton and $ denotes the spin angular 
momentum vector measured in units of #. However, since the off- 
diagonal elements of the angular momentum L (measured in units of h) 
between the ground state and the excited states are finite, even if the 
average value of L vanishes in the ground state, the orbital moment 
is also capable of contributing to the magnetic moment in the second 
order perturbation, if one takes the LS coupling and the Zeeman energy 
as perturbing Hamiltonians. The off-diagonal elements of L will be 
anisotropic because of the anisotropy of the energy states split by the 
crystalline field. Therefore the magnetic moment possessed by a magnetic 
ion, and consequently the g-value, will be anisotropic. 

Pryce (1950 a) has given a perturbational treatment of the above 
problem using the ingenious method of projection operators. Let us 
denote the energy levels split by the crystalline field as Hy and #; and 
take the perturbing Hamiltonian as 


V=A (L.S)4e,(L-H}. > See 


where A is the coefficient of the LS coupling and H represents the external 
field. Since the unperturbed states are independent of the spin state 
in the present case, we are able to treat spin S as a parameter in the 
perturbational calculation. Accordingly, we can immediately get the 
second order perturbation energy as follows : 


1 ; : , 
E,=— * Bah, [2 {A | L,, |) S, +p (| L, |0) 3], 
Vs t ie 


where indices , and v represent Cartesian coordinates x, y and z. This 
can be written 


EL,=—2dpy 2A. S,1,— Pe Ae S,S,—pe P94 ie fH eee) 
where * ca a 


0|/ LZ, /2)(|L,]0 

A = 2S el {te eels ae 
The first term in (9.2) represents the contribution of the orbital magnetic 
moment and the second term the anisotropy energy which arises from the 
combination of the crystalline field and LS coupling. The third term is 
spin-independent and gives rise to a temperature independent susceptibility 
which is generally anisotropic. If we add the unperturbed energy 
Ey +-2u, SH to (9.2), the effective Zeeman energy becomes 


2p 2 (Sy, AA) ge ee ee 


wy 
which means that the effective g-value is generally a tensor of the form 


2(1—AA), where 1 and A represent the matrices (6,,) and (4,,,), and is 
real and symmetric because we can choose real unperturbed wave functions. 
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If we take the principal axes of this tensor as 2-, y- and z-axes, the 
Zeeman energy and the anisotropy energy take, respectively, the forms 
LUZ (Geer. J yoy, 1 98H.) Con, aires (9.5) 
and 
—?MAS 2—A,S,2—V2AS,2. se Atak oe COG) 
Thus, the anisotropy constant per ion arising from the crystalline field 
is of the order of times (gy—2). If the crystalline field has cubic sym- 
metry, the principal values of A become equal to each other and (9.5) 
and (9.6) become isotropic. Even then, anisotropy energy arises from a 
higher order perturbation as 


A (S,4--S,'+S,"). ig Say Caen cya) 
It should be noted that when the magnitude of S is equal to 4, (9.6) and 


(9.7) simply reduce to constants; in this case the anisotropy energy does 
not arise from the anisotropy of the crystalline field. 


9.2. Magnetic Dipolar Interaction and the Anisotropic Exchange 
Interaction 
Next, we shall turn to the considerations on the anisotropic interactions 
between two spins. According to Kramers (1934) and van Vleck (1937 b), 
the interaction between spins / and m can be written 
Bs ge sae ee rs hl toe) (9.8) 
where C,,,” is a real constant. As van Vleck has shown, this expression 
is the most general one when S is 4, while for S larger than } higher order 
interactions such as an interaction of the quadrupole type must be taken 
into account. This form of spin-spin interaction involves isotropic 
exchange interaction. ‘To remove the isotropic part, the condition 
DHE Wp) 
is necessary. a 
The magnetic dipolar interaction between magnetic moments M, and 


M i-4 
ee Me) elie ely, (DAG, «Figs 


m? 


can be expressed in the form of (9.8), noticing that 
M/=81 : Ss; and M,,=-38m : Sin 
by (9.4). That is, we get 


Gin Oe tt. lee (2), u Sm, es Ben Sa a : Fin) (Ban ‘die Pain s (9.9.a) 
where r,,, is the radius vector joining the /th to the mth ion and g,_,, and 
Sm,» are the vectors whose components are (8; 42> 81,uv 81, yz) and. (Sn. ii, 
Em,ry Bm,rz)» respectively. The relation (9.9 a) is equivalent to that derived 
by Opechowski (1948). For the case of an isotropic g-value, it naturally 
reduces to 

CA 9 eT On 3LimYim)> Seine ae (9.9 b) 
where 2, », and v,, are the direction cosines of r,,, with respect to the p- 
and y-axes. 
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The anisotropic exchange interaction comes from the LS coupling and 
from an isotropic exchange interaction through the excited states and 
can be derived by a perturbational calculation. We shall take the system 
of two ions J and m and consider the energy states split by the crystalline 
field as the unperturbed states. The perturbing Hamiltonian of this 
two ion system will be given by 

V=A(L,. §,)+A(L,-S,,)+Vex. - - - = (9-10) 
where the first and second terms are the LS couplings of ions / and m and 
V.x means the isotropic exchange interaction. For simplicity, let us 
assume that V,, has no off-diagonal element but has only diagonal 
elements in the excited states, which can be expressed as 


(Va) diagonal element =—2J5, (S, Y Sn)» . oy oi laeee (9.11) 


where J, is the exchange integral in that excited state. The perturbing 
Hamiltonian (9.10) has only off-diagonal elements connecting the ground 
state to those excited states of the two-ion system in which either of the 
two ions is excited to its higher level #; and the remaining one is left in 
its ground state. In the present case, the unperturbed states are inde- 
pendent of spin orientations. Therefore, we may treat the spin operators 
in (9.10) and (9.11) as if they were mere parameters. 

As we observe from the expressions for the perturbing Hamiltonian, 
(9.10) and (9.11), the anisotropic spin coupling between spins S, and S,,, 
occurs only in the third order. From (8.3) and (9.10), it can be seen that 
the most important terms included in the expression of the third order 
perturbation energy #, can be written 


Fy 301 Mba 8) LE Vox 106 Ar S19) 
(L,;—£, = 
4+ym (0| A(L,, - S,,) | 2) [Vex [AACE SLY PO) 
‘ (4 ni — Emo)” ; 


where Hy, H,; and E,,9, H,,; represent the energies of the ground state 
and of the excited state of ions / and m, respectively, and 2 and 2” are 


taken over the excited states of ions 7 and m. . 
By use of (9.11), this expression can be written as 
—2r2 2) Tees (Sas Sin)Spy— 20? 2 LD Sy (Si °Sm)Smy - (9.12) 
yv MY 
where 


j= 3! RRL (0 | Lu | (a | Li, 0) : 
i ales 
(0 [Ln | )G | Ems | 9) 
(Eni — Emo) 
J,’ means the exchange integral between ion J in its ith excited level and 
ion m in its eee state, and J,” has a similar meaning with / and m 


interchanged. J”,,' is a real symmetric tensor, since we can choose the 
unperturbed wave “funotions to be real. 


(9.13) 
My, ye " J; "ihe 


mi 
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For S=3, the following relations hold : 
as Sm)S12= 28 1 gSima—4E(S; + Sn), 
S,(S,. Se =F Smat PSicsmy— YS)S 
S,,(S; S,,)8 le = 48 meas A) eS P Peaed C 1/8 )S 
(9.12) can be simplified to ae 
—2 2D ae Bus = Pit WSSnv— 3 {P,,—46,, we Py }S,Smy (9-14) 


me? 
etc. 


me? 


if we exclude the ate de term. Thoretore we arn the following result: 
for S==- : 
Cy = {48,, 2 Pyf—L,, } AHS, TD yh—Tyy"}. «  (9-15) 
A 


This result is equivalent to that originally derived by van Vleck (1951). 

As we see from (9.3) and (9.13), I’, is of the same order of magnitude 
as J/,,°, since the off-diagonal elements of L are of the order of unity, 
so that J(g—2)? gives the order of magnitude of the anisotropy energy per 
ion arising from the anisotropic exchange interaction. This fact was 
pointed out by Kittel (see van Vleck 1951). Thus we see that the aniso- 
tropy of the g-value is related to the anisotropy energy arising both from 
the crystalline field and from the anisotropic exchange interaction. 


9.3. Anisotropy Constant below the Néel Temperature on the Basis of the 
Molecular Field Approximation 


We have considered above three origins of the anisotropy in anti- 
ferromagnetic substances. These give the anisotropy energy of the spin 
system in the antiferromagnetic state as well as the anisotropy of the 
susceptibility above the Néel temperature. 

In the antiferromagnetic state, the spin lattice is divided into two 
sublattices, one being occupied by plus spins and the other by minus 
spins. The direction of the common axis of these plus and minus spins 
with respect to the crystal axes cannot be determined by the isotropic 
exchange interaction. It is the anisotropic energy discussed above which 
determines this direction. 

We denote the direction cosines of the plus and minus spins, respec- 
tively, as («;, B,, y) and («, 8, y.). We have a,—=—« =a, 
B.=—fh_=B and y,=—y_=y when there is no external field, but 
these relations still hold approximately when the external field is present 
so long as the field is weak compared with the exchange field, which we 
shall assume. Now let us focus our attention on one spin only and employ 
the molecular field approximation. This spin is subjected to the average 
exchange field arising from the surrounding spins and is quantized in 
that direction. In the statistical average, the spin has a finite component 
in this direction but no component perpendicular to it. We introduce a 
new coordinate system €, 7, ¢ where the ¢-axis coincides with the direction 
of the exchange field. A part of the anisotropic energy given by (9.6), 
which can be written in the form 

(A+B) S2—AS,?—BS,?, paeieewe t os fart 9.16) 
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“= 


when the isotropic part is ignored, then becomes 
{(A+B)n,2—AlL?2—Bm?}82+{(A + B)n,n,—Al,l,—Bmym3} (S88, +-S,S:) 
4{(A+B)n2—Al2—Bm2}8,2-+{(A +B)nyy—Alya—Bm B}\(S8-+-8,8;) 
4-{(A-4+B)y?— Ao? — BB }S2+{(A+B)nyy—Aal,—BBms}(8;8,+5,2), 
where (1,, 1, 21) and (Jy, ma, N2) mean respectively the direction cosines of 


the é- and 7-axes with respect to the x, y, 2-coordinate system. In 
calculating the statistical average of this expression, we can put 


(S2 = Ben y=4F(S(S+ ij iS 7 >); 
(S,S-)= (SS, )=9, ete., 


the symbol ¢) representing the statistical average. Thus we get for the 
average anisotropic energy of one spin 


1{(A+B)y?—Aa?— BB{3(S,2)—S(S+1)}, . . . (9.17) 


where (S.) is given by 


s s 
x mexp (mgpy,/kT)/ 2X exp (mgupzH,/kT), 
m=—S m=—S 
H, being the magnitude of the exchange field. The sum of (9.17) over all 
the spins is the free energy arising from the orthorhombic part of the 
crystalline electric field. «, 8, y refer to the principal axes x, y, z of this 
field and A and B are given by (9.16). 

From (9.17) we see that the temperature dependence of the anisotropy 
constant in this case is given by the factor of 3¢S.2)—S(S+1), which 
becomes S(2S—1) at the absolute zero of temperature and zero at the 
Néel temperature. It vanishes identically for S=4 because S,? is always 
equal to } for this case. The anisotropy energy coming from the cubic 
part can also be derived in a similar way. 

Next we shall turn to the derivation of the anisotropy energy arising 
from the anisotropic interaction between spins, likewise on the basis of the 
molecular field approximation. This part of the anisotropy energy is the 
only non-vanishing part for S=4. 

The energy of this interaction for the /-th spin on the plus sublattice 
is given by the sum of (9.8) over m, namely, 

ZL C,,8,78 


m wy? 


Re heh yg CS 


lu 
If we replace S,,, by its average value, this becomes 


Li. Gay Syst k St) yo 6,55 a CSa ye eee 
uy 


ag 


where 

y= ZC, and at a SER yee gh 
The summations 2',,t and 2’, are taken overall the lattice points of the plus 
and minus sublattices, respectively. It can be seen that this energy is 


EEE 


OO EE aaaEaEOEOOeOooo 
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equivalent to the magnetic energy of the plus spin subjected to an 
anisotropic molecular field whose pth component, H,,*, is obtained by 


solvin 
3 HB Guy ae Ayy< St De 2 Cy S— a Sa e Y (9.21 a) 
vad “ i“ 


If the plus sublattice is equivalent to the minus sublattice, the anisotropic 
field for the spin of the minus sublattice is given similarly by 


(eg = 2 kS*) a(S)... (9.218) 
i lu 7 


The anisotropy energy per spin pair arising from the interactions between 
Spins is given by half the sum of the average values of (9.19) and the 
corresponding expression for the negative spin, namely, 


FX alg (St he yes ‘Se eS aco Co Ko” oy - (9.22) 

Noticing that the direction cosines of (S*) are (-La, =+-8, ty), this 
expression can be written as 

<S yt (On = Calon (Qa Cy) (Cea Ca)y 

+ 2(4 py —C py) OP ++ 2(Ay,—Cye) BY +2(Az2—Co) ya}, Pies (9.23) 

where we have neglected the square of the difference between ($+) and 

{S~) and put (S)=3((S+)—(S~)). <S) can be considered to be equal 

to the values of ¢S+) and—{S~) for the case of no external field. The 

anisotropy constant in this case is thus proportional to the square of the 


spontaneous magnetization of each sublattice and shows therefore the 
same temperature dependence. 


§ 10. Anisotropy ENERGY IN CuCl, . 2H,O 

The anisotropy constants of antiferromagnetic substances can be 
obtained from the measurements of antiferromagnetic resonance absorp- 
tion, the field dependence of the Néel temperature or the value of the 
critical field for the flopping of spins. However, these observations are 
possible only when the external field is stronger than the geometrical 
mean of the anisotropy field and the exchange field. This geometrical 
mean amounts to more than 10° oersteds in ordinary antiferromagnetics 
and therefore the measurements of the anisotropy constants are usually 
impossible. CuCl,.2H,O is the only substance whose magnetic pro- 
perties have been thoroughly investigated and whose anisotropy constants 
have been determined. It has a Néel temperature of 4:3°K and conse- 
quently the exchange field is comparatively low ; in fact, the observed 
critical field is in the neighbourhood of 7000 to 8500 oersteds. 
MnCl, . 4H,O, investigated quite recently at Leiden and at Washington, 
has a Néel temperature of 1-6°K, so that in this substance the determin- 
ation of the anisotropy constants may also be possible. 

Another method of deducing the value of the anisotropy constants, 
though indirect, is to observe the anisotropy of the susceptibility above the 
Néel temperature. This arises also from the three causes mentioned in 
the preceding section and affords certain information about the aniso- 
tropy constants below the Néel temperature. Stout, Griffel and Matarrese 
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have measured the anisotropy of the susceptibility above and below the 
Néel temperature in the rutile type antiferromagnetics, MnF,, FeF, and 
CoF,. From their results we can obtain information about the anisotropy 
energy below the Néel temperature, as we shall see in § 11. 

In this section we shall consider the anisotropy energy of CuCl, . 2H,0. 

The crystal structure of CuCl, .2H,O is orthorhombic, as shown in 
fig. 5(a), 5(b). The magnetic ions are subjected to the crystalline field 
arising from the surrounding charges and electric dipole moments, whose 
configuration for the corner ions is different from that for the face-centre 
ions, the former corresponding to the mirror image of the latter with 
respect to the ab-plane. This crystalline field seems to remove completely 
the directional degeneracy of the D-state of each free Cu®+ ion and 
consequently the orbital moment seems to be completely quenched in the 
ground state. 

For each Cu? the crystalline electrostatic potential, which is. not 
rigorously orthorhombic, can be expressed in ascending powers of x, y and 
2, Le., 


Ag?—(A+B)y?+B2+Caz+ D(at+yt+24)+ HE (a4-+ 6y?2") 
+G(yt+ 622x?) + (Ha?+ Ke2*—3(H+K)y*)uz,. . . . (10.1) 


where plus and minus signs are taken for the corner ions and the face- 
centre ions, respectively. Terms higher than the fourth degree are not 
effective for the splitting of energy levels arising from d-orbitals, in which 
we are interested. From this expression for the crystalline potential we 


see that it does not mix real wave functions (¢,—d¢_,)/(4/27i) and ~ 


(b2—¢_2)/(+/2¢) with other real wave functions (¢,+-¢_»)/+/2, (6, +¢_4)/1/2 
and ¢9, where ¢o, 6.1, 6.2 are the wave functions with L=2, m=0, +1 
and +2, m being the z-component of the orbital angular momentum L. 
Therefore the orbital wave functions of the Cu®* ion in this crystalline 
field can be expressed as follows :— 


Pyo=ApgotArs ($ot+¢$_2)//2+A41 (¢1+¢- iV 2, 
Pi =Bobot+ Bye ($o+_9)//2+B,4 ($1+-¢4_1)/v 2; 


P=Cobot+Cy» (bot 2)/V2+C,, (1+¢1)/V2, - + + (10.2) 
P3=D_, (¢.—¢ »/V2t+D_» (62—$_2)/V 2, 
Y4=L_, ($1—$-3)/V/2t+H_, ($2—$_9)/V/ 2%. 4) 


A,, B;, C,, D; and EF, are real constants to be determined by the coefficients 
of the crystalline potential. We now regard such non-degenerate wave 
functions as unperturbed wave functions and their energies, Hy, £,,... E4, 
as unperturbed energies. Then, introducing LS coupling and exchange 
interactions between ions as perturbing Hamiltonian, we obtain the aniso- 
tropy of the g-value and the coefficients of the anisotropic exchange 
interaction by the procedure developed in the preceding section. The 


anisotropy energy discussed in §9.1 doesnot appear in this case, because the 
Cu*+ ion has a spin of 4. 


EEE ee 
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We shall first discuss the anisotropy of the g-value, assuming that the 
ground state corresponds to Y%, as predicted from the cases of other 
crystals containing Cu?+ ions. The matrix elements of L,,, L, and L, 
poe the unperturbed wave functions (10.2) can easily be calculated as 

ollows : 


(0) Z, | )=(0 | LZ, | 2)=0, ] 
(0[L,|3)=—i(4,.D_,+24,.D_.), (10.3 a) 
(0 | L, | 4)=—2(A,,H_,+2A,,H 4), J 
(0|L.|1)=(0|L,|2)=0, 

(0 [Lz | 3)=i(849D_.+4,_D_1+44,D-»), | | (10.35) 
(0 | Ty | 4)=1(3A,H_,+A,.H +A, 1H»), 

(0 | L, | I= —t(349B,1—A, 2B, 1+4,1B 42), 

( 


0| L, | 2)=—1(34,C,,—A+.C,,+A4,1C0,9), ee LOLS C) 

(0[L,|3)=(0| LZ, | 4)=0. 
Inserting these values of the matrix elements into (9-3), we get 4,,. In 
particular, we have 1,,—A,,—0. Since the g-tensor is 2 (1—AA), as we 
have seen in (9.4), we see that one of the principal axes, 7, of this tensor 
coincides with the 6 or y-axis, and the remaining € and ¢-axes make an 
angle « with the a and c-axes, respectively.. We also see that the angle « 
for the corner ions has the opposite sign to that for the face-centre ions. 
The principal values g, and g, of the g-tensor and the angle « between the 
€-axis and the a-axis are given by 


92=WY cx tGee—l (Gee —IeeV- +49 ao? 31, 0 A ae Eo (10.4) 
92=HY cat Gea th (Ic2—IJec) +49 227} 1, Cn See ott OR at (10.5) 
sin (20) = + 29 ool{ (9 2e—IJee)° +49 are? }*. Gi TPES Oe anh (10.6) 


The above calculations are similar to those which Polder (1942) has done 
in his paper on the anisotropy of the susceptibility of CuK ,(SO,), . 6H,0. 

The Zeeman energy for an anisotropic g-value is represented by np SgH 
and so its eigenvalues for S=4 can be obtained by solving the following 
secular equation : 


Poel ot Gagtly te Gell 2) E dunt el eer sully ied eall o) 
OY yl oT IyytytG ues) } 


yiant (Gat gal Gay gt Gigelte) —$ Hp Ger otGevtt ytGootts) Prag ; 
+6G yc atGuyHytIGuet)} 
(10.7) 
which gives 
e=t4ye{ Gael et+Gavd yt9 ot) + (Gul at Iuy yt Goel)” 
Gest et Seg yirdalta) ee =) ~ « (10.8) 


Now if we write the direction cosines of H as ag, By, yq, (10.8) becomes 
o +4upH{ GaetHt9 avPHt9 aca) + JyatatIuPu +9 yx)" 
a5 (Jea%u+Geubnt Gee¥n)}* ; 


76 T. Nagamiya, K. Yosida and R. Kubo on 


therefore we obtain for the effective g-value 
Jost ={ YautHtG evPat9 xv) + (Jy 2tatIyyPHtGysV a) 
+ (Gex%H+GevBu+9evu)?}- BS te fogs (10.9) 


For the present case, g,,—9,-—0, 80 that the effective g-values for the 
cases where the external field is applied along each of the crystal axes are 
Ja= (Gra? +Iex°)' = (92 cos? a9, sin? «)?, 

W=Ivy ; (10.10) 

I= (ue? +922) = (92 sin? atg?? cos? a} 
These values can be obtained either from experiments on paramagnetic 
resonance absorption or from susceptibility measurements at temperatures 
sufficiently above the Néel temperature. The former experiments have 
been done by Itoh, Fujimoto and Ibamoto (1951) and the latter by van den 
Handel, Gijsman and Poulis (1952). The values they obtained are cited 
in §6. To obtain these values theoretically, we must know the numerical 
details of the crystalline field, for which we have at present no reliable 
data. 

The ordinary dipolar energy can be obtained easily by use of the above 
experimental g-values, whereas the evaluation of the anisotropic exchange 
energy is rather difficult, because we have to evaluate J”,, given by (9.13) 
which requires detailed knowledge of the crystalline field. For this 
reason, Moriya and Yosida (1953) and independently Keffer (1953) 
estimated the anisotropy constants of the substance under consideration 
on the basis of much simplified assumptions. 

Moriya and Yosida assumed that the é-axis coincides with the line 
joining a Cu?* ion to its nearest Cl- ion. This assumption gives cos? «= 
0:6425. Using this value, and the experimental g-value, they estimated 
the principal g-values to be 

ge=2-111, g,=2:075 and g-=2-339. 2 7 t1GsEs 
They further assumed that the crystalline field has an orthorhombic 
symmetry with axes coinciding with the principal axes of the g-tensor, 
namely, €,7 and ¢. Then the following splitting of the energy level can be 


Serived: WA Gat ball 24 Aube Be 
P,=Be ($2 ¢- 2)/V2+Bobo; Ey, 
Pa=(b2—$_2)/ 2% elias oeae (hGoes 
P= (b1+-¢1)/V/2 setae 3 
P5=(¢1—$¢_3)/V/ 2% 5 Virdigs J 


For this simplified case, (9.3) gives the following g-values : 
9s=2{1—y5(Ag+ /3Ap)?}, 
In=2A1—y,(A.— V/3A9)}, | 
9:=2(1—4y,A,?), J 


(10.13) 


where 
Vi=A/(E,—E,). 5 BRT eat aly hs ek et ad LO 
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Equating (10.11) with (10.13), we obtain the following relations : 
4yq?Ao4=7-174X 10-8, -y,2(A,+.4/3A,)#=3-056 x 10-3, 
¥22(Ay—1/3Ay)*= 1-406 x 10-3. } 


Since the crystalline field is not known, the problem was further simplified 
by putting A,~1 and A,~0, which seems to be reasonable, because it is 
equivalent to putting g-=g, (see (10.11)). Finally, to evaluate the coeffi- 
cient I”, represented by (9.13), all the exchange integrals for one ion in the 
excited state and the other in the ground state were put equal to the 
exchange integrals for both in the ground state. The exchange integrals 
then appeared in a combined form 2*+ J,, —2X'+~ J,_, where plus and 
minus represent the species of the sublattices and the summations are 
extended over the interacting neighbours. This corresponds to the 
coefficient of the Weiss field and so could be estimated from the 
Néel temperature. Having thus estimated the values of I’,,, they 
calculated the anisotropy energy arising from the anisotropic exchange 
interaction, using (9.15), (9.20) and (9.23) and referring to the spin 
superstructure obtained by Poulis and Hardeman (1952 a) from proton 
resonance absorption (fig. 5(a)). The result for the anisotropy energy was 


i pemefon set KY enl-d 9st. 2h, V5 0s) Sart eos (£0.16) 


(10.15) 


with 
=a lO Cm sand) x5—=— 1-7 x 10-? em=* .- ,. (10.17) 
per Cu ion at the absolute zero of temperature. 

To the result above we have to add the anisotropy constants arising 
from the dipolar interactions between spins, which is given by (9.9) and 
(9.23). For these, we must sum the dipolar part of C ee over all the lattice 
points m. The summation can be carried out by use of Ewald—Kornfeld’s 
method, as has been done by Keffer (1952) for the crystal of MnF,, to 
which we shall refer later. The result obtained by Moriya and Yosida 
(1953), assuming an isotropic g-value of 2, was 


Ky—=—0:27X10%cm-* and «,=29:1x10%cem™}. . (10.18) 
However, when the anisotropy of the g-value, obtained experimentally, 
was introduced, these values were altered to 

K,=0-93X10-8em-! and x,=36-4x10-8em-1, . (10.19) 
The total anisotropy constants are given by the sum of (10.19) and (10.17), 


namely, 
K,=5:-43x10-8em-!?_ and «,=34-7x10-8em-1._. (10.20) 


‘The experimental values to be compared with these calculated values are 
K,=5:1X10-8em™ and x,=16-8x10-%cm-. . (10.21) 


(10.21) has been derived as follows. The analysis of antiferromagnetic 
resonance observations described in § 6 gives the values of 24K, and. 
24K,, where A means the molecular field constant arising from the 
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isotropic exchange interaction between the plus and minus spins and K, 
and K, are, respectively, the anisotropy constants per unit volume in 
the ab and ac-planes. Their values are 
2AK,=5:11X107 00%, K,/K,;=3-3. . . . (10.22) 

Combining these values with the value of A estimated from the perpen- 
dicular susceptibility at a temperature sufficiently below the Néel point 
(A-1=y,=1-40 x 10-4), we have (10.21). We might remark here that the 
value of A estimated from the Néel temperature, neglecting I” (see (4.13)), 
is considerably smaller than that evaluated from y,. i 

Comparing (10.20) with (10.21) we see a good agreement between the 
experimental and theoretical values of «, but a rather poor agreement for 
ks, the theoretical value being nearly twice as large as the experimental 
one. This discrepancy may be due to the crude assumptions made in the 
theoretical derivation of the anisotropy constants originating from aniso- 
tropic exchange interaction. But it should be remarked that the value of 
kK, arises mainly from the latter, while for x, the magnetic dipole-dipole 
interaction contributes the greater part. 


§ 11. Antsorropy ENERGY IN MnF,, FeF,, anp CoF, 

Stout and Griffel (1949), Griffel and Stout (1950) and Stout and 
Matarrese (1953) have measured the anisotropic susceptibility above the 
Néel temperature in single crystals of MnF,, FeF, and CoF,, which have 
a structure of the rutile type. From these measurements we can draw 
some information about their anisotropy energies. 

In crystals of the rutile type the metallic ions constitute a body-centred 
tetragonal lattice and each is surrounded nearly octahedrally by six 
negative ions. The crystalline field at each metallic ion is therefore 
orthorhombic, and the principal axes, €, 7 and ¢ of the field at each of the 
corner ions are in the directions [110], [110] and [001], respectively, 
whilst at each of the body-centre ions they are in the directions [110], 
[110] and [001], respectively. The electrostatic potential of such a 
crystalline field can be expressed in powers of €, 7 and ¢ as follows : 


Ag? By?—(A+B)C+D(Et+-9t +08) + FE 692L2) + Gg 622), 
(11.1) 


where we have written the terms up to the fourth power, since only these 
are effective for the splitting of the energy levels of the d-electrons. This 
crystalline field produces the anisotropy of the g-value, as well as the aniso- 
tropic energy, which is expressible as a function of the spin $ of one ion and 
the anisotropic exchange interaction, as mentioned in § 9. 


MnF, 
The anisotropy of MnF, has been studied by Yosida (1951) and Keffer 
(1952). In particular, Keffer concluded that the main part of the aniso- 


tropy energy of this substance is accounted for by the magnetic dipolar 
energy. 
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The free state of Mn?+ is %S(3d)5. Therefore, there is no orbital 
degeneracy. However, van Vleck and Penney (1934) have shown that the 
configurational interaction due to the combined action of the crystalline 
field (11.1) and the spin-orbit coupling gives rise to anisotropy energy of 
the form of (9.6), while Pryce (1950 b) and Abragam and Pryce (1951) 
pointed out that in paramagnetic salts the dipolar interaction between 
electronic spins within the same Mn?* ion is more important than the 
spin-orbit coupling. In either case we can expect in this crystal the 
following type of anisotropic energy : 


DS2-+-D,S8,2—(D:-+D,)S2 for the corner ion, ] 
é ie 
D,S?+DS8,7—(D.+D,)8?2 for the body-centre ion. j ( ) 


The sum of these two expressions is actually important for the anisotropy 
energy and the value of D =3(D;+D,) has been estimated by Keffer 
(1952) by assuming the relation 


D(MnF,) | A(MnF,)+B(MnF,) 
D(Tutton) ~~ A(Tutton)+B(Tutton) ’ 


where A and B mean the coefficients of the quadratic terms of (11.1). He 
used the rough values of A and B calculated by himself for MnF, and those 
calculated by Polder (1942) for the Tutton salt and also the value of D for 
manganese ammonium Tutton salt deduced from paramagnetic resonance 
absorption. The value of D estimated in this way for MnF, was found to 
be 0-012 cm~!, which was much smaller than that necessary to account 
for the experimental value of the anisotropy of the susceptibility. 

The fact that the observed g-value is almost equal to 2 indicates that 
the anisotropic exchange interaction, whose order of magnitude is given by 
J(g—2)?, is also small. 

The anisotropy energy arising from the magnetic dipole-dipole inter- 
action has been calculated by Keffer (1952) using the superstructure of 
spins determined by neutron diffraction study by Erickson and Shull (1951) 
and Erickson (1953). Apparently the corner sites are occupied by plus 
spins and the body-centre sites by minus spins. The internal anisotropic 
field represented by (9.21 a) and (9.21 6) is, therefore, 


Aft =2O, /M ALO," M,*,-. . . » (113) 


@,,' and ®,,” being given by the following lattice sums : 
Di, =—(2/N) din Pim (One — 3 tmY im)» ] 
®,,’ = —(2/N) 2 Vim (Spr 3H mY im): J 


Since the spin lattice is tetragonal in the present case, the non-diagonal 
components of these @’s vanish when referred to the crystal axes. The 
values of the diagonal components calculated by Keffer are as follows : 


Gi O/= 925, 9G, 0,=—4:85, <5, (11-5) 


P.M. SUPPL.—JAN. 1955 G 


(11.4) 


80 T. Nagamiya, K. Yosida and R. Kubo on 


If we introduce (11.3) as a part of the internal field, the magnetic moments 
of the two sublattices above the Néel temperature can be obtained from 
M+=(C/T)(H—AM*—I-M#+ @'M++@6"MtT), . . (11.6) 
where C' means the Curie-Weiss constant (see § 4 for the molecular field 
formalism). From these equations the susceptibilities parallel and per- 
pendicular to the c-axis are obtained, i.e., 
X= O/{1+40(A+l—®,'—G,"")}, | 
Xe=C/{T+4C(A4+I—G,'—,")}. J 
These equations lead to 
Xe—Xa= Xara Pe — Po 1- Py —Pz ); a om Lae 
which is positive. Keffer calculated this difference as a function of 
temperature by using (11.5) and replacing x,x, by the square of the powder 


susceptibility. The latter has been given experimentally by Bizette and 
Tsai (1939). The results are shown in fig. 27 together with experimental 


Fig. 27 
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The anisotropy in the molal susceptibility of MnF, above the Néel temperature. 
ae after Stout and Griffel (1950) and calculation after Keffer 
(1952). 

values for Xe—Xa obtained by Griffel and Stout (1950). This figure shows 

that the major part of the anisotropy can be accounted for in terms of the 

magnetic dipolar energy. The small difference between these two curves 

may be due to the anisotropy energy of (11.2). 

The anisotropy energy arising from this magnetic dipolar interaction 
can be expressed, as seen from (9.23), as 


Hy=—M,{(®,'—G,')—(®,""—@,"")}y?. oat oS as (11.9) 
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The values given in (11.5) show that the coefficient of y? ig negative, so 
that the preferred axis coincides with the c-axis. This is in accord with 
the results of neutron diffraction. Therefore, the susceptibility along this 
direction, x,, corresponds to x, of the van Vleck theory and that perpen- 
dicular to this direction, y,, corresponds to y,. Thus, below the Néel 
temperature, x, becomes smaller than y,. When the temperature is — 
lowered, the experimental curve in fig. 27 deviates from the theoretical 
curve near the Néel temperature and cuts the temperature axis just 
above the Néel temperature. This behaviour must be due to the short 
range order which should exist above the Néel point but which has been 
neglected in the present molecular field approximation. 

At absolute zero, the magnetization of each sublattice, M,, becomes 
equal to the absolute saturation magnetization given by (N/2) gu,S. 
Using the latter value, Keffer estimated the magnitude of the dipolar 
anisotropy field at absolute zero to be 8300 oersteds. Combining this 
with the contributions from other sources then Keffer obtained a total 
anisotropy field of 8800 oersteds. This gives rise to a value of the critical 
field, H,, of 1-0 x 10° oersteds (see also the value on p. 40). 

It is interesting to compare the anisotropy energy of MnF, with that of 
MnO, which is also of the rutile type. The spin superstructure of MnO, © 
studied by Erickson (1952) was referred toin§ 2. For that superstructure, 
the contribution of the magnetic dipolar interaction to the anisotropy 
energy is entirely different from that for MnF,. The results obtained by 
Yosida (1952 c) are as follows : 


ce 05 6 
Xe—Xa=2'86 XeXa \ ee t10, 


H,=21-46 M,? (y?+y3") 


where y, and y; mean the direction cosines to the c-axis of the common 
axes of spins, for the corner lattice and for the body-centre lattice respec- 
tively. In this case the dipolar energy gives a positive value of x,—y, 
above the Néel temperature, while below the Néel temperature it confines 
the spins to the plane perpendicular to the c-axis, so that again y,—y, is 
positive. This result agrees with the behaviour of the susceptibility of the 
single crystal of MnO,, given in Bizette’s article (1951). However, the 
manganese ions in MnO, are Mn*+ and the free Mn** ion is in the state 
4F' so that the anisotropy energy of the type of (9.6), which arises from the 
anisotropy of the crystalline field, would give a considerable part of the 
anisotropy energy of this substance. This energy can be expressed as 


—A S2—2A,S,?2—VAS? for the corner ion, 

_2A,8 P—VA, 8,27—MVAS? for the body-centre ion. 
Then, corresponding to the experimental fact that y,—x, is positive both 
above and below the Néel temperature, it is expected that 4, should be 
larger than 4(A,+-A,) and less than the larger of A, and A,. Thus, below 


the Néel temperature the spins would point in the direction of the principal 
axis that corresponds to the largest value of A in the plane perpendicular to 


G2 


82 T. Nagamiya, K. Yosida and R. Kubo on 


the c-axis and the direction of the corner spins would be perpendicular to 
that of the body-centre spins. 
FeF, 


tropy of the susceptibility of FeF, is similar to that of MnF,, and y,.—xq 
changes sign near the Néel temperature (see fig. 11). This is consistent 
with the results of a neutron diffraction experiment by Erickson (1953) 
who showed that the superstructure of this substance is the same as that 
of MnF,, the preferred axis of spins coinciding with the c-axis. However, 
the difference of y, and y,, above the Néel temperature is greater than that 
for MnF, by a factor of 100. Stout and Matarrese suggested that this large 
anisotropy originates from the crystalline field, because a free Fe? ion is 
in the state °D. 

Niira and Oguchi (1954) analysed the experimental data of Stout and 
Matarrese and confirmed that the magnetic anisotropy of this substance 
can be accounted for by the anisotropy of the crystalline field acting 
on the Fe?* ions. 

The directional degeneracy of the orbital state of a single Fe?* ion is 
completely removed by the orthorhombic crystalline field (11.1). 
Assuming the energy differences between two levels to be large compared 
with kT’, the effective Hamiltonian of the ground state can be obtained 
from (9.4), (9.5) and (9.6), as 
| py (1-AA,)S pH e+ 2p3(1—A,,) SH, + 2prn(1—AA,)S eH 

—MA,S2—7A,S,?— AAS? (for the cornerion). . . (11.1la) 
and 
2p (1—AA,)S,-H.+ 2ug (1—AA,)S,,H,, + 2uz (1—AA,)S-H- 
—WA,S2—MPAS,2—AS, (for the body-centre ion) . (11.116) 
where €, 7 and ¢ denote the principal axes of the crystalline field acting on 


the cornerion. The terms quadratic in S;, S, and S, in these expressions 
can be written in another form, apart from a constant term, as 


DS?2+D,S,?—(D;+D,)82, ee (11.12) 
D;=—(2A;—A,—A,)/3=X(2g:—9,—9)/6, \ (11 13) 
D,=—(2A,—A,—A,)/3=A(29,—9-—9)/6, 


for the corner ion and the same expression with D. exchanged with D 
for the body-centre ion. Both the susceptibilities y, and y, obey the 
Curie-Weiss law fairly well and so we can estimate the effective g-value from 
the tangent of the graph of the reciprocal susceptibility versus temperature. 
The values estimated by Niira and Oguchi are 


g2=245 and {(g?,+9,2)/2}#=2-24. . . . (11.14) 


Since it is impossible to estimate the values of g. and J, separately, without 
knowing the crystalline field, one may put J:~9, Then the value of 
D-+D,, can be estimated to be 


D,;+-D,~ — NGe—G«)/BV6-7 cm—}, 


According to the experiments of Stout and Matarrese (1953), the aniso-— 


ee aI at 
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This value is very large compared to the corresponding value for MnF, and 
therefore we can expect a large anisotropy to arise from the Be tallies 
field. 


The total Hamiltonian of the spin system of this substance is written 
He to¢ =2I Zp, 5 Sp SAL (up (G Sj e+Gy/Sjn +9 SiH ¢) 
J 


+.DS8;2+-D, Sj? sent aear (Hp (G/See +9 Sp, +9 SpcH -) 


+D,Sp2P+DSp?+DS,2}, (De-=—D;—D,) Ce pe gee eho) 
where indices j and k denote the corner and body-centre ions, respectively. 
The susceptibilities x, and y, above the Néel temperature are obtained 
most simply by the method of diagonal sums developed by van Vleck 
(1937 a) and Opechowski (1937). The result of the expansion up to the 
second power of 1/7 is as follows : 


Nuy?S(S+1)g?? 2Jz8(S+1)  3(D;+D,) [28(S+1) al 
i ecoe3k Tan 1-H) var ee ae ie ais 
-_. (11.16) 
So Npp’S(S+1)(92+ 9,7) + 4J2S(S+1)9 Gy 
CS Seay, rr 3kL'(9e2-+9,2) 
3(9°D +9,°D,) at 
kL (9e?+9,") 15 10f |? 
(11.17) 


and the difference of these two expressions can be written 


4NupS(S+1) 2J28(S-+1) (g2—In\* 
Nome No aT |t Weegee Geet) ers re | 
: 3 (28(S+1 1 

+ 5pr Ge — To} ADEA, + HGeD:+9,2D») | 


(11.18) 


Putting (11.14) and g;=g, into this expression, we obtain a large value of 
the anisotropy in the susceptibility, which agrees in order of magnitude 
with the experimental result. The contribution from the anisotropic 
exchange interaction should be considerably smaller than this, since the 
order of its magnitude is given by J(g—2)? while that of the anisotropic 
energy arising from the crystalline field is given by \(y—2) as seen from 
9.2). 
te the Néel temperature the anisotropy energy arising from the 
crystalline field is given by (9.17). The value of D;+D, is about 6.7 cm“, 
as estimated above. Therefore the anisotropy energy of this substance 
should be very large at low temperatures, the anisotropy field at the 
absolute zero of temperature amounting to 10*~10° oe. hana ee 
In CoF,, the difference between the two principal susceptibilities is 
also very large. However, this difference, y,.— x,, shows a somewhat 
different behaviour from that in FeF, ; it is positive at high temperatures 
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but changes its sign at 150°K, that is, at a temperature considerably 
higher than the Néel temperature of 37-7°K at which an anomaly of the 
specific heat is observed (see fig. 12). 

Recently Nakamura and Taketa (1954) explained the magnetic 
behaviour of CoF,, taking into consideration only the lowest orbital level 
which has a four-fold spin degeneracy. This level is split by the crystalline 
anisotropy energy (9.6) into two Kramers doublets. They assumed that 
the energy separation between these two doublets, which is denoted by 
kO, is greater than the exchange energy. (In (9.3) the exchange energy 
was assumed to be much greater than the anisotropy energy.) In the 
high temperature region (7’>@) these two doublets contribute with 
equal weight to the susceptibility, so that it follows from y,>y, that 
A, is larger than }(A;+A,). In the low temperature region (7'’<@) 
the lower doublet contributes the main part of the susceptibility, and, 
according to Nakamura and Taketa, the effective exchange interaction 
in this lower doublet is anisotropic and has a form 2J,S,-’S9-'+ 27 84,8 on’ 
+2J8,-"S5-’, where $,'’ and S,’ are the effective spins of magnitude 
one half. This anisotropic exchange interaction makes y, less than yx, 
in the temperature range 7',<7'< © because J, is smaller than J,—a 
tendency encouraged by the anisotropy of the effective g-value in this 
doublet. Below the Néel temperature, this anisotropic exchange inter- 
action causes an antiferromagnetic alignment along the c-axis, so that 
Xe—Xa remains negative. Furthermore, the theory gives a non-vanishing 
parallel susceptibility even at absolute zero on account of the existence 
of the higher doublet ; this result agrees qualitatively with the experi- 
ments of Stout and Matarrese. 

We have frequently emphasized the necessity of calculating the 
crystalline field in antiferromagnetic crystals, in order to obtain an 
ultimate theoretical conclusion about the anisotropy energy. Numerical 
calculations of this sort for FeF, and CoF, are in progress at Osaka 
University, and we hope that the results will be published before long. 
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Part III. Sratisrican THEORY OF ANTIFERROMAGNETISM 


§ 12. GENERAL REMARKS 


Like ferromagnetism, antiferromagnetism is a typical example of a 
‘cooperative phenomenon’ and has therefore been one of the favourite 
subjects of statistical-mechanical theories. Unfortunately, however, 
difficulties in the way of a rigorous treatment are even greater for anti- 
ferromagnetism than for ferromagnetism, so that none of the existent 
theories is satisfactory. A statistical theory usually starts from a certain 
molecular model without inquiring much into its deeper justification. 
Thus, it is customary, in the statistical treatments of antiferromagnetism, 
to assume the Heisenberg model (see § 4), in which the spins are coupled 
with each other by exchange interactions expressible by the Dirac ‘ vector 
model’ formula. That is, the Hamiltonian of the antiferromagnetic 
system is assumed to be of the form 

Eee thie. oe! pe tas Cae 
jl 
where S, and S, are the vector spin operators associated with the jth and 
th constituent atoms of the crystal. (The sign of J,, is opposite to that 
used in Parts I and IT.) In addition to this exchange interaction, we need 
sometimes to take account of other kinds of energy, for instance anisotropy 
energy and dipolar interactions. 

Although it is not certain if this formula correctly represents the mole- 
cular interactions in any real physical system, such a model is physically 
plausible and seems to have the properties characteristic of antiferro- 
magnetism, so that it is the general hope that the statistical theory of this 
model will explain the observed behaviour of antiferromagnetic substances, 
at least in a semi-quantitative way. This approach is essentially pheno- 
menological, since the interaction constant J is taken to be more or less 
arbitrary, and may be chosen to explain the behaviour of a given real 
system. — 

In the past ten years, great progress has been made in the general theory 
of cooperative phenomena. Most remarkable was the success of the 
rigorous theory of the two-dimensional Ising model (see the recent review 
of Newell and Montroll, 1953). But the theory has not been very successful 
for three-dimensional systems, particularly for those which require 
quantum-mechanical treatment rather than classical or semi-classical 
treatments. The fact that the spins in the Heisenberg model are quantum 
mechanical operators’ necessarily introduces great difficulties in the 
mathematical manipulations, which prevent us from any rigorous treat- 
ment. Thus, the quantum-statistical theory of the Heisenberg model, 
whether ferromagnetic or antiferromagnetic, lags far behind the classical 
statistical theory of cooperative systems. Of course, we may apply to our 
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problem the various approximation methods used for classical systems, but 
the criteria for their validity usually become much more obscure. It is 
still impossible, for instance to give any reliable estimate of the critical 
temperature in terms of the exchange constant. 

In the statistical theory of antiferromagnetism there is another difficulty, 
which concerns the introduction of the long-range order parameter. Both 
from a theoretical and an experimental point of view, the most natural 
definition of this would be the partial magnetization of the crystal. Let 
us take the simplest case, where the whole crystal is divided into two 
interpenetrating sublattices, of which the magnetizations are M+ and M-. 
These are the statistical expectation values of the dynamical quantities, 


M+ =p 3, Sip, My =—ig 2a Su ee 


where j and / refer to atoms on the two different sublattices, and 
HoS (4o=9"p) Means the magnetic moment vector associated with the 
spin S. Generally speaking, there would be no confusion in defining the 
order parameter if the corresponding dynamical quantity were a constant 
of the motion. We could then pick out those eigenstates which belong to 
a given value of the order parameter and construct the partition function 
in this subspace. This is the case for the idealized Heisenberg model of a 
ferromagnet, where the order parameter is the total magnetization, and 
the corresponding operator commutes with the exchange Hamiltonian. 

On the other hand, the partial magnetizations are not constants of the 
motion of the spins. That is, in each of the eigenstates of the antiferro- 
magnetic system, the partial magnetizations must have broadened distri- 
butions. In the course of the natural motion of the system, the antiferro- 
magnetic pattern can never be conserved. Thus, it is impossible to set up 
a statistical distribution of the system among the eigenstates to represent 
3 ensemble with sharply defined values of the order parameters M+ and 

It is, however, possible to calculate the probability distribution of the 
order parameters even if they are not constants of the motion (Kubo 1952). 
The probability density function for M+ and M~ is given by 


P(M+, M-)= — | | Tr {exp (—B#)exp (A*Mop* +A Mop-)} 
(2ar2)? 'I'r exp (—BH) 
xexp(—At+M+—X-M-)dMtdk- www. (12,8) 


where the integral means the inverse of the Laplace integral, and A+ and 
A~ are selector variables similar to those used in the well-known 
Fowler—Darwin formulation of statistical mechanics. Thus, the usual 
definition of the partition function may be generalized to 


Z(M+, M~)=Tr p (M+, M-)=P(M*+, M-) Tr exp (—B#H). (12.4) 
Naturally we have 


Z(8)=Tr exp (—B%)=| | Z(M+, M-)dM*dM-. . (12.5) 


op 


———— 
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The partition function Z(M+, M-) represents the relative probability of 
the realization of given values of the partial magnetizations. In contrast 
to the case where the order parameters are constants of the motion the 
order parameters here will assume various values, not only because of the 
thermal agitation, as expressed by the canonical distribution, but also 
because of the natural motion. Except for this complication, there is now 
no difference between constant and non-constant order parameters. 
If the partial magnetizations are really adequate to describe the physical 
situation, we may suppose that any radical change of the ordering pattern, 
for instance the reversal of all the spins, will seldom be caused by the 
natural motion, and therefore will have hardly any important effect on the 
thermodynamic behaviour of the antiferromagnet. Thus, we may assume 
the free energy of the system to be 


F(M*+, M-)=—kT log Z(Mt+, M-) ee (1 


which is a function of the order parameters. 
We may also define a generalized partition function, 
BATA )— Tr jexp\(— Be) exp (At M..+-+-A-M,,)}. . . (12.7) 
Then we can see at once from eqns. (12.3), (12.4) and (12.6) that 
F( M+, M-)=—kT log E(At+, A~)+kT(AtMt++A-M-)  .SC«w:ssC(12.8) 
and 
0 log & 0 log & 
oAt 0A 
provided that the integral (12.3) is reasonably approximated by means of 
the method of steepest descents. The transformations, (12.6), (12.8) and 
(12.9), are exactly those well known for the Helmholtz and the Gibbs free 
energies, A+kT' and A-kT being the intensity variables corresponding to the 
hypothetical magnetic fields H+ and H~ acting on each of the sublattices. 
It should be remembered that H+ and H~ are hypothetical, because, if they 
were real, we should have included them in the Hamiltonian, so that the 
partition function would have been 
Tr exp {—8 (#—H+tM,,+—H-M,, )} 
rather than 4, eqn. (12.7). 


= Mt, sate Be ate tiers (12.9) 


12.1. The Ising Model 


The complications of the quantum-mechanical theory of antiferro- 
magnetism are avoided if one decides to adopt the Ising model rather than 
the Heisenberg model. It is true that the Ising model, as a physical 
system, is far less satisfactory than the Heisenberg model, because it 
- fails to describe many of the interesting features of antiferromagnetism. 
The Ising model may, however, be considered as an extreme case of the 
Heisenberg model, where an additional field is introduced, that is, an 
infinitely strong anisotropic field with uniaxial symmetry. This 
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simplification has certain advantages, in that we may use more straight- 
forward mathematics to clarify certain statistical features of the system. 
It is important to note the symmetry between ferro- and antiferro- 
magnetic Ising models. If the lattice is divided into two sublattices and 
if the antiferromagnetic coupling exists only between points on different 
sublattices, the statistical properties of the antiferromagnet are the same 
as those of the ferromagnet with the same structure. This is seen at once 
by reversing the signs of the Ising spin variables on one of the sublattices, 
when the Ising spin Hamiltonian remains invariant if the sign of the 
exchange constant is also reversed. If the lattice structure is such that it 
involves closed circuits on the lattice with odd numbers of interacting 
bonds, then this invariant transformation is not allowed and the symmetry 
relation no longer holds. The triangular and hexagonal nets and the 
face-centred cubic structure are of this type, so that simple antiferro- 
magnetic ordering, due to nearest neighbour interactions only, is impossible 
in them. The symmetry relation means symmetry of the energy level 
spectrum with respect to the sign of the energy. This will never be true 
for the Heisenberg model. The eigenvalue spectrum of the exchange 
Hamiltonian, eqn. (12.1), will generally be asymmetric, which means that 


the thermodynamic functions of the ferromagnetic and antiferromagnetic . 


models are not related in any simple way. 

Another way of looking at the Ising model is to notice that its energy 
expression is just the diagonal element of the Heisenberg Hamiltonian, 
with each spin equal to one half, in the representation diagonalizing the 
z-components of the spins. By the general theorem due to Peierls (1938), 
the free energy for the Ising model must always be higher than that for the 
corresponding Heisenberg model. At high temperatures, the free energies 
of the Ising model and the ferromagnetic and antiferromagnetic Heisenberg 
models approach the same asymptote. At the absolute zero, the ferro- 
magnetic lowest energy is the same as that for a system of Ising spins but 
the antiferromagnetic lowest energy lies lower (see §14). This seems to 
suggest that the entropy, at a given temperature, is higher for the ferro- 
magnetic Heisenberg model, and lower for the antiferromagnetic 
Heisenberg model, than for the corresponding Ising model. It is 
unnecessary here to go into a more detailed discussion of the Ising model. 


12.2. Survey of the Approximation Methods 

The crudest approximation in the theory of cooperative phenomena is 
provided by the so-called molecular field theory, which is equivalent to the 
Bragg—Williams method, well known for alloys and regular assemblies 
(see, for instance, Fowler and Guggenheim 1939). When applied to 
antiferromagnetism, this gives the Van Vleck theory, discussed in § 4. 
Our general experience of cooperative phenomena suggests at once two 
ways of making closer approximations, which we may call the high- 
temperature approximation and the low-temperature approximation 
respectively. The moment expansion method of Heisenberg or the 1/7' 
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expansion method developed by Opechowski (1937, 1939) is an analogue of 
the Bethe—Kirkwood expansion method for a regular assembly and belongs 
to the first category. Another typical high-temperature approximation is 
provided by a quantum-mechanical version of the Bethe—Peierls method, 
which was first applied to ferromagnetism by P. R. Weiss (1948) and later 
to antiferromagnetism by Y. Y. Li (1951). The most important method 
for low temperatures is that of spin-waves. This was first invented by 
Bloch (1930) to apply to ferromagnetism and it led to the celebrated T?? 
law for the saturation magnetization. In recent years, the application 
of this method to antiferromagnetism has been considerably developed. 


§13. High-TEMPERATURE APPROXIMATIONS 

The most systematic approximation from the high temperature side is 
the 1/7’ expansion method. Let us illustrate the method for a simple 
model where the spin Hamiltonian is given by an exchange Hamiltonian 
H .. and the order parameters are taken to be the partial magnetizations 
in the z-direction. 

Generally, the equilibrium values of the order parameters, M+ and M-, 
are determined by the minimum condition on F(M+, M—) with respect to 
M+* and M-, so that the equations 

oF oF 
oM+ OM- 
give possible solutions for the ordering phenomena. By eqns. (12.8) and 
(12.9), eqn. (13.1) is equivalent to the conditions 


PRS, tear hr (18.2) 
Thus, we have first to solve eqn. (12.9) for A+ and A-, obtaining the 
expressions A+(M+, M-) and d-(M+, M-). Theneqn. (13.2) has solutions 
M*+(T) and M~(T), which will be the equilibrium values if they actually 
give a minimum of the free energy F(M*, M~). 
Now suppose that we have found an approximation for log & up to the 
second order terms in A* and A-, in the form 
log H=4 (AA+*#+2BAtA-+ Ad *) 
(which must be symmetrical in A+ and A~ since we have assumed that the 
sublattices are equivalent). Then we easily find, by eqn. (12.9), that 
Mt++M-=(A+B)(At+r), 
M+—M-=(A—B)(At—~A-). 
Therefore, the ordinary susceptibility y is 
= A +B, 


ae eas ee ake « (13:1) 


whereas 
ia A185 


might be called the hypothetical susceptibility corresponding to the anti- 
parallel polarization. Spontaneous antiparallel polarization becomes 
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possible when x’ becomes infinitely large. In other words, the critical 
temperature 7’, is to be determined by the condition 
lim 1/y’(7)=0. ea yy Ate RES 
TT ,+0 
Above this temperature, M+—M-=0 gives the absolute minimum of the 
free energy F (M+, M-). 

Now, the principle of the 1/7’ expansion method is to calculate the 
generalized partition function 4, eqn. (12.7), by expanding exp (—B#) in 
powers of 8. The calculation then reduces to the computation of the 
traces of various products of spin operators, which are quite elementary, 
though extremely tedious for higher order terms. Thus, for instance, 
x’(8) can be calculated in the powers of 8. The determination of the 
critical point is necessarily an extrapolation. That is, one assumes that 
1/x’(B) can also be expanded in a power series in 8, which is to be obtained 
from the calculated y’(8). Terminating the series for 1/y’ at the highest 
power of 8 to which y’(f) has been calculated, one can obtain the approxi- 
mation values for 7',. The first order approximation is of course that of 
van Vleck. The approximations up to the fourth were calculated by 
Kubo, Obata and Ohno (1952), whose results are listed in the table for the 
particular case of S=4. These are to be compared with the ferromagnetic 


Table 3. The Antiferromagnetic Transition Points of Different Lattices 
(Each atom has spin 4.) The numerical values are those of kT y/J. 


Number fehs Expansion methods* Bethe-Peierls 
of Ao (Kubo, Obata and Ohno) method 
Lattice nearest SSE SS GE Ge 
neigh- oe al: S88. | Oguchi- 
bours |™°%°Canprox} 22d | 3rd | 4th | Li | Obata 
‘ : 1:00 | 1:00 | None} 0°313 
2 
Linear chain 1-00 i) Toss Nakel Nae 1-387| None| None 
Quadratic an 2-00 | 2-00 | None} 0-439 ; 
layer * | 70 | 9.361! 1-243] None| 2-073] None} 1-582 
! 
‘: 5 . : 3-00 | 3-00 | None} 0-654 
Simpl 3. 2. 2. 
imple cubic 6 300 | 3.403] 9-746 None| 2:710| 2004} 2-618 
Body-centred ar 4:00 | 4-00 | None} 0-613 ‘ 
cubic ° | #00 | 4.495] 3-827| None| 3:278| 218 | 3-636 


* Upper figures are calculated by the expansion in 1/7’, and lower values by 
the expansion in terms of a=4 tanh J/k7/(2—tanh J/kT). 


case. The most noticeable fact is that the 1/7’ expansion gives very poorly 
convergent results for the antiferromagnetic case. In the ferromagnetic 
case, the second approximation was very bad, but it gave reasonable results 
in the third and fourth approximations, at least for three dimensional 
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lattices, as proved by the corrected calculation due to Zehler (Opechowski 
1937, 1939, Zehler 1950). These are consistent with the results of P. R. 
Weiss (1948) based on a generalization of the Bethe—Peierls method. For 
the antiferromagnetic case, however, the third approximation gave no 
transition for any of the lattices considered, while the fourth approximation 
gave transition points which are considerably different from those given 
by Li’s calculation, which will be discussed below. 

In principle, the 1/7’ expansion method can also be applied to the 
approximate calculation of the partial magnetizations below the critical 
point. To do this, one calculates 4(A+, A~) to higher order terms in A+ and 
d-, obtaining eqns. (13.2) as non-linear equations in M+ and M-. If the 
solution of these equations is proved to be the minimum of F(M+, M-), 
this will give the equilibrium state. The underlying assumption of this 
method is that the free energy F(M+, M—) is essentially analytic in f, in 
spite of the presence of the critical point. This assumption is necessary 
in order to justify the extrapolation which is involved in any of the high 
temperature approximations. One might object to this, pointing out 
that F(M*, M~) cannot be analytic below the critical point. The 
problem is related to the general theory of condensation ; here we only 
mention the fact that the high temperature approximations always 
involve implicitly the assumption that the whole system constitutes a single 
phase. 

13.1. The Bethe—Peierls Method 


Yin-Yuan Li applied the Bethe—Peierls-Weiss method to antiferro- 
magnetism. The basic assumption of the method is that a cluster of n-+-1 
atoms, i.e., a central atom and its neighbours, is statistically represented 
by the statistical operator, 


n 


Pcluster— &XP {(—2J 2 SS; HS oo —A yo - Ret c (13.4) 
j=1 j 


where H, is an external field and H, the internal field representing the 
effect of outer atoms, which is to be determined in a self-consistent 
manner. Taking two kinds of cluster, one being centred on an atom of the 
sublattice a, and the other on the sublattice b, and assuming different 
internal fields H,, and H , in place of H,, the statistical average value of 
the spin on a given type of atom is calculated in two ways. Equating 
these expressions, one obtains two equations to be solved for H, and H ,. 
For vanishing external field, one of the possible solutions is identically 
zero, which corresponds to the disordered state, but another non-trivial 
solution appears below the critical point. The values of the transition 
points calculated by Li are also listed in table 3, and they seem quite 
reasonable. In particular, Li succeeded in showing that the one- and two- 
dimensional lattices will not be ordered. For simple cubic and body- 
centred structures, the critical temperatures are somewhat higher than 
those for the corresponding ferromagnetic lattices. This may possibly be 
connected with the larger entropy of the ferromagnetic system. 
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13.2. Other Approximations at High Temperatures 


Several modifications of the two standard methods discussed above 
have been tried by several investigators. Kubo, Obata and Ohno (1951, 
1952) tried an expansion method using the parameter 


a=4 tanh J/kT {2—tanh J/kT} 


instead of J/k7’. This is analogous to the expansion in power series of 
tanh J/k7’ which is known to converge well for the Ising spins. Although 
this expansion is good for ferromagnetism, it failed to improve the 
convergence of the series for the antiferromagnetic case. Oguchi and 
Obata (1952) proposed a modification of the Bethe—Peierls—Weiss method 
taking the statistical operator 


Potuster=eXP (—2F 2 S)S,/k7') exp (H; 2 S52/kT) lS Ay 
j=1 j=1 


instead of (13.4). There is no obvious reason for preferring one of these 
alternatives to the other. In the formulation of Weiss and Li, H,; was 


interpreted as a local field, whereas in (13.5) it is a sort of selector variable . 


similar to A+ and A~ in eqn. (12.7). After reducing the complete density 
matrix of the whole system, we shall not get either formula exactly. 
Thus, one may, if one likes, use (13.5) which has the advantage that the 
calculation can be carried out without using any perturbation calculation, 
which was necessary in the former formulation. Oguchi and Obata 
obtained somewhat higher critical points than those given by Li. But 
the failure of the method for ferromagnetism makes one uneasy. 
Nakamura (1953) examined another modification of the Bethe—Peierls 
method. He introduced some simplifications and reduced the problem 
to that of clusters with two spins. The results were nearly the same as 
those of Li’s calculation. He was particularly interested in the tempera- 
ture dependence of the perpendicular susceptibility below the transition 
point. Unfortunately, there is evidence that this, and similar calculations, 
are unreliable at low temperatures. Even in ferromagnetism, it has been 
shown that the Bethe-Weiss method leads to an anti-Curie point 
(Anderson 1950 c); quite ‘ unphysical’ behaviour is also predicted for 
the antiferromagnetic Ising model, if the Bethe model is used (Ziman 
1951). The reason seems to be that the model is quite unadapted to 
long range ordering, since it takes no account of the topology of the 
lattice, and really only treats a dozen or so spins coupled by the fictitious 
internal fields. The appearance of a new solution to the equations as 
the temperature is lowered is no argument that a new phase must occur 
there, and any method of checking, say, the stability of the new phase, must 
inevitably lead to consideration of a much more complicated model of 
the lattice as a whole. It does seem, however, that one may have reason- 
able confidence in the method when it treats short range order, where 
only a few spins are interacting. Li showed that the ordinary suscepti- 
bility calculated by his method is consistent with that given by the 1/7’ 
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expansion (we cannot check the singularity because this occurs in x’(B) 
not in x(8)) and it is plausible to extrapolate this curve right down to 
whatever transition point one might discover by other means. It is 
worth noting that Li’s value at 7’, is only 0-871 x, where y, is the constant 
value of the perpendicular susceptibility at low temperatures. The 
temperature should thus fall further, and rise further before a transition 
is really possible, if the perpendicular susceptibility is nearly constant 
in the ordered state. 

Nevertheless, the theoretical situation is quite unsatisfactory and a 
statistical treatment giving a reliable description of antiferromagnetism 
around the transition point would be very welcome. Thus, there are 
some interesting experimental facts for which the short range order 
(about which the molecular field theory can tell us nothing) plays an 
essential role. For instance, the entropy increase in antiferromagnets 
above the transition point is often found to be considerable fraction of 
the total entropy increase. It is not clear at present how much of the 
residual entropy is theoretically expected at the transition point, although 
it seems that the statistical theory would predict a rather small amount, 
not more than 10% of the residual entropy, if the structure and 
the interaction are simple. The large residual entropy is possibly to 
be connected with the complicated structure of the antiferromagnetic 
ordering which is found to be present in many antiferromagnetic sub- 
stances. For instance, Stout suggested that in an MnF, crystal the coupling 
of spins among one-dimensional arrays in the direction of the c-axis is 
much stronger than that between neighbouring arrays. This means that 
the antiferromagnetic ordering is essentially created by the weaker 
coupling and that the stronger coupling makes the short range order 
persist in each of the one-dimensional chains at high temperatures above 
the Néel point. A simple calculation shows that this model gives a 
reasonable explanation of the residual entropy, though it leaves a doubt 
about the absolute magnitude of the susceptibility which is also much 
influenced by the short range order (Stout and Kubo, unpublished). 
The same model was applied by Oguchito the case of CuCl, . 2H,0 crystal 
(Oguchi, unpublished). 


§ 14, Low-TEMPERATURE APPROXIMATION AND THE 
ANTIFERROMAGNETIC GROUND STATES 


At high temperatures, the important quantum states are those which 
belong to the group of levels contributing the largest amount of entropy, 
which means that these quantum states can be regarded, in a sense, 
as being highly degenerate. This allows us to treat the constituent units 
as nearly independent. The molecular field theories, both the crudest 
and the most refined, are methods of extrapolation from the high tempera- 
ture regions, where they are primarily justified, to lower temperatures. 
On the other hand, the statistical properties at low temperatures are, of 
course, determined by the quantum states at the bottom of the energy 
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spectrum, so that the low-temperature theory always requires more 
detailed analysis of the quantum states of the whole system. Thus we 
meet here essentially a many-body problem, which can be treated usually 
only when some way has been found of reducing the problem to a separable 
system. One of the best examples of a successful low-temperature 
theory is the spin-wave theory of Bloch (1930) for ferromagnetism. 

Since antiferromagnetism is just the opposite of ferromagnetism, the 
signs of the exchange integrals being reversed, and since this type of 
interaction is actually more common than the ferromagnetic interaction, 
the quantum-mechanical treatment of an antiferromagnetic system was 
taken up as early as 1930 by Slater. In 1931, Bethe gave a detailed 
mathematical theory of the spin-wave in a linear chain, in which he 
succeeded in treating states involving any number of reversed spins. 
This general theory gave, essentially, a rigorous solution for the ground 
state of an antiferromagnetic one-dimensional array. Some errors in 
Bethe’s calculation were later corrected by Hulthén (1938), who solved 
the problem explicitly and gave an exact value for the energy of the 
antiferromagnetic ground state. Hulthén (1936) was also the first to 
apply to antiferromagnetism the semi-classical method of spin-waves 
due to Kramers and Heller (1934). It was, however, not until recent 
years that the importance of these investigations and the usefulness of 
the spin-wave method for antiferromagnetism became fully recognized. 
In 1952, Anderson re-examined the Kramers—Heller—Hulthén theory and 
showed that the spin-wave method can be used to give a reasonable 
approximation to the ground state of an antiferromagnetic array. Since 
then, several investigators have applied the method to the thermodynamic 
and dynamic properties (such as micro-wave resonance) of antiferro- 
magnetic substances. Although the theory is still not very satisfactory, 
the spin-wave method is, at present, the only reasonable treatment feasible 
for low temperatures. 


14.1. The Bethe-Hulthén Theory for the Antiferromagnetic Linear Array 
Before going into the approximate theory of spin-waves, we shall 
briefly describe here the Bethe—Hulthén theory of an antiferromagnetic 
linear array, which is the only case ever treated in a rigorous way. Let 
us take the Hamiltonian 
HE = 2S Xj 58544 


and assume the spins to be of magnitude one half. The eigenfunctions ‘ 


of #,, can be expanded in the form 


P=2a,, ad de eetaee NQ++e Ny > . id (14.1) 
where ¥,,, ng +++ n, represents the spin function of N atoms, of which r, 
that is, the atoms n, <n, <...< n,, have plus spins, the rest of the 


atoms being orientated in the minus direction. Then, the secular equation 
for #., can be written 


2€4ni noses n= = (a 


fy’ Ns oe Ny’ AL Marks np n'y N' dere n'y) . ° (14.2) 
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where 2 < is an eigenvalue of 
Hl T= (HL gg—4NIV|I 
= —2 2', (S; $;,,—4) oe ne rare (ees) 
and the summation on the right-hand side is over all sets 7,1’, n’.. Dak 
which arise from n,, 3, ...7, by interchanging two neighbouring atoms 
with unequal spins. If the set n,, n.,...n, does not involve any 
neighbouring pair, eqn. (14.2) takes the form 


260, NOs +s Np 


ary a UE eee ne Mpa Mgt ++. mae ot aS 


oe (14.4) 
which at once suggests a particular solution of the type, 


exp [¢ (4,7, +h, n.+ ... +k, 7,)]- 
Therefore the general solution can be assumed to be 


any (POS OG ae exp (iQkp, N;+$ oe Pp; Pi)» cl . 0 o (14.5) 
i 


namely a linear combination of those functions obtained by permutation 
of k,, k,...k,. In eqn. (14.5) Xp means a summation over all r!, 
different permutations of k,...k,, and Pj is the number which is carried 
to the jth place by P. In order to conserve the formal eqn. (14.4) for 
those cases where some of n,...7, are neighbours, we impose conditions 
such as 


= = =0. (14.6 
Le ee Natis:: Mp DR eMn Bars a ® Vig os Ry LM yey 6p 0 ( ) 


The phases ¢,, in the expression (14.5) are determined by these conditions. 
In addition, the condition 
(14.7) 


an, Ng. np ong. ++ Rp, n4+N 
is imposed if the chain is assumed to be a ring. Then, the eigenvalue is 
given by P 

e= 2) (1—cos k,) 
j=l 

with the wave numbers determined by the conditions (14.6) and (14.7). 
Bethe and Hulthén showed that the ground state of the antiferromagnet 
is described by the presence of r=N/2 (N even) spin waves, the wave 
numbers being real. Other possible solutions, including those which 
were called the ‘spin complex’ states by Bethe, belong to the higher 
states. When the number of the atoms J is very large, eqns. (14.6) and 
(14.7) can be written 


Ik 
La)=2na-+b| $v, 9) dy 
cot £4(v, y=} {cob F(a) —cot 4&(y)}. 
The eigenvalue 2 ¢ is given by 
1 . 
2 c=N[ (1—cos h(a) dx. 
0 


and 


P.M. SUPPL.—JAN. 1955 H 
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These equations were solved exactly by Hulthén, who obtained the 
remarkably simple result 
2 «=2N log 2. 
Thus, the lowest eigenvalue of %,,. 1s 
(#,,) min= —2NJ log 2+3NJS 
— —0-8863 NJ=—3NJ X1:7726. . . (148) 
Note that the rigorous value is 1-7726 times lower than —2NJ(3)?, 


which is the expectation value of the Hamiltonian when the spins are 
aligned in an antiferromagnetic pattern in the sense of the Ising model. 


14.2. Upper and Lower Bounds for the Ground Level 
The lowest energy of an antiferromagnetic array is always lower than 
the lowest energy of the corresponding Ising model, which is the lowest 
of the diagonal elements of ”,,. in the representation diagonalizing the 
z-components of the spins. Thus we know that 


(4 ox) min < —$N2 X 2582. clues pened 
Here z is the number of the interacting neighbours, and S is the magnitude 
of the spins. We have assumed that the spins are all equivalent and that 
the lattice is divided into two sublattices. Anderson (1951) also pointed 
out that 


1 
(ox) min > —}N2 258? (+35) Saget 


gives a lower bound. This is easily seen from the fact that %,, can be 
divided into 

H x= Hi; (2F S$; 2 S,) 
and that the lowest eigenvalue of each of the summands, 2S; 2, S,, 
is —2JS (zS-+1). From the variation principle (14.10) follows. If we 
write the lowest energy as 


( .<)min==—4N2 2SS? (1+ *) » | ie oe 
then y must lie between the limits 


0<yal, 


The rigorous value of y for the linear chain is 0-7726. 


14.3. Other Treatments of the Antiferromagnetic Chain 

The spin functions ¥,,,,...n, in eqn. (14.1) may be classified 
according to the numbers of plus and minus spins and the number of 
antiparallel pairs in the chain. If one assumes a constant coefficient for 
the spin functions belonging to each group specified by these numbers, 
the eigenvalue problem of %,, is greatly simplified. This method was 
first introduced by Slater (1930), and was refined later by Hulthén (1938), 
who called this the ‘ statistical method’. Hulthén showed that the 
approximation is fairly good, obtaining y=0-632 in a first approximation 


EO EEE EEE 


ce eNO ese 
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and y=0-698 in the second approximation. He also made use of this 
method to calculate the (perpendicular) susceptibility, which turned out 
to be 

Xo= 0-473 Ng?/2S=0-473 Nu 2/8J. . . . . (14.12) 
More recently, however, Ledinegg and Urban (1953) obtained 


Xo=(Nyep/H) exp {—(2cT/upH)*} 
by solving Hulthén’s integral equation for values of energy close to the 
lowest level. Here c is a constant. 

An interesting application of the Slater-Hulthén method was tried by 
Kasteleijn (1952). By a variation calculation, he treated the ground state 
of a linear chain with anisotropic coupling. He showed that, at the 
absolute zero, the antiferromagnetic ordering becomes possible only when 
the anisotropy exceeds a certain value. This is to be expected, since 
the Ising model is ordered at the absolute zero. 

Syozi (1951) examined another method, which consists in writing the 
Hamiltonian in terms of anti-commuting operators, which, after Fourier 
transformation, represent the wave motion of spins. This is similar to 
the spinor method used for the treatment of two-dimensional Ising models 
(see, for instance, Kaufman 1949, Nambu 1950). The value of y obtained 
by Syozi is 0-476, which indicates that the approximation is not very good. 


§ 15. THe Semi-CiassicaL THEORY OF THE SPIN-WAVE IN 
ANTIFERROMAGNETIC LATTICES 
If the spins are coupled by exchange forces, eqn. (12.1), the equation 
of motion for the spin angular momentum is given by 


| 
i S==[S,, UoJ=—S,X2,2Fy$, - - (15.1) 


the right-hand side being the torque acting on the jth spin, from the 
other spins, through the exchange forces. Thus, the motion of the spin 
system is determined by a set of non-linear equations, which hold both 
classically and quantum-mechanically. If there exists a stable equili- 
brium configuration where the torques all vanish, then by linearizing the 
equations we may treat the small oscillations around equilibrium, which 
presumably represent the low energy dynamical states of the system. 
This is the general principle of the spin-wave method, stated in semi- 
classical terms. 

We shall illustrate this with some simple examples. Let us take the 
two-sublattice model of antiferromagnetism and assume only nearest 
neighbour interaction. These assumptions are easily removed or 
generalized, if necessary. The torques on the spins are zero for the 
antiferromagnetic ordering, that is, for $,=S, on the sublattice a and 
S,=S, on the sublattice b. (Of course, $,-+$,=0 holds for antiferro- 
magnetism in the absence of external fields, but this is not necessary for 
the following treatment. In fact, we shall also consider the ferrimagnetic 


case.) 
H2 
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Iniserting the expressions 
S,=S,+6S, S$,—$,+6S; 


into eqn. (15.1) and retaining only the first order terms, we get 


K88,——Z,2I,, (88;X $,+S,XdS,), en 
48S,=—2Z,2J3 ,, (88, X Sg+S,X4S))- 
This describes wave motion, so that we can put 
—(2/N)t Dei 8S, 
A, =(2/N)! Ze" 8S, RS 
B= (2/N)! Z,e'*! 8S,, 


where N is the total number of atoms and K is a wave number vector. 
By (15.3), eqn. (15.2) is transformed into 

hAg=—2,2J Axx S$ r+2p2J 6%? Be X Sq, | 

ABy=—2Z,2J Bx X Sq+2,2J 64? Ag XS, 
with p representing the vectors connecting any one lattice point to its 
neighbours. Now let us choose the z-axis in the direction of S, and take 
the spin-deviations, 5S, and 5S, as the w- and y-components of the spin 
vectors. Using the variables 

Ag*=Ag,tiAg,, By*=By,tiBy, 


we can write eqn. (15.4) as 


(15.4) 


Axt=iwg(S v4 gt—ygS a2Bx*), 
: (15.5 a) 
Byt=tw9(S q-Bgt—yx8 yA x"), 
Ax =—Iw(S p45 —yxSoBx), 
As 2 - (15.5 b) 
By = —two(SazBr- —y x8 p-Ax ) 
where the abbreviations 
Wy= 22d ,/h, 
Parle 2 es \ (15.6) 
VERA A ped eR P/ 22, 
are used. Equations (15.5) are easily solved as follows. 
15.1. Antiferromagnetic Case 
In this case, we may put 
S,,.=—S8,,=8 
so that the frequencies w, are given by 
w= S(1—yz*)t=w(1—yg*)t . . L (15.7) 
with w, defined by 
Wp= WS. 


This is an important formula for the spin-wave theory of antiferro- 
magnetism. Note that every normal mode is doubly degenerate. For 
small wave numbers, we may put 


yr=1—}K ~-%. K-+0 (K?) 


ee a ee 
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where 
a=2 pp/2J p 
is a constant, being generally a tensor. Therefore the relation 
wo |K| EN eh ek tubes PTE LCS) 


holds at the bottom of the frequency spectrum. This is to be compared. 
with the dispersion law 


woe he Pi a. er sie LOD) 


of the spin-waves in ferromagnetism. 


15.2. Ferrimagnetic Case 


Although ferrimagnetism is not the subject of this review, it is of some 
interest to note the simple results of the spin-wave theory for this case. 
Here the magnitudes of the spins are different for different sublattices, 
which we assume to be S, and S, respectively. Thus we take 


Saz=Sa S,,=—S, 
in eqns. (15.5), which give the following frequencies : 
w= tug! {( Sigh ») es. oR (84 —h 5) © si (15.10) 


The degeneracy is removed here and the normal modes are separated into 
two branches. For small wave numbers, the dispersion law of the lower 
branch is (15.9), just as in ordinary ferromagnetism. This will lead to the 
T°’2 jaw for the saturation magnetization for ferrimagnetism (Kaplan 
1952). 


15.3. Effects of Anisotropy and External Fields 


The above treatment of spin-waves can be generalized in various ways. 
We may also take account of the second neighbour interactions, long range 
forces like dipolar interactions, the crystalline anisotropy and the effect 
of an external field. We can start from some ordered configuration, 
linearizing the equation of motion for small spin-deviations from this 
equilibrium state, and calculate the normal modes. Here we shall not go 
into such a general formulation, but will mention only a simplified case 

where the effect of the anisotropy is approximated by the anisotropy 
field H_, acting on each of the spins as if it were an effective magnetic field 
in the direction of easy magnetization. This is allowable when the 
deviations from the easy axis are small. Then the torques, 6S, x H,? and 
/2,S,XH, are added to the right-hand side of each of eqns. (15.2). By 
H_,? and H, , we mean the effective anisotropy fields for the spins on the 
prulittices a and b, which we assume to be antiparallel to each other: 


For the antiferromagnetic ordering, we may take H,*=—H,’, and 
S,=S,=S. Thus, we find easily 
Wg=[(wyo,)?—wPy ct ©. 2 «2 ~ (15.11) 


instead of (15.7). Here w4=4)H_,/h means the frequency corresponding 
to the anisotropy field H,. This can be again generalized to the case 
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where a constant magnetic field H, is present in the z-direction to give the 


frequencies : 
q w R=WR +z: } P f F (15.12) 


ow” ,=WR—Wy, 

where w,, is the Larmor frequency corresponding to Hy. In particular, for 
long waves with K — 0, we find 

w=[w,(w,+2w,)Jttoy, - - + + - (15.13) 
which are the resonance frequencies when the constant magnetic field is 
applied along the axis of easy magnetization. This is in accord with the 
formula derived by Kittel (see (6.15)). This classical model has been dis- 
cussed by Keffer, Kaplan and Yafet (1953), who have given a kinematical 
description of the spin-wave motion. 


§ 16. THE QuaNTUM-MECHANICAL FORMULATION OF THE 
Spry-Wave THEORY 

In the preceding section, we have used the classical treatment only for 
the calculation of the normal frequencies. It is easy, in this way, to find 
the normal modes and their canonical coordinates and momenta, which 
make it possible to quantize the spin-wave normal modes. But we shall 
here go back to the original Hamiltonian and show that it is equivalent to a 
set of harmonic oscillators in a certain approximation. There are two ways 
of doing this, which are actually equivalent, one being the Kramers—Heller 
method and the other the Holstein—Primakoff formulation. In the first 
method, one observes that the commutation rule 


[Se S]=t8,-~48, bw see 


is, if S, is nearly constant, approximately equivalent to the relation which 
would hold between S,, and S, if they were canonically conjugate variables. 
Also, the expansion 
y 2 oo 
S,=[S(S+1)—S,?—S,7]}}~S— Bee S 

is valid in the same approximation (Anderson 1952). 

The Holstein—Primakoff method consists in writing the spin operators in 
the following forms : 


$ 
st= $,+i8,=@28) (1— 5.) at, 
os Ny ee ee 
£=5,-i8,-e90"(1-2), 
S,=S—n, 


where the operators a and a*, called the creation and destruction operators 


of the ‘spin-deviation ’, or briefly the spin-deviation operators, satisfy 
the commutation rule 


(16.15) © 


a a*—a*a=1, wee ee ee we ss) 
and the spin-deviation » is defined by 


a*a=n. SON es ee Err he: 


' Antiferromagnetism 101 


The equivalence of (16.2) and (16.1) is easily seen if one approximates 
(16.2) by 


St=(28)ta and S-~(28)ta*. 


The definitions (16.2) are consistent with the usual representation of spin 


operators as one can see by writing down the matrix elements, using the 
well-known expressions 


(n+1|a*|n)=(n+1)?, (n—1]|a|n)=n?. 

However, it is to be remembered that the representations (16.2) are not 
quite the same as the original definition of the spin operators. The spin- 
deviation operators operate primarily in a space with an infinite dimension 
rather than in the space with the dimension (2S+-1) of the spin functions. 
This does not matter as long as one is carrying out rigorous calculations, 
because the operators St+ and S~ never connect the proper subspace of 
dimension (2S-+-1) with the rest of the infinite space. But this apparent 
infinity of the dimension of the spin-deviation operators seems to give rise 
to some difficulties when one makes use of certain approximations (Kubo 
1952, 1953). 

In the problem of antiferromagnetism, we introduce two kinds of spin- 
deviation, one for the sublattice a through the definition (16.2) and the 
other for the sublattice b by 


(16.5) 


Thus, defining the spin-deviation operators, a,;*, a,, b;* and b,, for each of 
the spins, the Hamiltonian can be expressed in terms of these operators. 
For the two sublattice models with nearest neighbour interaction, we have 


H =—N2zJIS+22d8 (2X nj+2, 2) 
+258 Digi lf (105); f(%1)0, +4," f (0;)bYf (N,)}—2F 24, HN, . (16.6) 


where f(”) means 
pe 
flny=(1- 5) eee oe. (16.7) 


Now the spin-wave operators are defined by the Fourier components of the 
spin deviation operators 


ag=(2/N) 2, a, ag*=(2/N)' 2, e"% a,*, \ (16.8) 
b,=(2/N)t X, e*™ 6, byt = (2182, of) b,*, 

where the wave number vector K runs over N/2 points in the first zone 

of the reciprocal space of the lattice. 


When these expressions are inserted into eqn. (16.6) and terms of order 
higher than the second power of the spin-deviation operators are omitted, 
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we obtain an expression representing spin-waves. With a little generaliza- 
tion corresponding to the case which leads to eqn. (15.12), we write down 
the spin-wave Hamiltonian in the following form : 
H = —N2JS°—NA yuyS+h(w.+wytwy) 2 Ag” Ox 
By the canonical transformation, 
An=«, Cosh 8y+yx* sinhd,, dg*=a,* cosh 6x+Y, sinh dx, 
by=2,* sinh 8+yx cosh8,z,  bg*=a, sinh 6x+Yx,* cosh dx, 
with 5, defined by 
Pe aK 
WetWs / 
the Hamiltonian (16.9) is easily diagonalized to 
H = —NzJS(S+1)—Nug/S+4)H +2 {Ticog (rag! +4) iad (ng +4)} 
(ng =2,* ex, Ne =YR*Y x) . . . . . (16.10) 
where the normal frequencies are exactly equal to those given by eqn. 
(15.12) (Ziman 1952, Nakamura 1952). 


tanh 26,;=— 


§17. CONSEQUENCES OF THE SPIN-WAVE THEORY 
17.1. The Approximate Ground Levels of Antiferromagnets 
Anderson made use of the spin-wave method for the calculation of the 
lowest energy of an antiferromagnetic array. In eqn. (16.10) the first 
term, —NzJS(S-+1), can be regarded as the energy at the perfect align- 
ment of spins, if the magnitude of a spin at rest is taken as {S(S+1)}* 
instead of S. There is, of course, zero point motion of the spins, due to the 
uncertainty coming from the non-commutability of spin components, so 
that a certain zero point energy has to be added to the naively obtained 
lowest energy. This zero point motion is now organized into spin-waves, 
which give rise to 
Erero=t2K (hwg’+hog”)=Lighwg, . . . . (17.1) 


where wx is given by eqn. (15.11). Table 4 shows the estimates for y in 
eqn. (14.11) calculated by this method. This looks quite satisfactory. In 
particular, the close agreement with the rigorous value, y=0-7726, for the 


Table 4. Numerical Values of y, eqn. (14.11), Calculated by the 
Spin-wave Method 


Simple spin-wave Correction by first 


Lattice Z calculation order perturbation 
Linear chain 2 0-726 0-066 S-! a 
Quadratic layer 4 0-632 0-025 S-1 
Simple cubic 6 0-58 0-014 S-1 
Body-centred 8 0:58 0-010 S-1 


ee) 


Antiferromagnetism 103 


linear chain is remarkable. All the calculated values lie about midway 
between the limits, 0 and 1. This seems to suggest that the spin-wave 
method can be used as a reasonable approximation, although one might 
feel somewhat sceptical about its validity considering the fact that the 
spin-wave method here can give only an approximation even for the ground 
state, while it can start from an exact ground state in the case of ferro- 
magnetism. In fact, the eigenfunction obtained by this method lacks the 
symmetry with regard to the reversal of the spins, which ought to be 
possessed by a rigorous solution (as that of the Bethe—-Hulthén ground 
level). But, as Slater (1930) pointed out, the connection between the 
opposite types of ordering pattern is so indirect that the ground energy can 
be calculated to a good approximation by taking only those configurations 
around one of the ordered patterns. This means that the time required for 
the system to swing from one arrangement to its opposite is extremely 
large. Anderson (1952) estimated this to be of the order of years. 
Physically, this will also mean that the approximate eigenstates, rather 
than the rigorous states, are just what we are interested in, if the antiferro- 
magnetic arrangement is ever observed. 


17.2. The Internal Energy, the Entropy and the Specific Heat 


At low temperatures, the spin-waves of low frequencies are excited. If 
the effect of the anisotropy is neglected, one can easily find the temperature 
dependence of thermodynamic quantities from the approximate dispersion 
law, eqn. (15.8), which is analogous to that for the Debye model of solids. 
Thus, for three-dimensional antiferromagnetic lattices, the internal energy 
E will be proportional to 7%, the entropy to 7°, and the specific heat to 7. 
These predictions are, of course, quite different from those of the molecular 
field theory. These relationships hold only for temperatures in the range 


(22I Sup _,)? <kT <2J8. a ty enemas C4) 


For very low temperatures, the dependence is more rapid than any power 
of T (Kubo 1952). 


17.3. The Sublattice Magnetization 


The antiferromagnetic order is measured by the algebraic difference of 
the partial magnetizations, which may be shown to be 


ike Dae 4 1 (ww, , ” 
Mt—M-=,8N| (1+ 55)— ype — Hg” Tat Sea >} 


(17.3) 


by the spin-wave method. The second term in the bracket is the decrease 
of the magnetizations due to the zero point motion. This over compensates 
the correction 1/28 in the first term, which corresponds to the last term in 
(16.15). In the one-dimensional case, this zero point correction will 
diverge, if the anisotropy field H, tends to zero. This is in agreement with 
the general expectation that the linear chain cannot be ordered even at the 
absolute zero, which is actually proved by Bethe and Hulthén’s calculation 
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and is also guessed from Bloch’s theory of ferromagnetism. In the two- 
and three-dimensional lattices the corrections are convergent in the limit 
of vanishing anisotropy. Thus the sublattice magnetizations are expressed 


in the form 
M==ip,N(S—«) ae (ee 


with «0-197 for the square net, <=0-078 for the NaCl-type, and «=0-075 
for the CsCl-type (Anderson 1952, Kubo 1952). 

The decrease in magnetization with temperature is also found from 
eqn. (17.3). For three-dimensional lattices, the decrease will be propor- 
tional to 7? in the temperature region (17.2) as one sees easily with use of 
(15.8). On the other hand, this temperature dependent correction 
diverges for two-dimensional lattices, which, therefore cannot be ordered 
except at the absolute zero. This is again similar to the consequence of 
the spin-wave theory for ferromagnetism. 


17.4. The Magnetic Susceptibilities 


The parallel susceptibility y,, can be calculated by the spin-wave method 
(Ziman 1952, Kubo 1952, Nakamura 1952, Tessman 1952). It is given by 


2E9" exp (hw, /kT’) 


0 LK (exp (erg) T)— 1 aes 


This is divergent, in the limit of H,—0, for one- and two-dimensional 
cases, but is convergent for the three-dimensional lattices. Thus, in the 
temperature range (17.2), x, is proportional to 7, for example, 

— 2N po" ( (eae ) 2 


AU~ 32S \22d5 
for body-centred structures. At lower temperatures the anisotropy field 
gives an exponential dependence, actually in the form 

xu Tt exp (—A(H,J)! S/kT) 
where A isa constant. At higher temperatures, eqn. (17.5) shows that yy 
becomes proportional to 7’ rather than to 7?. But this is an extrapolation 
beyond the limit of the theory. 

When an external field is applied in a direction perpendicular to the 
spontaneous magnetizations, the magnetic energy, —po(2)S,, +28 .)H » 
is to be considered. In the spin-wave formulation, this energy involves 
only the waves of zero wave-number which means that the effect is only 
to displace the equilibrium point of these oscillators. Therefore, the 
decrease of the total energy is entirely static in this approximation, so that 
the perpendicular susceptibility is obviously independent of temperature. 
Thus we can show at once that 
Ny” 
42J 


bt he (17.6) 


which is the same as that which follows from Néel-van Vleck theory. 


—_— 
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Figure 28 shows the 72-dependence of X1—xXu Observed by Griffel and 
Stout (1950). For other crystals, like FeF,, CoF, and CuCl,2H,O, there 
is evidence that the temperature dependence is faster than T2 (Stout and 
Matarrese 1953, van den Handel, Gijsman and Poulis 1952, Poulis and 


eis 1953). This is perhaps to be attributed to the anisotropy 
effect. 
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Molal magnetic anisotropy y,—y,, of MnF,. Experimental data by Griffel and 
Stout (1950). 


§ 18. FurTHER REMARKS AND SUMMARY 
18.1. Higher Approximations 

The spin-wave theory so far discussed is based on the simplification that 
it neglects the collisions, or the interactions between spin-waves. 
Obviously, this is not valid if the concentration of spin-wave quanta is 
increased. Thus, the approximation of Bloch’s theory of ferromagnetism 
becomes poor when the temperature is raised to about one-tenth of the 
Curie point, which will presumably be true also for the spin-wave theory of 
antiferromagnetism. Furthermore, the approximate nature of the 
antiferromagnetic ground state of the spin-wave method seems to indicate 
that the interaction between spin-waves is even more important for 
antiferromagnetism than for ferromagnetism. 

The effect of this interaction, which is represented by the higher order 
terms in the spin-deviation operators, may be examined by the perturba- 
tion method. The first difficulty here is to find out what these higher 
order terms should be. The discrepancy between the spaces in which the 
spin operators and the spin-deviation operators are operating necessarily 
introduces a certain ambiguity in the definition of the perturbation terms. 
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The simplest, and perhaps most natural way, is to choose the higher order 
terms of the binomial expansion of f(”), eqn. (16.7), 

f(n)=1—n/4S+ .... oe Sly Set Olen 
as the correction terms. One easily sees that this procedure corresponds 
to an expansion in powers of 1/8. That is, the classical treatment of 
spins in the limit of S— oo is the zeroth order approximation, the first 
order approximation being afforded by the spin-wave theory expressed by 
the Hamiltonian (16.9). In the second order approximation, one has to 
consider the correction terms coming from the second term of (18.1) and 
also from the last term of eqn. (16.6). 

Kubo (1952) showed that reasonable approximations are obtained by 
this procedure in the second order. For instance, the corrections to the 
ground energies come out to be less than 10% of the preceding corrections 
(see table 4). The perpendicular susceptibility y, at the absolute zero 
decreases by an amount depending on the magnitudes of the spins but of 
the order of 10%. This is consistent with Hulthén’s result for the one- 
dimensional case, and also with the sign of the change of x, with temperature 
predicted by Hulthén (1936) in the general case. The temperature 
dependence of the perpendicular susceptibility is of some interest. If the 
decrease is real, it means that the susceptibility at the Néel point must 
be smaller than that predicted by the van Vleck theory. The observed 
decrease of y, in MnF, might be correlated with this effect but a recent 
measurement by Bizette (1954) shows that it is almost independent of 
temperature ; the strong temperature dependence of y, observed, for 
instance, in FeF, and CoF, is more likely to be attributed to the anisotropy 
effect (Stout and Matarrese 1953). 

An alternative approach to this problem has been given by Ziman 
(1952, 1953). Instead of expanding the (operator) function f(”) in a power 
series in n, and using perturbation theory, he has replaced the operator 
by its average expectation value, calculated in a ‘ spin-wave state ’. 
The method is not justified rigorously but is, in a sense, not unlike a 
molecular field method, in that a factor is introduced which weakens the 
coupling between the spins as they become more disordered, and this 
factor is then calculated by a self consistency argument. The correction 
to the ground energies is much the same as that found by Kubo, but Ziman 
has concluded that y, should remain constant, even when the interactions 
have been taken into account, and he has shown (unpublished thesis) that 
the difference between this result and the result obtained by Hulthén is 
due to the neglect by the latter of the effect of disorder on the commutation 
properties of the (approximate) spin-wave operators. 

The general effect of the interactions is to crowd the spin-wave levels 
closer and closer as they become more and more excited, so that Xu 
tends to rise fairly steeply above the 72 curve. Ziman has suggested 
that this might give a fair approximation to the Néel temperature, 
defined as the point where y,=y,. It is admittedly a crude extrapola- 
tion, but some of the results are interesting. Thus, it is found that the 
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Néel temperature is proportional to JS(2S+1) and in all cases lies much 
below that calculated by the molecular field and Bethe-Weiss methods. 
This is not surprising, since in this method the ‘ zero point disorder ’ of 
the ground state is included, so that less further disorder, i.e. less tempera- 
ture, is needed for the antiferromagnetism to break down entirely. 

It was also possible to calculate the transition points for the body- 
centred and face-centred structures, taking account of the next nearest 
neighbour interactions. If r represents the ratio of these to the nearest 
neighbour interaction, then it can be shown that for the face-centred lattice, 
when r — 0-5, the critical temperature tends to zero. This value of r 
corresponds to a change from one type of ordering to another, and at exactly 
this special value the lattice cannot decide between the two types, and 
remains disordered. It is concluded that the value r=0-5 predicted by 
Anderson (1950 b—see eqn. (4.22) above) for MnO, on the basis of 
observed values of @/T'y, cannot be correct. The larger values of 1, 
near to unity, are in accord with the theory of interactions, where super 
exchange plays an even more important role than the direct exchange. 
The case of the body-centred structures is similar. 


18.2. Difficulties of the Spin-Wave Method 


The perturbation calculation, however, runs into difficulties over 
convergence, especially when one uses any other method than the 
binomial expansion of f(n), eqn. (18.1). For instance, if f(n) is defined 
by f(n)=(1—n/28S)? for 0 <n < 28 and f(n)=0 for n>>28, then the 
first order perturbation diverges in many cases. The perturbation to 
the ground state energy becomes infinite for a one-dimensional chain, 
although it remains small for three-dimensional lattices. The situation 
becomes worse if one goes to the calculation of the magnetization or the 
susceptibilities. For these, the perturbation will diverge even for 
three-dimensional lattices, in the limit of vanishing anisotropy. This is 
certainly very discouraging to attempts to improve the approximations 
of the spin-wave method along the lines of a perturbation calculation, 
but further investigations are needed. 

This divergence in the perturbation method is closely connected with 
another difficulty of the spin-wave theory, which is revealed by the fact 
that the fluctuations in the magnetizations calculated by this method are 
quite anomalous (Kubo 1952). The temperature dependent parts of the 
fluctuation of the partial magnetizations are divergent for lattices of all 
dimensions in the limit of H,>0. The same is also true for the ferro- 
magnetic case. Closer examination shows that these fluctuations are of 
the order of N“?, D being the dimensions. If the anisotropy is taken into 
account, then they become normal, that is of the order of N. This 
abnormal fluctuation might, at first sight, be correlated with the free 
rotation of the spins in the absence of the anisotropy field. But this 
possibility is eliminated, at least for the calculation of the magnetization, 
by the condition that the quantum number of the spin-wave with infinitely 
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long wave length be zero. It seems more likely that this kind of 
divergence is related to the apparent infinity of the dimension of the space 
for the spin-deviation operators and to the inadequacy of the spin-wave 
method for such calculations. Perhaps it should be noticed here that the 
fluctuation is physically significant, because it is essentially the sum 
total of the short range order which is present in ordered states in excess 
of the long range order. Thus it should be observable, at least in principle, 
by neutron experiments. 

The abnormal fluctuation is also inconsistent with the current picture 
of the spin ordering, which describes the spins as being aligned primarily 
by strong exchange forces, the resultant moving around in the anisotropy 
and external fields. This means that the anisotropy and external fields 
can be thought of as perturbations. The divergence of the spin-wave 
theory is, in fact, avoided if one always takes the anisotropy field into 
account. But, it then means that the importance of the anisotropy field 
is far greater than one has ever thought and the effect of anisotropy 
can never be treated separately. The calculated fluctuation of the 
partial magnetizations is finite at the absolute zero, because of the zero 
point motion. The temperature dependent part is proportional to the 
inverse square root of the anisotropy constant, if a finite anisotropy is 
assumed. But we cannot be confident of this, for the reasons we have 
just discussed. On the other hand, it must be admitted that the actual 
amount of anisotropy energy required to reduce the fluctuations to normal 
is very small, so that for many purposes it does not appear explicitly in 
the final formulae, but has only been used as a ‘ convergence factor ’, 
which has cancelled out during the calculation. It is inconceivable that 
a real physical antiferromagnetic array should not have some anisotropy 
energy available to stabilize it, and the apparent divergence of certain 
quantities as H,— 0 may merely be a formal mathematical description 
of the instability of the system in that unusual case. 

A possible way of avoiding the difficulty was suggested by Kubo (1953), 
who applied the spin-wave theory as a variational method. This method 
is, in principle, analogous to the Hartree-Fock method known for many- 
electron systems. The coupled spin-wave system is now approximated by 
a set of independent spin-wave oscillators which are to be chosen in a 
self-consistent manner to minimize the energy or the free energy of the 
whole system. The collisions of spin-waves result in a certain modification 
of the frequency spectrum, which removes the divergence difficulties in 
the limit H,=0. 


18.3. Summary 


Although the simple spin-wave method has some difficulties, it has been 
successful in giving some physical insight into antiferromagnetism. 
The approximation is generally better for lattices of higher dimensions 
and for larger magnitudes of the spins. Thus, for the actual crystals of 
three-dimensional structure, it is able to give fairly reliable predictions 
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about the free energy, the entropy, the partial magnetizations, and the 
magnetic susceptibilities at low temperatures. The great advantage of 
this method is that it can easily be generalized. For instance, second 
neighbour interactions and dipolar forces can be taken into account 
(Ziman 1952, Tessman 1952). The quantization axes of the spins may 
also be chosen arbitrarily. Thus one can start by assuming some 
directions to be those of the antiferromagnetic ordering and construct the 
spin-wave Hamiltonian to represent small deviations from this order. 
The assumed directions can be determined by the variation principle to 
minimize the total energy or the free energy. This can be applied to the 
general consideration of anisotropy and external fields (Kubo 1953). 

But the mathematical difficulties are great for such refinements of the 
theory as the application to higher temperatures and the calculation of 
higher order quantities such as the fluctuations or the short range order. 
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